资源描述：
《r语言进行方差分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、方差分析在试验科学中有重要的地位，今天谈谈如何用R做方差分析前提假设：独立、正态、方差齐次（各水平间）例子：x<-c(25.6,22.2,28.0,29.8,24.4,30.0,29.0,27.5,25.0,27.7,23.0,32.2,28.8,28.0,31.5,25.9,20.6,21.2,22.0,21.2)数据集用5个因子水平测量，问是否存在差异光是这样是无法进行分析的，对数据x进行格式转化b<-data.frame(x,a=gl(5,4,20))b得到结果如下（gl指定因子，5是水平，4是重复次数）xa1
2、25.61222.21328.01429.81524.42630.02729.02827.52925.031027.731123.031232.231328.841428.041531.541625.941720.651821.251922.052021.25在进行方差分析之前先对几条假设进行检验，由于随机抽取，假设总体满足独立、正态，考察方差齐次性（用bartlett检验）>bartlett.test(x~a,data=b)Bartletttestofhomogeneityofvariancesdata:xbyaB
3、artlett'sK-squared=7.0966,df=4,p-value=0.1309方差齐次性符合下面进行方差分析m1<-aov(x~a,data=b)summary(m1)DfSumSqMeanSqFvaluePr(>F)a4132.032.994.3060.0162*Residuals15114.97.66—Signif.codes:0‘***’0.001‘**’0.01‘*’0.05‘.’0.1‘’1从这个结果看出差别显著接下来考察具体的差异（多重比较），>TukeyHSD(m1)Tukeymultipl
4、ecomparisonsofmeans95%family-wiseconfidencelevelFit:aov(formula=x~a,data=b)$adifflwruprpadj2-11.325-4.7185827.36858180.95845663-10.575-5.4685826.61858180.99818154-12.150-3.8935828.19358180.80466445-1-5.150-11.1935820.89358180.11405373-2-0.750-6.7935825.29358180
5、.99491814-20.825-5.2185826.86858180.99269055-2-6.475-12.518582-0.43141820.03302404-31.575-4.4685827.61858180.92513375-3-5.725-11.7685820.31858180.06751525-4-7.300-13.343582-1.25641820.0146983除了5、2和5、4间外，其他之间的差异是不显著的。方差分析之二具有交互效应的双因素方差分析昨天谈了单因素方差分析，今天谈谈具有交互效应的双因
6、素分析，举一个实例，例子来源于汤银才所著《R语言与统计分析》。有一个关于检验毒品强弱的试验，给48只白鼠注射3种毒药（因素A），同时有4种治疗方案（因素B），这样的试验在每一种因素组合下都重复四次测试老鼠的存活时间，试分析毒药和治疗方案以及它们的交互作用对老鼠存活时间有无显著影响。数据如下time<-c(0.31,0.45,0.46,0.43,0.82,1.1,0.88,0.72,0.43,0.45,0.63,0.76,0.45,0.71,0.66,0.62,0.36,0.29,0.4,0.23,0.92,0.61,
7、0.49,1.24,0.44,0.35,0.31,0.4,0.56,1.02,0.71,0.38,0.22,0.21,0.18,0.23,0.3,0.37,0.38,0.29,0.23,0.25,0.24,0.22,0.3,0.36,0.31,0.33)tox<-gl(3,16,48)cure<-gl(4,4,48)第二个4代表每种因素组合下的重复试验次数b<-data.frame(time,tox,cure)b结果如下timetoxcure10.311120.451130.461140.431150.821261.
8、101270.881280.721290.4313100.4513110.6313120.7613130.4514140.7114150.6614160.6214170.3621180.2921190.4021200.2321210.9222220.6122230.4922241.2422250.4423260.3523270.3123280.
