1

分散分析に使用するデータポイントの影響を除外して確認したいと思います。

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

、それは私が欲しいものを与えません。誰かが私を助けることができますか？

4

2 に答える 2

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

2 番目のバージョンは、結果を保存して分析できるため、より便利です...

于 2012-02-01T22:46:37.703 に答える
1

ベース Rの関数が既に目的を果たしているかどうかを確認しましたinfluence.measuresか? リーブワンアウト診断の範囲を計算します...

于 2012-02-01T22:51:26.237 に答える