使用R语言做ANOVA后提取 P值

Posted

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了使用R语言做ANOVA后提取 P值相关的知识,希望对你有一定的参考价值。

如图,我试了用p$的方法,因为有(>F)会报错,求大神指教

参考技术A p$`Pr(>F)`
或者得到线性模型p之后运行
anova(p)追问

谢谢你的回答,我已经解决了,可以用p[[1]][5]

追答

给个采纳吧~我应该没有写错吧(除了没想到你之前不知道为什么多了一层list之外)。。。嘤嘤嘤。。。

本回答被提问者和网友采纳
参考技术B F.value<-p[[1]][4][1,];
p.value<-p[[1]][5][1,];

R语言—方差分析

方差分析简称ANOVA,也称“变异数分析”,用于两个及两个以上样本均数差别的显著性检验。与线性回归不同的是,自变量是因子,研究关注的重点通常会从预测转向组别间差异的分析。基本类型包括:
  • 单因素方差分析ANOVA(组间、组内)
  • 双因素方差分析ANOVA
  • 协方差分析ANCOVA
  • 多元方差分析MANOVA
  • 多元协方差分析MANOVA
(一)模型的拟合
我们通常使用aov()函数对ANOVA模型进行拟合。表中列举了表达式中可以使用的特殊符号。其中y是因变量,字母A、B、C表示因子。

下表中是一些常见的研究设计表达式,其中小写字母表示定量变量,大写字母表示组别因子,Subject 是对被试者独有的标识变量。

R语言—方差分析(1)

(二)表达式中各项的顺序
表达式中效应的顺序在两种情况下会造成影响:(a)因子不止一个,并且是非平衡设计;(b)存在协变量。出现任意一种情况时,等式右边的变量都与其他的变量相关。此时,我们无法清晰的划分它们对因变量的影响。例如,对于双因素方差分析,若不同处理方式中的观测数不同,那么模型y~A*B与模型y~B*A的结果不同。

R语言—方差分析(1)

一般来说,越基础性的效应越需要放在表达式前面。具体来讲,首先是协变量,然后是主效应,接着是双因素的交互项,其次是三因素的交互项,以此类推。
下期我们将使用aov()函数对常见的实验设计进行方差分析,比如说治疗焦虑症、降低胆固醇水平、帮助怀孕小鼠生下胖宝宝、确保豚鼠的牙齿长长、促进植物呼吸、学习如何摆放货架等。我们下期见!


扫描二维码

获取更多精彩

大康的笔记

R语言—方差分析(1)

以上是关于使用R语言做ANOVA后提取 P值的主要内容,如果未能解决你的问题,请参考以下文章

R语言使用cph函数和rcs函数构建限制性立方样条cox回归模型使用anova函数进行方差分析通过p值确认指定连续变量和风险值HR之间是否存在非线性关系

R语言单因素、多因素方差分析ANOVA analysis of variance

R语言使用pwr包的pwr.anova.test函数对单因素方差分析进行效用分析(power analysis)在已知效应量(effect size)显著性水平效用值的情况下计算需要的样本量

R语言嵌套方差分析(Nested ANOVA)实战

R语言基本数据分析

R语言使用anova函数进行方差分析比较两个回归分析模型的差异从而决定是否删除某些预测变量(Comparing nested models using the anova function)