Я хочу исключить значения в ANOVA с помощью цикла

Question

Я хотел бы проверить влияние точек данных, которые я использую для ановы, исключив их:

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

Вместо исключения 1:10 вручную, я хотел бы сделать какую-то петлю, чтобы все ановы текли одновременно. Я новичок и попробовал следующее:

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 в базе R то, что вам нужно? Он рассчитывает диапазон диагностики «один-на-один» ...