2012-02-01 16 views
1

我想通過排除他們檢查我用方差分析的數據點的影響:R.我想排除在ANOVA值與循環

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)) 
    } 

,它不給我想要的。有人可以幫我嗎?

回答

2

我想這個問題是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 

第二個版本是比較有用的,因爲你可以存儲結果,並分析它...

+0

感謝ü!!!!!!!! – user1183859 2012-02-01 22:50:16

1

你檢查,看看是否在基礎R的influence.measures功能已經沒有你想要什麼?它計算了一系列留一出診斷...