,s)且各xij 相互独立.不考虑两因素的交互作用,因此模型可以归结为其中,为总平均;αi 为因素A 第i个水平的效应;βj 为因素B 第j 个水平的效应.在线性模型下,方差分析的主要任务是系统分析因素A 和因素B 对试验指标影响的大小.因此,在给定显著性水平α 下,提出以下统计假设:对于因素A,“因素A 对试验指标影响不显著”等价于H01:α1=α2=…......
2025-09-30
有一个关于检验毒药强弱的实验,给48只老鼠注射I,II,III三种毒药(因素A),同时有A,B,C和D 四种治疗方案(因素B),这样的试验在每一种因素组合下都重复四次测试老鼠的存活时间,数据见表6-9.试分析毒药和治疗方案以及它们的交互作用对老鼠存活时间有无显著影响.
表6-9 老鼠存活时间(年)的实验数据

(1)根据表6-9以数据框形式导入数据,并用函数“plot( )”作图
>rats<-data.frame(
+Time=c(0.31,0.45,0.46,0.43,0.82,1.10,0.88,0.72,0.43,0.45,
+0.63,0.76,0.45,0.71,0.66,0.62,0.38,0.29,0.40,0.23,
+0.92,0.61,0.49,1.24,0.44,0.35,0.31,0.40,0.56,1.02,
+0.71,0.38,0.22,0.21,0.18,0.23,0.30,0.37,0.38,0.29,
+0.23,0.25,0.24,0.22,0.30,0.36,0.31,0.33),
+Toxicant=gl(3,16,48,labels=c("I","II","III")),
+Cure=gl(4,4,48,labels=c("A","B","C","D"))
+)
>op<-par(mfrow=c(1,2))
>plot(Time~Toxicant+Cure,data=rats)
结果如图6-4所示.

图6-4 毒药和治疗方案两因素的各自效应分析
图6-4显示两因素的各水平均存在较大差异.(https://www.chuimin.cn)
(2)再用函数“interaction.plot( )”作出交互效应图,以考查因素之间交互作用是否存在
>with(rats,interaction.plot(Toxicant,Cure,Time,trace.label="Cure"))
结果如图6-5和图6-6所示.

图6-5 以治疗方案为跟踪变量

图6-6 以毒药为跟踪变量
(3)有交互作用的方差分析
>rats.aov<-aov(Time~Toxicant∗Cure,data=rats)
>summary(rats.aov)
结果如下:

根据p 值知,因素Toxicant和Cure对Time的影响是高度显著的,而交互作用对Time的影响却是不显著的.
(4)再进一步使用前面的Bartlett和Levene两种方法检验因素Toxicant和Cure下的数据是否满足方差齐性的要求


从以上结果可以看到,各p 值均小于0.05,这表明在显著性水平0.05下两因素下的方差不满足齐性的要求,这与图6-4是一致的.
相关文章
 
                                
                                ,s)且各xij 相互独立.不考虑两因素的交互作用,因此模型可以归结为其中,为总平均;αi 为因素A 第i个水平的效应;βj 为因素B 第j 个水平的效应.在线性模型下,方差分析的主要任务是系统分析因素A 和因素B 对试验指标影响的大小.因此,在给定显著性水平α 下,提出以下统计假设:对于因素A,“因素A 对试验指标影响不显著”等价于H01:α1=α2=…......
2025-09-30
 
                                
                                ,βp的线性函数.因此,对μ 作变换,则可得到下面几种分布的连接函数的形式:正态分布m(μ)=μ=∑βi xi.二项分布.Poisson分布m(μ)=lnμ=∑βi xi.上述推广体现在以下两个方面:通过一个连接函数,将响应变量的期望与解释变量建立线性关系m[E]=β0+β1 x1+β2 x2+…......
2025-09-30
 
                                
                                在实验2.3.1中,我们对mtcars数据集进行了展示和描述.在实验2.3.3中,我们对mtcars数据集进行了可视化.现在我们对该数据集中的变量进行回归分析.(1)不考虑变量交互项对该数据集,如果把mpg(汽车每加仑公里数)作为因变量,自变量为hp(马力)和wt(汽车重量)进行回归.>fit<-lm(mpg~hp+wt,data=mtcars)>summary(fit)结果如下:得到的回归方程为......
2025-09-30
 
                                
                                如果是,建立变量之间的定量关系式,并用于预测或控制——回归分析.变量之间的相互关系:分析两组变量之间的相互关系——典型相关分析.多元数据的统计推断这是关于参数估计和假设检验的问题.特别是多元正态分布的均值向量和协方差矩阵的估计和假设检验等问题.多元统计分析的理论基础多元统计分析的理论基础包括多维随机向量,以及由此定义的各种多元统计量,推导它们的分布并研究其性质,研究它们的抽样分布理论.......
2025-09-30
 
                                
                                随着大数据、人工智能在我们日常生活的渗透,学习多元统计分析的人越来越多.“多元统计分析”课程已经被越来越多高校列为相关专业的必修课或选修课.《多元统计分析》教材的特点各有不同,有的教材侧重理论的讲述,读者需要具备较深厚的数学基础;有的教材则注重模型的应用,理论和技术细节不是重点.作者认为,对于侧重“应用”多元统计方法进行数据处理和分析的读者,重点不在于理解多元统计方法的理论证明和公式推导,而是要应......
2025-09-30
 
                                
                                多元线性回归的数学模型可以用矩阵形式表示为1.参数β的最小二乘估计设b0,b1,…表4-1 回归正交设计计算表通过推导,可以得到方差分析表4-2。表4-2 方差分析表3.回归方程的显著性检验F>Fα 若式成立,我们可以认为在显著水平α下,线性回归方程(4-6)是有显著意义的。反之,则认为线性回归方程没有什么意义,这时需要查明原因,具体情况具体分析。......
2025-09-29
 
                                
                                原来检验果汁中含铅量有三种方法A1,A2,A3,现研究出另一种快速检验法A4,能否用A4 代替前三种方法,需要通过实验考察.观察的对象是果汁,不同的果汁当做不同的水平:B1 为苹果,B2 为葡萄汁,B3 为西红柿汁,B4 为苹果饮料汁,B5 桔子汁,B6 菠萝柠檬汁.现进行双因素交错搭配试验,即用四种方法同时检验每一种果汁,其检验结果见表6-8.问因素A和B对果汁的含铅量是否有显著影响?......
2025-09-30
 
                                
                                在实验2.3.2中对iris数据集进行了描述和展示,以下对iris数据集进行可视化.花萼长度、花萼宽度、花瓣长度、花瓣宽度的小提琴图小提琴图是箱线图的变种,因为形状酷似小提琴而得名.小提琴图是将箱线图与核密度图结合在一起,它在箱线图上以镜像方式叠加上核密度图.绘制小提琴图,可以使用vioplot包中的“vioplot( )”函数,但在第一次使用之前请先安装vioplot包.“vioplot( )”函数的调用格式为:vioplot(x1,x2,…......
2025-09-30
相关推荐