0

から p 値を抽出する方法についての投稿がたくさんあることは知っていますaov。ただし、数千のサンプルを含むリストがあります。2 つの異なる処理を比較するために各サンプルの aov を実行しましたが、1 つずつコピーできないため、すべての p 値を含むリストを取得する方法を探しています..

これは可能ですか?によって作成されたp値に対してこれを行うのに問題はありませんでしたttest

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

データはインポートした .csv であり、[1:3] は列 [4:6] と比較される列を示します

それは実際には問題ではありませんでしたが、aov に対して同様のことを行うことはできないようです:

results <- apply(data,1,function(x){aov(x[1:3]~x[4:6])})

すべての p 値 (Pr(>F) と呼ばれる) のリストを取得できません。

4

2 に答える 2

0
results <- apply(data,1,function(x){anova(aov(x[1:3]~x[4:6]))[['Pr(>F)']][1]})
于 2014-06-11T18:20:18.703 に答える