1、1,第三讲 实验设计方法,2,统计设计,实验设计 现场调查设计,3,实验设计方法,1.配对设计(matchedpairs design) 2.交叉试验设计(cross-over design) 3.完全随机组设计(complete random design) 4.配伍组设计(randomized block design) 5.均衡不完全配伍组设计(balanced incomplete blocks design),4,实验设计方法,6.拉丁方设计(Latin-Square design) 7.析因设计(factorial design) 8.正交设计(orthogonal design)
2、 9.序贯实验设计(sequential trial),5,一、配对设计 (matchedpairs design),配对设计是将受试对象配成对子,随机给予每对中的两个体以不同处理。配对条件为主要的非处理因素。 在动物试验中,常将种属、性别、年龄、体重相近的两动物配成对子;临床试验中,常将性别相同,年龄,生活条件、病情轻重等相似的两个病人配成对子。,6,配对设计的形式,自身配对 同一对象接受两种处理,如同一标本用两种方法进行检验,同一患者接受两种处理方法。 异体配对 将条件相近的实验对象配对,并分别给予两种处理。,7,配对设计的特点,配对设计的优点:比较理想地控制了非处理因素的影响,均衡性较好
3、,减少抽样误差。 缺点:观察对象要经过挑选,特别是临床试验中病例较少时,样本含量较少。,8,配对设计的统计分析,当服从正态分布时,用配对比较的t检验;当不服从正态分布或分布未知时,用配对符号秩和检验。,9,二、交叉试验设计 (cross-over design),将A、B两种处理先后施于同一批试验对象,随机地使半数对象先接受A,后接受B;另一半对象先接受B,后接受A。两种处理在全部试验过程中“交叉”进行,称为交叉试验。 由于A和B处于先后两个试验阶段的机会是相等的,因此平衡了试验顺序的影响,而且能把处理方法之间的差别与时间先后之间的差别分开来分析。,10,交叉试验设计,确定受试对象的例数必须是
4、偶数,并编号,尽量使相邻的第1、2号非处理因素近似,第3、4号非处理因素近似,余类推。随机地确定各单号的接受两种处理的顺序,而各双号的顺序,与其前一个单号的顺序相反。 统计分析方法:当服从正态分布时,用方差分析;当不服从正态分布时,用秩和检验。,11,表4-1 16只大自鼠的痛阈值,12,表4-3 表4-1资料的秩和检验,13,交叉试验设计的特点,交叉设计的优点是节省样本数,且均衡性较好。但要求两种观察时间不能过长,处理没有持久效果,使两种处理能够很快区分开。如比较两种抗菌素治疗肺炎的效果,病人一旦痊愈就不能再用另一种抗生素,就不适用。当药物有蓄积作用或排泄缓慢,使得第二次处理不能很快排除第一
5、次的干扰,而影响第二次处理,也不适于应用。,14,三、完全随机组设计 (complete random design),完全随机组设计是将受试对象,随机地分配到各个处理组中进行实验观察,或者从不同总体中随机抽样进行对比观察的一种试验设计。在实验研究中可先将实验对象编号,再用随机数字表或随机排列表把它们随机地分成两组或多组,分别用各种处理进行实验观察,各组实验对象例数可以相等,也可以不相等。,15,完全随机组设计的特点,此种设计的优点是设计和统计分析比较简单。 缺点是试验效率较低,且只能分析一个因素。,16,完全随机组设计的统计分析,一、两处理组 1.大样本(当nl及n2均大于50时) 当资料为
6、计量时可用u检验,当资料为计数时可用2检验。 2.小样本 计量资料:当服从正态分布且方差齐时,采用t检验;当服从正态分布但方差不齐时,采用t检验,当不服从正态分布或分布未知时,采用秩和检验;当资料为计数时可用确切概率法。,17,完全随机组设计的统计分析,二、多个处理组为计量资料时:当服从正态分布且方差齐时,采用方差分析;当不服从正态分布或分布未知时,采用秩和检验。当资料为计数时采用2检验。,18,四、配伍组设计 (randomized block design),配伍设计是配对设计的扩大。当处理组数为三个或三个以上时,将受试对象的非处理因素相同或相似者组成配伍组,每个配伍组的受试对象随机地分配
7、到各个处理组。,19,配伍组设计的特点,各个处理组中的受试对象不仅数量相同,而且比较均衡。这种试验设计,既减小了抽样误差,还可以分析出处理组及配伍组两个因素的影响。 其缺点和配对设计相同,受试对象要经过挑选。,20,缺项估计,配伍组试验的特点是各处理组的受试对象数相等,各配伍组的受试对象数也相等。如在试验过程中,因故造成某个数据丢失,例如试验过程中死去一只动物,或由于仪器故障。此时如果将此配伍组整个去掉,对信息是一种损失。但由于缺少一个数据,无法进行方差分析。可以采用缺项估计加以解决。,21,22,缺项估计自由度的处理,补入的X的值不占自由度,总的自由度应减1,即bt-2。误差自由度也应减1,
8、而处理及配伍的误差自由度不变。,23,五、均衡不完全配伍组设计 (balanced incomplete blocks design),一般每个配伍组的受试对象个数k应等于处理组数v,但有时处理组数v多于配伍组所能容纳的受试对象个数,即vk。此时每个配伍组不能把所有的处埋都安排进去,可采用均衡不完全配伍组设计,简称BIB设计。,24,均衡不完全配伍组设计,设有A、B、C、D四种处理,每个配伍组只能按排3个处理,不能安排所有的4个处理,因此是不完全的。如果按照表51设计,则每个处理因素出现的次数都相同 (3次);且任意两个因素在同一配伍组内的次数也是相同的(2次),因此设计是均衡的。,25,表5
9、1 均衡不完全配伍组设计,26,均衡不完全配伍组设计要求,设v为处理组数,k为每个配伍组受试对象的个数,b为配伍组组数,r为每种处理重复数,为每两种处理同 时出现的配伍组组数。12 必须是整数,27,均衡不完全配伍组设计表,书中附表达式6给出了各种不同情况下的均衡不完全配伍组设计表,选择的原则是由:v : 实验设计所确定的处理因素的水平数。k: 实验条件所能提供配伍组内实验的对象(动物)数或实验同时进行的实验的对象(动物)数,28,六、拉丁方设计 (Latin-Square design),用r个拉丁字母排成r行r列的方阵,使每行及每列中每个字母都只出现一次,这样的方阵称为r阶拉丁方,或rr拉
10、丁方。例:33拉丁方I 1 A B C 2 B C A 3 C A B,29,拉丁方的基本型,第1行及第列都是按某自然的或标准的次序,如按拉丁字母的次序排成的。rr基本型的拉丁方个数Nr为r 2 3 4 5 6Nr 1 1 4 56 9408,30,44拉丁方I 1 A B C D2 B A D C3 C D A B 4 D C B A,31,正交拉丁方设计,A B C B C A C A B A B C B C AC A B,32,拉丁方设计使用与统计分析,拉丁方实验设计用于三个因素,每个因素的水平数相同,且因素间没有交互作用。优点:试验次数减少到最小,并且均衡。 统计分析:方差分析,33,
11、444 拉丁方设计,34,方差分析表,35,拉丁方设计缺项估计,36,七、析因设计 (factorial design),析因设计是多因素的交叉分组设计,即各个因素的所有水平相互组合在一起,形成多个交叉分组。交叉组组数是各因素水平数的乘积。 析因设计不仅可以对每个因数各水平间进行比较,还可以分析因素间的交互作用。,37,交互作用,因素间如存在交互作用,表示各因素不是各自独立的,而是一个因素的水平有改变时,其它因素的效应也随之改变,反之,如不存在交互作用,则表示各因素具有独立性,一个因素的水平有改变时,其它因素的效应不受影响。,38,表7-1 治疗缺铁性贫血四种疗法的红细胞数,39,交互作用,不
12、用B药时(即B保持在Bl的水平):用A药可平均增加红细胞:1.2-0.8=O.4用B药时(即保持在B2的水平):用A药可平均增加红细胞2.1-1.O=1.1 可看出B因素在不同水平下,用A药比不用A药,平均增加红细胞数是不同的,其相差为 1.1-O.4=O.7这就是两药的交互作用,也就是两药具有协同作用。,40,交互作用,一般来说,对于22析因设计(A2B2效应-AlB2效应)-(A2Bl效应-A1B1效应)0 为正交互作用,即协同作用。(A2B2效应-A1B2效应)-(A2B1效应-A1B1效应)0 为负交互作用,即拮抗作用。,41,交互作用,(A2B2效应-A1B2效应)-(A2B1效应-
13、A1B1效应)=0为无交互作用。也可用下面条件表示(A2B2效应+A1B1效应)-(A1B2效应+A2Bl效应)=0为无交互作用;0(0)为正(负)交互作用。当然由于抽样误差的影响,即使上式不为零,也不一定存在交互作用,需要进行假设检验。,42,析因设计的统计分析,22析因设计的统计分析 SS总=SSAB+SS误差= SSA+SSB+SSAB+SS误差 222析因分析设计 SS总=SSABC+SS误差=SSA+SSB+SSC+SSAB+SSAC+SSBC+SSABC +SS误差 3222析因分析设计 SS总=SSABCD+SS误差 =SSA+SSB+SSC+SSD+ SSAB+ SSAC+SS
14、AD+SSBC+ SSBD+ SSCD+ SSABC+ SSABD+SSACD+SSBCD + SSABCD + SS误差,43,22析因设计方差分析,44,22析因设计方差分析,SS总=SSAB+SS误差= SSA+SSB+SSAB+SS误差,45,222析因设计方差分析,46,222析因设计方差分析,手工计算:按书45页生成3个副表计算 SS总=SSABC+SS误差 =SSA+SSB+SSC+SSAB+SSAC+SSBC+SSABC +SS误差方差分析表见书45页,47,222析因设计方差分析,计算机程序计算:,48,八、正交设计 (orthogonal design),正交设计是一种高效
15、、快速的多因素试验方法。它是利用一套规格化的正交表,使每次试验的因素及水平得到合理安排,通过试验结果的分析,获得有用的信息。除了分析主因素外,还可分析交互作用。,49,正交表,正交表是实验设计中合理安排实验并对数据进行统计分析的主要工具。正交表用LN(km)来表示,其中:N表示试验次数k表示安排因素的水平数m表示最多允许安排因素的个数,含主因素与交互因素L代表正交表,50,正交表,例如 L8(27)正交表,51,正交表的特点,每列各水平出现的次数都相同。L8(27)表,各出现4次。 任意两列各种水平的搭配出现的次数也相同。L8(27)表,(1,1)(1,2)(2,1)(2,2)各出现两次。 这
16、种性质称为两列间的正交性,所以称为正交表。,52,设计方法,根据水平数 k 选用相应的表,再根据因素的个数及因素间是否有交互作用,选取多大的表 当交互作用不需考虑时,因素个数m。 当需考虑交互作用时,应选用可分析交互作用的正交表。 L N( km),当N为k的整幂时。将交互作用看成因素,选用合适大小的表。 因素水平不等应选用混合型的正交表。,53,正交表,L8(27)二列间的交互作用,54,L8(27)主效应与与交互作用混杂的 设计表,55,有交互作用的结果分析,直观分析:1. 根据各因素各水平的试验结果,确定各因素好的水平;2. 分析极差,确定因素的重要程度。3. 有时会出现表面上的矛盾。如
17、:A1A2,B2B1 A2B2A1B1,56,有交互作用的结果分析,方差分析 误差项的计算为两部分:1.不满列的空列2.有重复实验的重复 分析过程:1.首先确定主效应。2.其次看两级交互。3.最后看三级交互。,57,有交互作用的结果分析,如 A因素 P0.05 A1A2B、C两因素 PB1C1C2AB P0.05 无交互此时AB水平为A2B2C1 当 AB C P0.05 首先看AB 、 AC、 CB 中P0.05 的。,58,实例分析1,2222无交互(例4-9) 选用L8(27) 表头设计 试验安排 空列合计为误差项,59,60,实例分析2,333无交互(例4-10) 选用L9(34) 表
18、头设计 试验安排 空列合计为误差项,61,62,实例分析3,422无交互(例4-11) 选用L8(424) 表头设计 试验安排 空列合计为误差项,63,64,实例分析4,222有交互(例4-12) 选用L8(27) 表头设计 试验安排 空列合计为误差项,65,表头设计,试验安排,66,实例分析4,由于此设计为满列设计,只进行一次实验不能进行方差分析,因此要重复试验方可进行方差分析,因为此时的误差项是重复试验误差。重复试验次数越多,误差项的MS越小。,67,序贯实验设计,医学科研设计一般都是预先确定样本含量,结果分析在全部试验完成后作出。序贯实验(sequential trial)预先不需要确定
19、试验例数。事先根据实验标准绘好序贯实验图。将每次试验结果的数据点到图上,当达到要求的检验界限,便可停止试验,作出判断。,68,序贯实验设计,序贯设计的类型如下:,69,序贯实验设计,开放型单向一般用于检验药物是否有效;而其它类型一般用于配对比较。 开放型:最多样本数不预先肯定,在逐一试验过程中,究竟多少样本数可使试验终止,视逐一试验结果而定。 闭锁型:试验最多样本数可预先确定,能够试验样本数不超过一定数,必然会终止。单向和双向:是单侧检验还是双侧检验。质反应:二分类资料。量反应:计量资料。,70,序贯实验设计,71,习题,一、选择题 1.为评价某药的降压作用,某单位随机选择50名高血压患者,记
20、录每人服药前后的舒张压值,这种收集资料的过程,称为什么实验设计类型( )。A析因设计 B成组设计 C配对设计 D单组设计 2.对于一个试验而言,配伍组设计优于完全随机设计,这是因为( )。A配伍组设计简单易行 B降低了随机误差C降低了重要非试验因素的影响 D可以分析两因素之间的交互作用,72,3.某麻醉科医生研究催醒宁对氟呱啶的作用,选大鼠6只,按某些条件配成三对。均从腹腔注射氟呱啶40ml/kg,分别在入睡后15min,给每对中的1只腹腔注射催醒宁0.5ml/kg(用A表示)和另1只腹腔注射等量的生理盐水(用B表示),观察出现蹲和走动的时间。当大鼠完全清醒3d后,该医生又用同样的方法重复上述
21、试验,只是将原实验组和对照组的动物互换,结果如下,问:这是一种什么设计类型?A配对设计 B交叉设计 C成组设计 D拉丁方设计,73,4.某地在钩虫病普查中,对粪检阳性者与阴性者各随机抽查10例,测定其血红蛋白量,为了解钩虫粪检阳性者与阴性者的血红蛋白量之间的差别是否有显著性意义。首先需判定此资料取自什么试验设计类型( )。A配伍组设计 B成组设计 C配对设计 D交叉设计,74,5.已知A、B、C都是七水平因素,其中:A是实验因素,B、C都是区组因素,且因素之间的交互作用可忽略不计。试验者希望试验次数尽可能少,问最好选择什么设计类型?A 成组设计 B拉丁方设计 C析因设计 D正交设计,75,6.
22、某医院用中药复方治疗高胆固醇血症,把12例高胆固醇患者随机分为4组,用不同疗法治疗。第1组用一般疗法,第2组在一般疗法上外加甲药,第3组在一般疗法上外加乙药,第4组在一般疗法上外加甲药和乙药,1个月后观察胆固醇降低数(mmol/L),资料如下,问:这是一个什么设计类型?第1组 0.4144,0.6475,0.4662;第2组 1.4504,1.1396,1.0878;第3组 0.7252,0.8029,0.5957;第4组 1.6576,2.0202,2.0720;A单因素4水平 B成组设计 C配对设计 D两因素析因设计,76,7.在析因设计中,必须有( )才便于研究交互作用。A多个因素 B两
23、因素两水平 C三个以上因素 D重复实验 8.某试验涉及A、B等试验因素,经方差分析发现AB的作用有非常显著性意义,这意味着:( )A.因素A或因素B的各水平对观测指标的影响相差很大B.因素A与因素B的各水平都对观测指标的影响相差很大C.因素A与因素B的各水平组合对观测指标的影响相差很大D.因素B的各水平都对观测指标的影响随因素A的水平改变而改变,77,9.在进行实验设计时,需要受试对象分组,如何分组,组与组之间具有最好的可比性?A 多分几组 B将条件相近的分入同一组C将条件相近的分入不同组D将体质强的分入对照组、弱的分入实验组,78,10.许多人习惯用单因素设计取代多因素设计,这样做在什么情况
24、下易得出错误的结论?A样本含量较小时 B实验者操作不太熟练时C因素之间独立时 D因素之间有交互作用时,79,11.已知A、B、C都是三水平因素,且根据预实验得知:AB、CB不可忽视。由于经费紧张,研究者希望试验次数尽可能少一些。试验设计时最好选择( )A 拉丁方设计 B正交设计 C析因设计 D交叉设计,80,二、填空题,1.设在某项试验研究中需要考虑三个因素的影,请选择设计类型: (1)当该三个因素的水平数相等,又不考虑交互作用时,可选用( )。 (2)当该三个因素的水平数分别为2、3、4,又无合适的正交表可供选用时,可选用( ),此时必须有( )。总实验次数至少为( )。这种设计不仅可分析出
25、各因素是否有显著性意义,还可分析出( )是否有显著性意义。,81,2.有人为了说明A、B两台仪器测定的结果之间差别没有显著性意义,居然用一个健康人的4次重复测定值及其统计检验结果为依据,这实际上是曲解了实验设计中的( )原则。 3.为了通过试验了解A、B、C、D 4种药物对某病的疗效。现有7窝同龄、同种属的小鼠,每窝4只,均在同一个时间感染上该病,且病情相同。观察指标为血中MP物质的含量,MP值愈低,表明治疗效果愈好。请按要求填空。 (1)应选择的设计类型为( ); (2)每个药物组各有( )小鼠; (3)按( )将小鼠( )分配到各( )组中去。,82,3.研究苯撑硅橡胶的性能,需要考察A、B、C、D四个因子,各有3个水平,打算做5次重复试验,并希望正交表中试验号尽可能少。若选用L18(37)安排实验,并且交互作用全部忽略,此时SSe由( )部分组成。,