r - R. 我想用循环排除方差分析中的值

Question

我想通过排除它们来检查我用于方差分析的数据点的影响：

wanted =(1:10 !=10)
wanted
summary(aov(log(Sdrtl)~location,subset=wanted))

我不想手动排除 1:10，而是想制作某种循环，以便所有 anova 立即流出。我是初学者，我尝试了以下操作：

for (i in 1:10){
    wanted =(1:10 !=i)
        summary(aov(log(Sdrtl)~location,subset=wanted))
    }

，这并没有给出我想要的。有人可以帮我吗？

Answer 1

我想问题是它summary不打印任何东西 - 它返回一个对象（你忽略它）。

# This will print the summaries
for(i in 1:10) print(summary(aov(log(Sdrtl)~location, subset=1:10!=i)))

# This will return a list of all summaries
x <- lapply(1:10, function(i) summary(aov(log(Sdrtl)~location, subset=1:10!=i)))
x # prints the summaries

第二个版本更有用，因为您可以存储结果并对其进行分析......

Answer 2

您是否检查influence.measures过 base R 中的函数是否已经完成了您想要的操作？它计算一系列遗漏诊断...