常用实验设计方法(三)六.析因设计(f a c t o r i a l d e s i g n)◆析因设计是一种多因素试验设计。
◆可将两个或多个因素的各个水平进行排列组合,交叉分组进行全面实验。
◆总的实验方案(组合)是各因素水平的乘积。
例如:2×2析因设计(两个因素,每个因素均为2个水平,常可写成22析因设计)A因素(A1、A2)和B因素(B1、B2)共4种实验方案或组合(A1B1、A1B2、A2B1、A2B2)3×3析因设计(两个因素,每个因素均为3个水平,常可写成23析因设计)A因素(A1、A2、A3)和B因素(B1、B2、B3)共9种组合(A1B1、A1B2、A1B3、A2B1、A2B2A2B3、A3B1、A3B2A3B3)2×3×3析因设计(三个因素,一个因素为2个水平,余均为3个水平)共18种组合1.特点①研究的因素个数m≥2,各因素的水平数≥2;②各因素在实验中同时实施且所处的地位基本平等。
③每个因素水平相互组合的实验方案,至少进行2次及以上独立重复实验。
④因素间存在交互效应。
例如,一级(两个因素间)或二级交互(三个因素间)效应。
⑤统计学分析时,各因素及交互项所用误差项是相同的。
◆优点:⏹可分析各因素的主效应(m a i n e f f e c t s)(某因素各水平间的平均效应差异)⏹因素间的交互效应(i n t e r a c t i o n)(一个因素的水平改变会影响另一个因素的效应)⏹寻找最优方案或最佳组合⏹可允许数据缺失(完全随机分配情况下)◆缺点:⏹当因素较多或水平数较多时,所需实验次数过多。
⏹一般来说,因素数最好不要多于6个,水平数亦不要过多,一般为2或3个。
2.析因设计的类型可采用完全随机分配方法或随机区组的析因设计。
可安排两因素或多因素实验⑴2×2析因设计结果见下表:分析:设计类型?如何制定设计方案?如何进行统计学分析?①设计类型两个因素:甲药(不用、用),乙药(不用、用),交叉全面组合,各实验方案独立重复3次,为2×2析因设计。
②制定设计方案将12名贫血患者随机分配到4个实验方案组,同完全随机设计③统计学分析分析:不用甲药:乙药的效应(两个水平间差异=3.0-2.4=0.6)用甲药:乙药的效应(两个水平间差异=6.3-3.6=2.7)提示:可能存在交互效应。
均数图提示:可能存在交互效应(若无交互效应,两条线接近平行)数据录入格式(d a t a2/析因设计1.x l s)编号甲药乙药红细胞数1000.82000.93000.74101.35101.26101.17010.98011.19011.010112.111112.212112.0统计学分析采用析因设计的方差分析,总变异被划分为A因素的效应、B因素的效应、A B的交互效应及误差。
具体见下表:两因素析因设计方差分析表变异来源S S(离均差平方和)d f(自由度)M S(均方)M S=S S/d fF值P值A因素S S A a-1M S A M S A/M S误差B因素S S B b-1M S B M S B/M S误差A B交互S S A×B(a-1)(b-1)M S A×B M S A×B/M S误差误差S S误差N-1-(a-1)-(b-1)-(a-1)(b-1)M S误差总变异S S总N-1假设A因素水平数为a,B因素水平数为b,试验组合数为a×b,重复k次,样本量为N=a×b×k联合效应。
例2:甘蓝叶中测定核黄素含量,数据如下:数据见(d a t a2/析因设计2.x l s)⑵2×2×2析因设计例3:研究实验动物的性别和不同饲料(玉米和大豆粉)对体重增加的影响。
性别A(A1=雌,A2=雄)玉米B(B1=加乙氨酸、B2=不加乙氨酸)大豆粉C(C1=加蛋粉、C2=不加蛋粉)三因素共8种试验方案,重复8次,共需64只动物。
数据结果如下:试进行统计分析①数据录入格式(d a t a2/析因设计3.x l s)编号A B C增重量11110.5421110.5531110.7441110.7151110.6261110.5871110.5181110.56101120.61 111120.77 121120.6 131120.58 141120.79 151120.57 161120.72 171210.43 181210.5 191210.58 201210.65 211210.51 221210.57 231210.68 241210.66 251220.52 261220.49 271220.49 281220.49 291220.48 301220.61 311220.59 321220.62 332110.69 342110.54 352110.7 362110.61 372110.73 382110.7 392110.59 402110.61 412120.76 422120.73 432120.63 442120.61 452120.84 462120.62 472120.67 482120.66 492210.61 502210.57 512210.67 522210.71542210.6552210.63 562210.66 572220.55 582220.48 592220.54 602220.49 612220.42 622220.6632220.64 642220.48②统计分析结果结果:B 和C 存在交互作用,B 1C 2组合平均日增重量最高。
例4 :钩端螺旋体培养实验,数据如下,试进行统计分析⑶ 含区组因素的析因设计● 设计具有析因设计的全部特点,还比析因设计多考虑一个区组因素对观察结果的影响。
● 设计时,先列出析因设计实验因素各水平的全部组合,将受试对象按某些重要的混杂因素形成一个区组,再随机分配到各处理组合中。
● 随机分配方案:同随机区组设计例5 将肉瘤接种在大白鼠的3个不同部位(b 1、b 2、b 3),并用药物(a 1、a 2)处理,观察肉瘤的生长情况,共有6种处理组合,将24只大白鼠按体重配成4个区组,每组6只大白鼠,随机分配到6个处理方案中。
结果如下:药物(a1、a2)对大白鼠不同部位的肉瘤生长的影响区组瘤重a1b1a1b2a1b3a2b1a2b2a2b313.08.05.010.05.46.521.16.57.212.07.28.034.23.28.09.59.27.043.25.59.07.06.33.2试回答:设计类型?如何制定设计方案?如何进行统计学分析?①随机区组2×3析因设计②设计方案:将区组中的每个受试对象随机分配到6种不同的处理中,同随机区组设计方案。
③统计分析:采用随机区组的析因设计的方差分析,总变异被划分为区组效应、A因素效应、B因素的效应、A B交互效应和误差。
数据录入格式(d a t a2.x l s)区组A B瘤重11132111.13114.24113.211282126.53123.24125.511352137.23138413912110221123219.542171225.42227.23229.24226.31236.5223832374233.2统计分析结果结果:A B间存在交互作用,b1a2情况下瘤重最重,b1a1瘤重最轻。
(注意:无重复的析因设计只能分析主效应,不能分析交互效应。
)案例1现有3种降低转氨酶的药物A、B、C,为了考察它们对甲型肝炎和乙型肝炎患者转氨酶降低程度之间的差异是否存在统计学差异。
研究者从两型患者的总体中各随机抽取9例,然后随机均分到3个药物组中,实验数据如下:药物转氨酶降低值肝炎类型:甲型肝炎乙型肝炎A10085906575100B12090110453050C506040506045案例2某临床医师用两种药物治疗某种疾病患者各15人,同时观测到可能的影响因素还有患者的性别、年龄和体重,观测的定量指标有M N C(108/k g)、C D34+(106/k g)、C D3+(108/k g)、C D56+(108/k g),资料收集如下。
试问,仅考虑药物类型和性别分别对4个定量观测指标的影响,其设计类型是什么?编号药物性别年龄(岁)体重(k g)M N C C D34+C D3+C D56+ 1甲药131604.427.072.640.30 2甲药243582.671.391.990.10 3甲药155584.142.152.000.37 4甲药155583.231.58..5甲药235602.541.091.720.33 6甲药124582.371.421.440.12 7甲药137602.380.48..8甲药137602.581.55..9甲药243604.542.952.170.1810甲药126601.240.31..11甲药238682.433.43..12甲药229732.161.191.500.1613甲药146733.494.361.480.3814甲药143853.065.51..15甲药146852.652.411.220.32 1乙药238553.864.981.710.14 2乙药116466.005.881.611.66 3乙药228584.573.661.670.68 4乙药23063.021.961.270.53 5乙药232603.752.663.240.63 6乙药238605.419.201.911.46 7乙药138682.683.640.970.46 8乙药138682.733.06..9乙药146563.993.83..10乙药146563.841.151.51.11乙药120605.796.542.121.4312乙药120605.233.141.611.5713乙药249573.422.331.331.0714乙药136674.381.931.211.87 15乙药243757.608.362.651.66案例3某研究者将75只小鼠完全随机的均分为3组,每组25只,分别在5.5G y剂量照射前给I L-11药、照射后给I L-11药和照射对照(即不给I L-11药),然后将每组中的25只小鼠随机地均分为5组,每组5只,分别在照后6小时、12小时、24小时、3天、1周5个时间点上,各观察I L-11对5.5G y照射小鼠骨髓造血细胞周期(G0/G1期)的影响。
数据如下,试问这是一个什么实验设计?干预重复编号骨髓造血细胞周期(G0/G1期)(单位)时间:6h12h24h3d1w e e k照射前给药166.8086.6395.8866.4366.74266.7284.6993.6767.1660.87366.5984.6993.6772.7459.90462.8287.9294.9380.6953.58563.6986.1090.8477.0653.41照射后给药159.3080.6995.5370.5454.56257.3486.7091.1468.1064.26358.3179.6795.0567.0462.22460.1684.6793.2668.9958.81560.4683.3093.6774.0557.86照射对照165.5886.5695.2375.1965.20262.6384.4292.1474.2665.79364.0387.4198.8474.5460.10467.5983.5096.1174.3260.47564.8985.4795.6574.5862.89案例4将24只豚鼠作为受试对象,按初始体重由轻至重排序,组成6个区组,保证每组中的4只豚鼠体重接近,然后将其随机分配至4个处理(饲料)组,观测豚鼠平均日增重量(g),数据如下,分析这是什么实验设计?区组平均日增重量(g)A(4%蛋白质)B(普通饲料)C(4%蛋白质+6%己氨酸)D(6%己氨酸)11.381.521.221.11 21.091.480.871.03 31.081.451.130.97 41.091.221.000.97 51.401.271.341.09 61.471.531.160.99。