4

この種のデータフレームが2つあります

a    b    c
1    2    3
2    3    3
1    2    3

同じ列名で。私は使用しています

mapply(t.test, df1, df2)

2 つのデータフレーム間で a、b、および c 型の列を比較する t.text の結果を出力します。これにより（実際の結果）が得られます。

        LHM                                      
statistic   -11.5671                                 
parameter   90.46322                                 
p.value     1.575918e-19                             
conf.int    Numeric,2                                
estimate    Numeric,2                                
null.value  0                                        
alternative "two.sided"                              
method      "Welch Two Sample t-test"                
data.name   "dots[[1L]][[23L]] and dots[[2L]][[23L]]"
        RaM                                      
statistic   -18.66368                                
parameter   172.2032                                 
p.value     3.200675e-43                             
conf.int    Numeric,2                                
estimate    Numeric,2                                
null.value  0                                        
alternative "two.sided"                              
method      "Welch Two Sample t-test"                
data.name   "dots[[1L]][[24L]] and dots[[2L]][[24L]]"

などなど（各データフレームに〜180列のデータがあります）。列の名前とそれに対応する p 値を行列に格納する必要があります。また、どのデータフレームに高い値が含まれているかを別の列に保存することも役立ちます。ヘルプ？

4

1 に答える 1

5

これを試して：

mapply(function(x, y) t.test(x,y)$p.value, df1, df2)
# on my sample(/dummy) data.frames
#           a           b           c 
# 0.009963864 0.009963864 0.020204103

stack2 列のフォーマットが必要な場合は、次のようにまとめることができますdata.frame。

stack(mapply(function(x, y) t.test(x,y)$p.value, df1, df2))
#        values ind
# 1 0.009963864   a
# 2 0.009963864   b
# 3 0.020204103   c
于 2013-04-15T15:48:44.327 に答える