Pb设计中出现的各种问题论坛总结_第1页
Pb设计中出现的各种问题论坛总结_第2页
Pb设计中出现的各种问题论坛总结_第3页
Pb设计中出现的各种问题论坛总结_第4页
Pb设计中出现的各种问题论坛总结_第5页
已阅读5页,还剩33页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

1、 请教:PB-design,最陡爬坡实验,CCD等中的响应面分析的相关问题请教:PB-design中的相关问题   在实验设计的过程中,通过阅读文献了解到了很多的,但是同时也积累了很多问题解决不了,希望各位高手们指点一二不胜感激.1 在PB设计中出现的dummy variable的具体含义到底是什么?是对照组吗?如果不是应该遵循什么样的原则去设计呢?2 在PB中实验组数应该是变量数加1,那么在相关的文献中看到15个变量设计为:15+1+4,其中4是dummy variable,但是表格中这4个变量也是有高低水平的变化的,那么设计时是作为15个变量来考虑还是19个呢?3 想问

2、问在最陡爬坡实验中,步长的选择有什么要求吗?纯经验还是有公式的?4 在设计CCD试验那的时候是否要包括全因子实验设计?5 什么是中轴点?   各位高手帮帮忙啊,谢谢了小妹我也正在做这块试验,是培养基优化的刚做完单因素试验正在想下面该怎么设计呢?是PB?还是最陡爬坡?还是两个都要做?反正最后是要做响应面的希望大虾们多多多指点一下下另外,关于,我也在想,是不是必须要做空白项的呢?那么空白项里面的+1,-1是没有具体的水平值的亚,那么在实验中具体该怎么操作呢?谢谢各位不吝指教了我自己是怎么想的:单因子实验只是为了保险使PB实验的结果更加明显而进行的预实验,PB实验本身就是有筛选单

3、因子的功能,如果有把握是可以直接做PB的.而我的实验是先进行单因子,然后是PB,根据PB的实验分析数据做最陡爬坡实验,否则不能很好的确定爬坡的方向以及步长.最陡爬坡实验的步长的选择:根据前面PB实验的结果,做一阶方程的法线,法线方向就是爬坡方向,步长就根据回归系数和规范变量的比值在通过自然变量来换算,算到的结果在综合实际的情况就可以基本确定步长了这个是最近看文献理解到的一些,希望哪位高手指点一下dummy variable 不是对照组。如果你的实验次数为4的整数倍。那么看你考察的因素是多少个。比如你有9个,那么至少要使用 runs=12次 的PB设计。此时,有3个空列。一般的处理方法是,等间距

4、空列。实在不行,你就随机选吧。这并不妨碍你的实验结果。空列只是为了估算误差。2. 这样吧,如果我这样表述,看你是否能够理解。我刚好有15个因素,想做PB,正好 runs=16 符合我的要求。可是,这样一来,估算误差呢?因为没有多余空列了。那么,再往上加列,只有再加4列。刚好19因素,20次实验,其中有4个空列。你认为这多出来的一列(16-15=1),我想你是这样算的吧。那你运行一下软件,看看15因素的PB, runs=16,到底是15列,还是16列?列是因素,行是实验次数。并不是列数=行数。描述:如何确定步长:图片: 图片: 描述:4 在设计CCD试验那的时候是否要包括全因子实验设计?图片:

5、看了这张图,我想答案应该是不言而喻了!描述:CCD各个设计点的作用:图片: 综合最后两个问题.如有异议,欢迎继续发帖,大家共同学习!几篇用响应面优化培养基的文章!0 这几篇都是用响应面优化培养基的文章。响应面现在是比较常用的优化培养基的方法,但是好多人(包括我在内)都不是很清楚。由于最近我要进行培养基的优化,所以找了很多文章来研究。一般的文章都是一带而过,看了也不明白。这几篇过程还算详细,对理解有一定的帮助。希望对大家有所帮助。我在看文章的过程中,我总结出了优化培养基的步骤。拿出来跟大家分享一下,共同学习和交流一下。有错误的话希望大家及时提出。第一步:在众多实验因素中找出主要因素。在这步实验中

6、你可以应用正交试验(因素比较少)和PB(Plackett-Burman)实验。这两种方法都可以找出主要因素。尤其是PB实验,它可以在很多的因素中,用较少的实验筛选出主要因素(一般选取大于90%)。通过PB实验还可以看出各因素的作用效果,即是增加还是减少浓度会使响应值向最优移动。第二步:是主要因素的取值逼近中心点,最陡爬坡实验这步实验不是必须做的,如果你确定你的实验取值已经逼近中心点,那么你可以直接进行第三步的分析。但是你要是不能确定或不相信这些取值那你就要进行最陡爬坡实验。这步实验根据第一步实验进行。为了尽快逼近最优值,增加步长通常取最大。我现在还没找到取最大步长的公式,有谁知道的话麻烦发一下

7、。细履平沙版主,帮忙解决一下,谢谢。这样就能找到中心点。第三步:响应面分析现在常用的有中心复合法和BB法(Box-Behnken)。在这步实验时最好因素不要太多,因素太多直接影响到试验次数,现在经典的一般是三因素。通过这步分析可以的回归方程,进而得到最优培养基。并且还能得到因素相互作用对响应值的影响。在这里我对x,即+1、1的取值还不太清楚,的取值也不是很清晰,希望大家帮忙解答一下。这是我最近的学习结果,希望大家批评和指正。 此贴被细履平沙在2008-06-05 18:55重新编辑 附件: 采用响应面分析法优化吩嗪_1_羧酸的发酵条件.pdf (180 K) 下载次数:4 需要经验值:20 附

8、件: 发酵法产氢培养基的响应面分析优化.pdf (908 K) 下载次数:4 需要经验值:20 附件: 酵母内海藻糖积累条件的优化.pdf (213 K) 下载次数:4 需要经验值:20 附件: 响应面法优化纳豆激酶液体发酵.pdf (217 K) 下载次数:5 需要经验值:20关于培养基的优化试验套路0 正交试验和均匀设计方法进行培养基优化已取得诸多成功的例子。  正交试验适合因子较多而因子水平不多的试验设计,从试验次数上看,是至少为因子数的平方。  均匀设计适合于因子少,而水平多的试验,从试验次数看,至少是因子数的两倍。  两种方法虽然多从拉丁方设计衍生而来,

9、不过效率却更高。  现如今,大多流行响应曲面设计来优化培养基。  首先,我们要从众多培养基成分及影响的环境因素中筛选出具有主效应的因子。这时,通常采用筛选试验。主要有全因子因析设计和PlackettBurman设计。两种筛选试验,各有千秋,但都能以最少的试验次数筛选出主效应因子。其中全因子设计能够表现出因子的三级以上交互作用,而PlackettBurman设计由于是两水平设计,所以交互作用只在二级交互作用。另外还有部分因子因析设计。  筛选到了主效应因子,我们就可以开始进行下一步优化试验。此时,主要有中心复合设计和Box-Behnken设计。  中心组合

10、设计是一种国际上较为常用的响应面法,是一种5水平的实验设计法。采用该法能够在有限的实验次数下,对影响生物过程的因子及其交互作用进行评价,而且还能对各因子进行优化,以获得影响过程的最佳条件。  Box- Behnken设计是另一种国际上较为常用的响应面法,是一种3水平的实验设计法。同样具有响应面法的优点。近年来利用该法进行生物过程优化的文献比用中心组合设计法的明显地少。  通常以上说的响应曲面设计和数据分析,都可以通过一些统计软件来运行,十分简便。  为此,我将本人一直在使用的Mintab绿色软件上传,希望大家能够好好利用,搞好试验设计,节省人力无力,为经济带来腾飞

11、(别忘了日本很大程度上就是因为田口宏一的正交设计表,在二战后迅速崛起!)谢谢西风,看你的回复我很有同感,我当时刚开始做RSM也是摸不到头脑,没办法,只有看文献一步一步的摸索过来的。为了不让帖子沉了,也为了有需要的战友节省时间,心中对试验设计有个大概的框架,我先抛砖引玉,简单讲下我做培养基优化的一点体会。可能有些战友认为,试验设计优化培养基就是找几种培养基成分放在一起做几个组合进行试验就行了,其实这样做往往是得不到最适的结果的,尤其在做响应面分析(RSM)的时候,试验得到的图形可能不是一个完美的曲面图,我见过一些别人发的文章,它的曲面图是一个扭曲的图形,高点根本没有出现在图形中,那试验得到的根本

12、不是最适的培养基组分,当然这种文章一般不会出现在核心期刊上的。这种情况一般都是没有在做RSM前,做“爬坡试验”的结果,“爬坡试验”的目的在于找出RSM设计的中心点,保证结果的准确性。另外一点需要提出的是RSM的一个缺点就是分析的因素不能太多,经典的是三因素的试验设计,如果你的因素比较多,那不妨试试方开泰的均匀设计。要做RSM就要先找出有显著影响的因素,剔除一些对结果影响不显著的因素,所以在“爬坡试验”之前还需要一个PB试验,PB能从众多因素中很简单的找出对试验结果影响最为显著的因素,排除一些影响不显著的因素,降低后续试验的工作量,又能保证结果的准确性。PB和“爬坡试验”对于刚接触试验设计的战友

13、,比较陌生点,加上发表文章中很少,全面的把整个试验设计写出来的,一般的都是只写RSM这部分,也有几篇文章是PB+RSM的,还有先做RSM后做“爬坡试验”的。我个人认为PB-“爬坡试验”-RS,这种试验设计的思路比较符合逻辑,又比较全面。如果做毕业论文的培养基优化部分,可以使内容比较丰富。当然如果是做毕业论文的话,在综合试验设计之前,还有个更基础的工作,就是单因素试验,那种放弃单因素试验的想法也完全错误的,尤其是没有一个基础配方的时候。我做的菌种是我自己从本实验室菌种库筛出来,然后又做了诱变育种,所以是没有合适的培养基的,所以做了很多工作,用单因素试验首先大致确定C、N、P、无机盐、微量元素的种

14、类及水平。大致讲了一下,说的比较散乱,或许不太好明白。希望大家把问题提出来,我们一起针对具体问题具体分析。讲得不错,都是很好的心得。先用PB进行显著性分析是一个很好的建议。另外可以从历史数据中找到合适的数据,然后再用软件进行分析,就可以不用再做实验就可找到影响显著的因素和水平范围。响应面分析,其实用我们常见的软件就可以做出来。首先对各个单因素的试验结果进行归纳,后选择变量的变化区间还变化间隔就可以,常见的软件如DPS就可以了。用这个软件可以完成你所需的所有工作,要是想让你的响应面好看,用MATLAB也可以,三维,四维的都可以,在MATLAB区中有我写的现成的用MATLAB语句实现第四维的表示方

15、法,很好用。另外用SARS软件也是可以的,不过语句能比MATLAB复杂一点点。EXCEL也可以,不过做出的响应面,美观程度差一点点。我有个想法,能不能在初始时用正交实验确定出中心点的范围,然后将数据用SAS中的主成份分析,得到影响最显著的几个因素,然后用RSM做,这样主要影响因素有了,中心点也有了。做起来也好做,不会PB等软件的战友也可以轻松操作。不知这样可不可行?希望大家讨论讨论PB实验设计的数据处理我使用SAS进行试验设计的,11个因素+4个dummy variables,做16次实验。PB实验设计的数据已经出来了,但是怎么分析不了。用AnalysistStatisticsRegressi

16、onLinear分析出来的结果像图里面这样,怎么回事呢?后面几列的Standard Error, t  value 以及p值都没有。我给试验设计的表也列出来吧,请各位高手分析一下,到底是哪方面原因导致的。表中X2,X4,X8,X15这4列是dummy variables,没有赋值,其余11个是考查的因素。y是考查指标(转化率%)上面是图片,顺便给表格放在word里面,以附件传上来,请大家帮分析分析。想知道的几个问题:1.试验设计有没有问题?2.就按照表中设计的试验做了16次实验,没有将每次重复3遍取平均值,是不是可以?3.这些数据该怎么处理,怎么用SAS处理会出现上面的情况,没有标准

17、误差,t值以及p值?谢谢了先附件: Coded Design.rar (7 K) 下载次数:5很想问一下,那个虚拟变量如何设置啊 一般多少个还有最重要的是 实验过程中到底如何进行啊  小妹真的没有概念虚拟变量好像一般设个三个左右实验次数都为四的的倍数,如果你要测定8个因素,实验次数仍为12,那么就可以设三个虚拟变量。8311<12记住k个因素和实验次数之间的关系是  K小于或等于n1具体实验过程中,我个人觉得就是高低水平都为,就是什么都不加,只是增加了实验次数,控制误差的作用。实验过程中控制误差可以通过重复试验达到。而实验误差的计算可以通过增加中心点次数实现。如果不能

18、进行误差分析,那么就谈不上数据分析了。我也是刚刚起步,幸亏有发酵人,大家可以一起交流!一般来说,试验设计的时候k次实验的话,(考查因素+虚拟变量)=k1。比如,次实验,考查个因素,那么剩下的个变量就可以不赋值,仍为。不过这种dummy variables选择国内期刊和国外期刊有所不同:国内的:如版友们所说,随机选择的比较多。外文期刊:1.没有标出dummy variables,比如12次实验,设计表就列出了8个考查因素,而3个dummy variables没有标出;         2.标出了,但是给dummy variables都列在一块,比

19、如就放最后几列。    3.也看到就根本没有设 dummy variables的,比如12次实验就考查11个因素的。初学求教:什么是PB实验和爬坡实验啊什么是PB实验和爬坡实验啊刚开始做实验,什么都不知道,万望各位赐教!pb是plackett-Burman的缩写,是一种筛选试验。就是从众多影响因素中筛选出具有主效应的因子。至于最陡爬坡我也不太清楚。我也想知道。额要点名了,细履平沙斑竹,指点一下吧。谢谢了!斑竹,给大家讲解一下吧,我也特别想知道什么是最陡爬坡,我看过你的帖子,知道你是试验设计的行家。望斑竹不吝赐教。最陡爬坡实验它能最快逼近最大响应面区域,确定中心组合实验的中心

20、点,能保证响应面分析结果的准确有效性.我只知道这么一点吧.还望楼下的人能补充细化一下呢.一般都会作,但是如果实验结果由方差分析发现选取的水平已接近最大响应区域,那么就可以不再进行最陡爬坡实验,最后用中心复合设计及响应面分析来确定主要影响因子的最佳浓度。你可以看看有关资料.一般是先作部分因素实验(FFD),再做最陡爬坡实验.再利用中心组合实验确定中心点,如果采用PB实验,一般不作最陡爬坡实验,直接用中心组合实验确定中心点.具体原理我不记不清了.明天我上传一个这个方面的资料吧.上传一篇文章.可能对你有帮助. 此贴被细履平沙在2008-06-08 15:02重新编辑 附件: 生物过程的优化.pdf

21、(358 K) 下载次数:63 需要经验值:30 怎么又方差分析看所选水平是否已经接近最大响应区域呢?谢谢!PB试验可以用来筛选培养基中不同成分对试验结果的影响,通过PB试验得到影响较大的试验因素;爬坡试验就是从PB试验的中心点开始,根据PB试验的得出的每种因子的系数来确定爬坡方向以及步长,爬坡试验主要目的是逼近主要试验因子的最适添加水平;之后根据爬坡试验的结果就可以做相应面了,也就可以最后确定每种重要因子的添加量了。嘿嘿,我也是初学者,希望这些能对你有帮助请教一下:响应曲面1.414到-1.414中间5个水平的步长如何设计?还需要看回归系数吗?设计的时候凭经验还是有具体的公式?紧急求救0是中

22、间,1如果相当于中间偏大2.0水平,那么1.414就相当于中间偏大2.828水平.我是这么理解的,这个不能靠经验,而是根据程序要求。关于PB设计中的t检验和P检验请问各位大侠,PB设计中的t 检验和P检验的功能是否是一样的?如何通过t 检验来认定显著因素?另外,某因子的显著性大于90%或80%等是如何判断的?感谢赐教P的意思是概率,不是P检验,是根据T值查到的概率和F值查到的概率。楼主可能没有学习过数理统计。进行实验设计,一些数学基本知识还是要具备的。关于最陡爬坡(最速上升)的问题我下了发酵人上的“Path of steepest ascent /descent”。在学习中有些地方不明白,希望

23、大虾能给予帮助。在设计中需要编写一个“最速上升宏”命令,请问这个“最速上升宏”在Minitab15这个版本中有没有,如果有在哪里?(本人愚钝,如果翻译错误望大虾指出以便更正)还有我点击“编辑器(D)启动命令(E)”然后按照minitab技术支持文件上所给出的方法在会话窗口里输入以下命令:ascent y x.1-x.n;sore u.1-u.j;step s;base b;descent;runs r.然后回车出现“* 错误* 未知 Minitab 命令: ASCEMTB > ”我不知道这个命令应该在哪里输入。望高手指点啊!最陡爬坡必须要输入命令才能完成吗?可不可以像pb一样直接用min

24、itab软件现有的程序就可以完成设计和分析?最陡爬坡实验我觉得最陡爬坡实验不需要软件分析的,直接观察数据变化的趋势,找到极点值,附近找到因子参数范围再做下一步的响应面实验和分析就可以了。可以参考PB试验设计外文资料里面dummy variables要么没有列出来,有些就根本没有设(比如8次试验就做7个因素),要么所有的dummy variables就全部都列在一起,只有看到国内的文章才是将dummy variables随机的列出来。在版上也看到说要有dummy variables,并且要随机的选择。那么dummy variables到底要不要,要的话是不列出来,还是所有的就放最后几列呢?我认为

25、dummy variables的位置 不影响分析结果。但是这个参数应该有。有道理,另外,每一组试验只做一次,还是做个重复比如同一组做3次,然后取他们的平均值呢?要是只做一次会不会误差比较大,影响试验结果?重复多次可以减少某些异常点引入的误差,减少噪音,提高模型显著性。有两种方法,一种是在程序中设置3个重复,实验结果都输入实验表,也可以在程序中设1个duplication,然后将数值平均后输入实验表。 从pb到最陡爬坡,如何设计步长,求救啊pb做完了,也做了重复,现在真是不会确定步长真的是很笨,希望高手指点一下求救啊细心看看以前的贴子,有详细的介绍。一个参数(A)根据自己的经验确定,其

26、它参数根据PB试验公式中的系数来确定。具体找贴子吧。看下面的贴子,topmit有介绍。做了响应曲面分析,图形不对,麻烦大家看看哪里出了问题图片: 我是先单因素-pb-最陡爬坡-响应曲面的结果很不如意啊,想听听大家的意见,修改一下,继续向前挺进呵呵大家多提意见啊hope59你好,我想请教一下你爬坡实验怎么设计的,是用minitab设计的吗?关于爬坡步长我看了些文献,关于爬坡步长的说法不一。  有文献这样确定最陡爬坡的,步长公式=e×j ×bj, 其中bj为一次回归的回归系数,j为变化半径这样看来,回归系数越大,步长应该越大啊。  而有的文献又说回归系数越大

27、,步长应该越小。我觉得这种说法也对,因为回归系数大,说明对响应值的影响越大啊,响应值对因素变化非常灵敏,所以步长取值应该小点,以免错过最优点。不知我的理解对不对,请大家指教。求助:可不可以用pb方法进行非培养基优化的发酵实验啊小女子刚接触发酵实验,对很多东西不是太了解,最近查文献资料发现PB方法进行实验优化很方便,但是所查的资料都是关于培养基优化的,很想知道是不是可以用PB进行其他非培养基优化的实验?比如说有很多因素像温度、PH、钙离子浓度等。如果可以用PB法优化,那在进行“确定因素水平”阶段是不是也可以用“最陡爬坡实验”进行非培养基优化的实验呢?希望高手给予解决,小女子不胜感激没有问题,可以

28、用。但是进行试验之前,要考虑好试验的目的,不能只是为了试验而试验,而要考虑试验的结果是不是能够达到某种效果和实验条件能不能实现。优化包括:培养基优化(碳源,氮源,无机盐等),培养条件优化(温度、PH、转速,装液量等)实验设计:根据需要,自由组合,可以同时优化众多因素,不必拘泥于是培养基还是培养条件。响应面设计的三个阶段It is possible to separate an optimization study using RSM into three stages. The first stage is the preliminary work in which the determina

29、tion of the independent parameters and their levels are carried out. The second stage is the selection of the experimental design and the prediction and verification of the model equation. The lastone is obtaining the response surface plot and contour plot of the response as a function of the indepe

30、ndent parameters and determination of optimum points. 三个阶段。1,准备阶段,决定独立变量及其水平。2,实验设计的选择及模型的验证3,根据响应面来确定优化条件。求救:pb试验我做了好几次,总是找不到影响显著的因子我做了好几次pb试验,可是总也找不到影响显著的因子,就是大于95%的两个到三个因子,我的基本都在80%左右,很是郁闷啊备受打击,但是郁闷完了,还是想找找问题,调整浓度继续筛选不知道大家有没有遇到类似的问题?都从哪里入手调节?对了,我看很多文献上都说,高水平一般是低水平的1.25倍,这个很固定吗我都是按照1.25倍做的,不会是这里有问

31、题吧还是要focus到上面的问题,遇到显著性不高的问题,应该如何调整?希望高手多多指教!试试高水平是低水平的1.5倍看看生物过程存在严重的交互作用,而PB实验次数少,因此在交互作用解决上有一定的问题。所以在生物领域先进行单因素考察,然后再进行组合试验是有一定道理的。我感觉也可以先进行PB实验多考察一些因素,然后再减少因素并进一步探索因素间的交互作用。首先:你考察的是几个因素?如果存在交互作用,确实会存在上述情况。第二:你的响应量是什么,该数据是否可靠?  比如,我们做糖的,响应量就是胞外多糖,胞内多糖。测糖的方法是浓硫酸苯酚法。如果实验设计者对该方法掌握不好,就算是同一个样品,测定结

32、果差别甚大,那么,这就不是实验设计分析本身的错误,而是实验设计者的错误!只要大于就行了我是以生物量作为指标,用涂平板检测的我用了七个因子,用minitab设计的12次试验,经过分析,只有一个因子的可性度在90%以上在此之前,我做过比较全面的单因子试验,但是由于单因子试验用的是化合物能清楚的知道碳浓度、氮浓度等详细的参数现在做工业发酵,用的都是豆饼粉之类的,所以在量上还是需要摸索后期要做响应面,所以pb试验做的比较谨慎,怕后期出现更大的错误。我这周重复一下上次的试验,再把倍数扩大到1.5试试看非常感谢大家这么热心的帮助我,真是很感动哦不知道我现在的设计有没有什么地方不合适?希望大家多多批评指正!

33、不知道楼主,为什么7个因子要用12次实验设计呢?minitab上自动生成的啊还有就是每次生成的图表都不一样呢,很是奇怪版主是不是觉得我的设计有问题啊我刚准备瓶子,准备做了,呵呵,看样子是有问题了哦我回去再看看,晚上再做吧,版主加我qq吧 43787903可能需要和您好好讨论一下呢谢谢,我在线等首先:7个因子,做8次实验足够。不知道你设计过程中是否有哪一步没选好。其次:低水平和高水平的倍数,一般在1.5左右,多了少了并没太大关系。只要不是认为的胡乱加大水平值,PB筛选还是比较可靠的。再次:我实验室是不准安装QQ等聊天软件的。所以,只能我回家加你了。关于爬坡实验的疑问!爬坡试验需要软件设计吗?我刚

34、做完pb试验,筛选完重要因子,然后是不是应该以高水平为中心,上下筛选几个浓度(一般也要七八个浓度吧),这时候步长设计有什么要求吗?浓度个数设置有要求吗?需要软件设计吗?版主快来指导一下吧!最爬坡试验最后要达到什么要求才能进一步进行响应曲面试验?爬坡试验需要或者能用minitab设计吗? 我看到的论文上好像只是将得到的数据比较产量高低呢?puzzled and help!你做完了PB从分析的数据上就可以看出每个因子的影响是正还是负!所以设计的时候一定注意不要搞错了!步长的设置也要看你自己实验本身,不能套用!所选择几个比较好,能达到你的实验目的,就是逼近最有的响应区域!越逼近越有利于下一步的实验!

35、回归与相关回归与相关-一、直线回归分析直线回归是用直线回归方程表示两个数量变量间依存关系的统计分析方法,属双变量分析的范畴。1. 直线回归方程的求法 (1)回归方程的概念:直线回归方程的一般形式是Ý(音y hat)=a+bx,其中x为自变量,一般为资料中能精确测定和控制的量,Y为应变量,指在x规定范围内随机变化的量。a为截距,是回归直线与纵轴的交点,b为斜率,意为x每改变一个单位时,Ý的变化量。(2)直线回归方程的求法确定直线回归方程利用的是最小二乘法原理,基本步骤为:1)先求 b,基本公式为b=lxy/lxx=SSxy/SSxx ,其中lxy为X,Y的离均差积和,lxx为

36、X的离均差平方和;2)再求a,根据回归方程 a等于Y的均值减去x均值与b乘积的差值。(3)回归方程的图示:根据回归方程,在坐标轴上任意取相距较远的两点,连接上述两点就可得到回归方程的图示。应注意的是,连出的回归直线不应超过x的实测值范围.2. 回归关系的检验回归关系的检验又称回归方程的检验,其目的是检验求得的回归方程在总体中是否成立,即是否样本代表的总体也有直线回归关系。方法有以下两种:(1)方差分析 其基本思想是将总变异分解为SS回归和SS剩余,然后利用F检验来判断回归方程是否成立。(2)t检验其基本思想是利用样本回归系数b与总体均数回归系数ß进行比较来判断回归方程是否成立,实际应

37、用中因为回归系数b的检验过程较为复杂,而相关系数r的检验过程简单并与之等价,故一般用相关系数r的检验来代替回归系数b的检验。3. 直线回归方程的应用(1)描述两变量之间的依存关系; 利用直线回归方程即可定量描述两个变量间依存的数量关系(2)利用回归方程进行预测; 把预报因子(即自变量x)代入回归方程对预报量(即因变量Y)进行估计,即可得到个体Y值的容许区间。(3)利用回归方程进行统计控制 规定Y值的变化,通过控制x的范围来实现统计控制的目标。如已经得到了空气中NO2的浓度和汽车流量间的回归方程,即可通过控制汽车流量来控制空气中NO2的浓度。4. 应用直线回归的注意事项(1)做回归分析要有实际意

38、义;(2)回归分析前,最好先作出散点图;(3)回归直线不要外延。二、直线相关分析1. 直线相关的概念直线相关分析是描述两变量间是否有直线关系以及直线关系的方向和密切程度的分析方法。用以描述两变量间相关关系的指标是相关系数(常用r表示),两变量间相关关系的种类有正相关(0<r<1)、负相关(-1<r<0)、零相关(r=0)、完全相关(|r|=1)等。相关分析对资料的要求是两变量(x,y)均是符合正态分布的随机变量。2. 相关系数的计算相关系数是x,Y的离均差积和lxy除以X的离均差平方和lxx与Y的离均差平方和lyy之积的算术平方根的商。故此相关系数又被称为积差相关系数。

39、3.相关系数的假设检验相关系数检验的目的是判断两变量的总体是否有相关关系,方法有t检验和查表法,t检验法是样本与总体的比较,查表法是直接查相关系数界值表得到相应的概率p。三、直线相关与回归的区别与联系区别:1. 相关说明相关关系,回归说明依存关系;    2. r与b有区别;    3. 资料要求不同。联系:1. r与b值可相互换算; 2. r与b正负号一致; 3. r与b的假设检验等价;4. 回归可解释相关。相关系数的平方r2(又称决定系数)是回归平方和与总的离均差平方和之比,故回归平方和是引入相关变量后总平方和减少的部分。四、等级相关分析等级相关分析

40、适用于资料不是正态双变量或总体分布未知,数据一端或两端有不确定值的资料或等级资料。常用的Spearman等级相关系数rs是利用x,Y的秩次来进行直线相关分析的。因此当x,Y的相同秩次较多时,计算出的rs需矫正。同样的,等级相关系数rs也需要进行假设检验。五、相关分析应用中的注意事项1. 相关分析要有实际意义;2. 相关关系不一定都是“因果”关系;3. 相关系数r假设检验中p的大小不能说明相关的密切程度;4. 直线相关和等级相关有各自不同的适用条件方差分析一、方差分析的基本思想1. 方差分析的概念  方差分析(ANOVA)又称变异数分析或F检验,其目的是推断两组或多组资料的总体均数是否

41、相同,检验两个或多个样本均数的差异是否有统计学意义。我们要学习的主要内容包括单因素方差分析即完全随机设计或成组设计的方差分析和两因素方差分析即配伍组设计的方差分析。2. 方差分析的基本思想  下面我们用一个简单的例子来说明方差分析的基本思想:  如某克山病区测得11例克山病患者和13名健康人的血磷值(mmol/L)如下,患者:0.84 1.05 1.20 1.20 1.39 1.53 1.67 1.80 1.87 2.07 2.11健康人:0.54 0.64 0.64 0.75 0.76 0.81 1.16 1.20 1.34 1.35 1.48 1.56 1.87

42、0; 问该地克山病患者与健康人的血磷值是否不同?  从以上资料可以看出,24个患者与健康人的血磷值各不相同,如果用离均差平方和(SS)描述其围绕总均数的变异情况,则总变异有以下两个来源:(1)组内变异,即由于随机误差的原因使得各组内部的血磷值各不相等;(2)组间变异,即由于克山病的影响使得患者与健康人组的血磷值均数大小不等。  而且:SS总=SS组间+SS组内 v总=v组间+v组内  如果用均方(即自由度v去除离均差平方和的商)代替离均差平方和以消除各组样本数不同的影响,则方差分析就是用组内均方去除组间均方的商(即F值)与1相比较,若F值接近1,则说明各组均数间的

43、差异没有统计学意义,若F值远大于1,则说明各组均数间的差异有统计学意义。实际应用中检验假设成立条件下F值大于特定值的概率可通过查阅F界值表(方差分析用)获得。3. 方差分析的应用条件应用方差分析对资料进行统计推断之前应注意其使用条件,包括:(1)可比性,若资料中各组均数本身不具可比性则不适用方差分析。(2)正态性,即偏态分布资料不适用方差分析。对偏态分布的资料应考虑用对数变换、平方根变换、倒数变换、平方根反正弦变换等变量变换方法变为正态或接近正态后再进行方差分析。(3)方差齐性,即若组间方差不齐则不适用方差分析。多个方差的齐性检验可用Bartlett法,它用卡方值作为检验统计量,结果判断需查阅

44、卡方界值表。二、方差分析的主要内容根据资料设计类型的不同,有以下两种方差分析的方法:1. 对成组设计的多个样本均数比较,应采用完全随机设计的方差分析,即单因素方差分析。2. 对随机区组设计的多个样本均数比较,应采用配伍组设计的方差分析,即两因素方差分析。两类方差分析的基本步骤相同,只是变异的分解方式不同,对成组设计的资料,总变异分解为组内变异和组间变异(随机误差),即:SS总=SS组间+SS组内,而对配伍组设计的资料,总变异除了分解为处理组变异和随机误差外还包括配伍组变异,即:SS总=SS处理+SS配伍+SS误差。整个方差分析的基本步骤如下:(1) 建立检验假设;    H

45、0:多个样本总体均数相等。    H1:多个样本总体均数不相等或不全等。    检验水准为0.05。(2) 计算检验统计量F值;(3) 确定P值并作出推断结果。三、多个样本均数的两两比较  经过方差分析若拒绝了检验假设,只能说明多个样本总体均数不相等或不全相等。若要得到各组均数间更详细的信息,应在方差分析的基础上进行多个样本均数的两两比较。1. 多个样本均数间两两比较  多个样本均数间两两比较常用q检验的方法,即 Newman-kueuls法,其基本步骤为:  建立检验假设->样本均数排序->计算q值->查

46、q界值表判断结果。2. 多个实验组与一个对照组均数间两两比较  多个实验组与一个对照组均数间两两比较,若目的是减小第II类错误,最好选用最小显著差法(LSD法);若目的是减小第I类错误,最好选用新复极差法,前者查t界值表,后者查q界值表。t检验和u检验简而言之,t检验和u检验就是统计量为t,u的假设检验,两者均是常见的假设检验方法。当样本含量n较大时,样本均数符合正态分布,故可用u检验进行分析。当样本含量n小时,若观察值x符合正态分布,则用t检验(因此时样本均数符合t分布),当x为未知分布时应采用秩和检验。一、样本均数与总体均数比较的t检验样本均数与总体均数比较的t检验实际上是推断该

47、样本来自的总体均数µ与已知的某一总体均数µ0(常为理论值或标准值) 有无差别。如根据大量调查,已知健康成年男性的脉搏均数为72次/分,某医生在一山区随即抽查了25名健康男性,求得其脉搏均数为74.2次/分,标准差为6.0次/分,问是否能据此认为该山区成年男性的脉搏均数高于一般成年男性。上述两个均数不等既可能是抽样误差所致,也有可能真是环境差异的影响,为此,可用t检验进行判断,检验过程如下:1. 建立假设H0:µ=µ0=72次/分,H0:µ>µ0,检验水准为单侧0.05。2. 计算统计量进行样本均数与总体均数比较的t检验时t值为样

48、本均数与总体均数差值的绝对值除以标准误的商,其中标准误为标准差除以样本含量算术平方根的商。3. 确定概率,作出判断以自由度v(样本含量n减1)查t界值表,0.025<P<0.05,拒绝H0,接受H1,可认为该山区成年男性的脉搏均数高于一般成年男性。应注意的是,当样本含量n较大时,可用u检验代替t检验。二、配对设计的t检验配对设计是一种比较特殊的设计方式,能够很好地控制非实验因素对结果的影响,有自身配对和非自身配对之分。配对设计资料的t检验实际上是用配对差值与总体均数“0”进行比较,即推断差数的总体均数是否为“0”。故其检验过程与样本均数与总体均数比较的t检验类似,即:1. 建立假设

49、H0:µd=0,即差值的总体均数为“0”,H1:µd>0或µd<0,即差值的总体均数不为“0”,检验水准为0.05。2. 计算统计量进行配对设计t检验时 t值为差值均数与0之差的绝对值除以差值标准误的商,其中差值标准误为差值标准差除以样本含量算术平方根的商。3. 确定概率,作出判断以自由度v(对子数减1)查t界值表,若P<0.05,则拒绝H0,接受H1,若P>=0.05,则还不能拒绝H0。三、成组设计两样本均数比较的t检验成组设计两样本均数比较的t检验又称成组比较或完全随机设计的t检验,其目的是推断两个样本分别代表的总体均数是否相等。其检验

50、过程与上述两种t检验也没有大的差别,只是假设的表达和t值的计算公式不同。两样本均数比较的t检验,其假设一般为:H0:µ1=µ2,即两样本来自的总体均数相等,H1:µ1>µ2或µ1<µ2,即两样本来自的总体均数不相等,检验水准为0.05。计算t统计量时是用两样本均数差值的绝对值除以两样本均数差值的标准误。应注意的是当样本含量n较大时(如大于100时)可用u检验代替t检验,此时u值的计算公式较t值的计算公式要简单的多。四、t检验的应用条件和注意事项两个小样本均数比较的t检验有以下应用条件:(1)两样本来自的总体均符合正态分布,

51、(2)两样本来自的总体方差齐。故在进行两小样本均数比较的t检验之前,要用方差齐性检验来推断两样本代表的总体方差是否相等,方差齐性检验的方法使用F检验,其原理是看较大样本方差与较小样本方差的商是否接近“1”。若接近“1”,则可认为两样本代表的总体方差齐。判断两样本来自的总体是否符合正态分布,可用正态性检验的方法。若两样本来自的总体方差不齐,也不符合正态分布,对符合对数正态分布的资料可用其几何均数进行t检验,对其他资料可用t检验或秩和检验进行分析。实验设计基本概念一、实验设计的概念、特点和分类1.实验设计的概念  将一组随机抽取的实验对象随机分配到两种或多种处理组,观察比较不同处理的效应

52、,这种研究称为“实验研究”.  “实验设计”是为实验研究做的周密计划。2.实验研究的特点(1)能有效地控制误差,节省人、财、物力,提高效率。(2)研究者能人为地设置处理因素。(3)受试对象接受何种处理及处理因素的水平是随机分配的。3.实验研究的分类(1)动物实验:在动物身上进行的实验研究。 (2)临床实验:在临床患病人群中进行的试验研究。(3)社区干预实验:在某地区所有人群中进行的试验研究。二、实验设计的基本要素1.处理因素 又称研究因素,一般是外部施加的因素,确定处理因素时应注意:(1)抓住实(试)验中的主要因素。(2)确定和控制非处理因素。(3)处理因素要标准化。2.受试对象&#

53、160; 即研究的对象,一般有人和动物两类。  人的选择:要注意病人和正常人的正确区分,诊断要明确,受试对象依从性要好等。  动物选择:注意所选动物的种类、品系、年龄、性别、窝别、体重等。3.实验效应 即处理因素的效果,应注意:(1)正确选用观察指标来反映实验效应。所选指标要灵敏、精确、客观,最好选用定量指标。(2)指标观察时应避免偏性,可采用盲法。三、实验设计的基本原则  1.对照原则:即实(试)验要设立对照,使得除实验因素外,对照组与实验组其余因素保持一致,常用的对照有:空白对照、安慰剂对照、标准对照、实验对照、自身对照和历史对照等。  2.重复原则

54、:即研究对象要有一定的数量,或者说样本含量应足够。根据每个具体研究,可有不同的方法来进行样本含量估计。  3.随机化原则:即应保证每个实验对象都有同等机会进入实验或接受某种处理。常用方法有查随机数字表和随机排列表等。随机化是保证均衡性的重要手段。  4.均衡原则:即各处理组非实验因素的条件基本一致,以消除其影响。四、常用的实验设计方法1.完全随机设计将实验对象随机分配至两个或多个处理组去进行实验观察,又称单因素设计、成组设计。 优点:操作简单、应用广泛。 缺点:效率低,只能分析单因素的效应。 资料处理方法:t,u检验,方差分析、秩和检验、卡方检验等。  2.配对(

55、伍)设计 将受试对象配成对子或配伍组,以消除非实验因素的影响。配伍设计又称随机区组设计。配对有自身配对和不同个体配对,配伍实际上是配对的推广。优点:所需样本数和效率均高于成组设计,而且很好地控制了混杂因素的作用。缺点:配对条件不宜满足。资料处理方法:配对t, u检验,秩和检验、配伍组方差分析、配对四格表卡方检验等。  3.其它实验设计方法:(1)交叉设计:在配对设计基础上再加入时间因素,可分析不同阶段的效应。(2)析因设计、拉丁方设计和正交设计等。实验设计的意义、原则与基本内容 一、实验设计的意义 实验设计是科学研究计划内关于研究方法与步骤的一项内容。在医学科研工作中,无论实验室研究

56、、临床疗效观察或现场调查,在制订研究计划时,都应根据实验的目的和条例,结合统计学的要求,针对实验的全过程,认真考虑实验设计问题。一个周密而完善的实验设计,能合理地安排各种实验因素,严格地控制实验误差,从而用较少的人力、物力和时间,最大限度地获得丰富而可靠的资料。反之,如果实验设计存在着缺点,就可能造成不应有的浪费,且足以减损研究结果的价值。总之,实验设计是实验过程的依据,是实验数据处理的前提,也是提高科研成果质量的一个重要保证。二、实验设计的原则实验设计有属于专业方面的,有属于统计方面的。从统计方面说,主要应当考虑对照、重复、随机化等问题,这就是所谓实验设计的三原则。其具体内容我们将在第二、三

57、、四节介绍。三、实验设计的基本内容(一)拟定相互比较的处理所谓处理,指的是在实验研究中欲施加给受试对象的某些因素。如营养实验的各种饲料,治疗某病的几种疗法或药物,药理研究中某药的各种剂量等。在实验的全过程中,处理因素要始终如一保持不变,按一个标准进行实验。如果实验的处理因素是药物,那么药物的成份、含量、出厂批号等必须保持不变。如果实验的处理因素是手术,那么就不能开始时不熟练,而应该在实验之前使熟练程度稳定一致。(二)确定实验对象及数量这里指的是实验所用的动物或活体组织标本等。在实验设计中,要根据实验观察的目的与内容,明确规定采用什么样的实验对象,实验对象中的每个实验单位必须具备的条件与要求,以保证受试对象的一致性。实验对象需要有一定的数量,例数不能太少,也不宜过多。如何估计例数,详见第四节。(三)确定将各实验单位分配到各种处理中去的原则这主要是随机分配或随机化问题。第三节将介绍几种常用的随机分组方法。(四)拟定观察项目和登记表要根据研究目的和任务,

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论