实验七双因素方差分析

(整理)实验七双因素方差剖析.(整理)实验七双因素方差剖析.10/10(整理)实验七双因素方差剖析.优选文档实验七双因素方差剖析实验目的：掌握两因素方差剖析的基本理论、统计意义与方法、掌握利用SAS进行模型的建立与明显性检验．实验要求:编写程序，结果剖析．实验内容：写出两因素等重复实验下方差剖析模型，检验步骤；模型：Y—因素A和B,水平A1,A2,,Aa,B1,B2,,Bb,组合水平(Ai,Bj)察看值yij1,yij2,,yijcyijkijijkiiijijk,i1,2,,a,j1,2,，k1,2,,c,bijk~N(0,2),且诸ijk相互独立ababi0,j0,ij0,ij0i1j1i1j1SSTSSASSBSSABSSE检验步骤：（1）统计模型Y—因素A和B,水平A1,A2,,Aa,B1,B2,,Bb,组合水平(Ai,Bj)察看值yij1,yij2,,yijcyijkijijkiiijijk,i1,2,,a,j1,2,，k1,2,,c,bijk~N(0,2),且诸ijk相互独立ababi0,j0,ij0,ij0i1j1i1j1SSTSSASSBSSABSSE（2）明显检验假设：HA0:12a0HA1:最稀有某个i0HB0:12b0HB1:最稀有某个j0优选文档优选文档HAB0:ij0,i1,2,,a,j1,2,,bHAB0:最少ij0HA0真HB0真统计量MSAMSB~(1,(1))FA~F(a1,ab(c1))FBcMSEFbabMSEMSABHA0真FAB~F((a1)(b1),ab(c1))MSE检验p值：如p，拒绝相应假设pAPHA0(FAfA)P(F(a1,ab(c1))fA)pBPHA0(FBfB)P(F(b1,ab(c1))fB)pABPHAB0(FABfAB)P(F(a1)(b1),ab(c1))fAB)（3）无交互效应的各因素均值的估计与比较bci无偏估计iyiyijk/bcj1k1i1置信区间：(yit(ab(c1)MSE/bc)12i1i2置信区间：(yiyi2t(ab(c1)2MSE/bc)，12ma(a1)/2个i1i2的Bonferroni同时置信区间(yi1yi2t(ab(c1)2MSE/bc)1m4）有交互效应时因素各水平组合上的均值估计与比较(Ai,Bj)上均值ij无偏估计：ijyij(i1,2,,a,j1,2,,b)ij的置信区间：(yijt(ab(c1)MSE/c)，12i1ji2j置信度不小于1的Bonferroni同时置信区间(yijyijt(ab(c1)2MSE/c)1212m3.63.7（选作）优选文档优选文档程序以下：1）dataexamp3_6;inputtielizi$jldj$jldjdz@@;cards;a1b10.71a1b11.66a1b12.01a1b12.16a1b12.42a1b12.42a1b12.56a1b12.60a1b13.31a1b13.64a1b13.74a1b13.74a1b14.39a1b14.50a1b15.07a1b15.26a1b18.15a1b18.24a1b22.20a1b22.93a1b23.08a1b23.49a1b24.11a1b24.95a1b25.16a1b25.54a1b25.68a1b26.25a1b27.25a1b27.90a1b28.85a1b211.96a1b215.54a1b215.89a1b218.30a1b218.59a1b32.25a1b33.93a1b35.08优选文档优选文档a1b35.82a1b35.84a1b36.89a1b38.50a1b38.56a1b39.44a1b310.52a1b313.46a1b313.57a1b314.76a1b316.41a1b316.96a1b317.56a1b322.82a1b329.13a2b12.20a2b12.69a2b13.54a2b13.75a2b13.83a2b14.08a2b14.27a2b14.53a2b15.32a2b16.18a2b16.22a2b16.33a2b16.97a2b16.97a2b17.52a2b18.36a2b111.65a2b112.45a2b24.04a2b24.16a2b24.42a2b24.93a2b25.49a2b25.77a2b25.86a2b26.28a2b26.97a2b27.06a2b27.78优选文档优选文档a2b29.23a2b29.34a2b29.91a2b213.46a2b218.40a2b223.89a2b226.39a2b32.71a2b35.43a2b36.38a2b36.38a2b38.32a2b39.04a2b39.56a2b310.01a2b310.08a2b310.62a2b313.80a2b315.99a2b317.90a2b318.25a2b319.32a2b319.87a2b321.60a2b322.25;run;procanovadata=examp3_6;classtielizijldj;modeljldjdz=tielizijldjtielizi*jldj;meanstielizijldj;run;优选文档优选文档样本均值与标准值以下：n=abc=108a=2,b=3则c=18注：SSRSSASSBSSAB，因素变量为自变量，个数为自由度p-1=ab-1=5．先检验H0:ij0，统计量F(5,102)，察看值，检验p值优选文档优选文档由结果可知，,接受HAB0，认为铁离子与剂量等级的交互不明显，即两种形式的铁离子在不同样样剂量下在动物体内的存留量认为同样．而检验HA0和HB0，检验p值分别为pA=0.1446和pB=0.001，接受HA0，拒绝HB0．说明因素A对剂量等级的值的影响不明显。由输出结果可看出，高剂量水平的标准差与其他两个的标准差是有明显性差其他，中剂量和低剂量水平的标准差无明显差异，因此可以认为各组和水平上的标准差有明显相差异，而认为误差的等方差性不合理。2）程序以下：datachange3_6;setexamp3_6;=log(jldjdz);run;procanovadata=change3_6;classjldj;modelaaa=jldj;meansjldj;run;优选文档优选文档由以上输出结果可以看出，各组合水平上的标准差趋于一致。3）程序以下：procanovadata=change3_6;classtielizijldj;modelaaa=tielizijldjtielizi*jldj;meanstielizijldjtielizi*jldj;meanstielizijldj/boncldiffalpha=0.05;run;存留量的方差剖析表自由度方差本源1离子2.07382.07385.990.0161剂量15.58847.794222.5220.0001交互20.81030.40511.170.3143效应误差10235.29590.3460总和10753.7648优选文档优选文档由此结果可看出离子与剂量之间的交互效应是不明显的，也就是说离子的不同样样对剂量差异无影响；同样，剂量之间的差异与铁离子的差异没关。两种离子在小白鼠体内存留量百分比是有明显差其他，不同样样剂量在小白鼠体内存留量百分比也是有明显差其他。4）程序以下：meanstielizijldj/boncldiffalpha=0.05;run;对于专业因素A，由上述结果求得，，，由于仅有一个差12，故m1，，12的置信度为95%的置信区间为（）．因此认为1