R 需要帮助从aov中提取多个p值吗

我知道有很多关于如何从

aov

中提取p值的帖子。然而，我有一份有几千个样品的清单。我对每个样本进行了aov，以比较两种不同的治疗方法，现在我正在寻找一种方法来获得一个包含所有p值的列表，因为我无法逐个复制它们

这可能吗？ 对于由

ttest

创建的p值，我这样做没有问题：

results <- apply(data,1,function(x){t.test(x[1:3],x[4:6])$p.value})

---

results
results如果数据已经在列表中，您可能需要
lappy
。您可以使用
summary
从
aov

lapply(yourData, function(x){
    av <- aov(yourFormula, data = x)
    summary(av)[[1]][,5]
})

lappy（数据，函数（x）{
av嗯，apply正在工作……那么使用Lappy（关于结果）有什么区别吗通常，
apply
用于将函数应用于数组或矩阵的边距，而
lappy
在列表中迭代应用函数。我希望您能再次帮助我：当我比较列1-3和列4-6时效果很好，但我无法比较列1-3和列7-8。我真的需要相同的大小吗？因为我只有一个重复的o对于最后一个处理，我怀疑你的线性模型做得不对。你有三个类，前两个一式三份，最后一个一式两份。yu应该使用lm（x~as.因子（c（1,1,1,2,2,3,3））），或类似的词语。尽管如此，您无法估计最终类的类内错误。因此，我无法使用方差分析对这些处理进行比较，因为它们具有不同的重复次数？不是因为它们具有不同的重复次数。因为您不应尝试使用如此少的样本来估计sd。请尝试我刚才使用的lm代码给你，anova（lm（…）会给你你需要的分析，但是你的CI太宽了，以至于毫无意义。你不能复制你的实验/数据采集吗？（如果它解决了你原来的问题，接受答案）顺便说一句，基本统计问题不是材料
results <- apply(data,1,function(x){anova(aov(x[1:3]~x[4:6]))[['Pr(>F)']][1]})

lapply(yourData, function(x){
    av <- aov(yourFormula, data = x)
    summary(av)[[1]][,5]
})