-8

結果1、2、および3の箱ひげ図をマイナーアレル数で割ったものを作成するにはどうすればよいですか?

> data

     mouse.id treatment   outcome1  outcome2  outcome3 snp1 snp2 snp3 snp4 snp5 snp6 snp7 snp8
1        186         2   2427.395   240.635   526.250    0    1    0    0    0    1    0    1
2        186         3   7922.080  3355.925  1786.400    0    1    0    0    0    1    0    1
3        187         1   6114.500  1048.615  1375.990    0    0    1    0    0    0    0    0
4        187         2   2176.345   187.980   631.030    0    0    1    0    0    0    0    0
5        187         3   8523.140  6054.180  2932.915    0    0    1    0    0    0    0    0
4

1 に答える 1

2

ボックスプロットの質問については、構造が非常に過剰に使用されており、それを学習する努力をしていないことがわかったので、わかりません....しかし、集計ステップはおそらく次のようになります。

aggregate(data[, c("outcome1","outcome2","outcome3")], 
           list( rowSums(data[, 6:13]) ), #construct the by-list
           FUN=mean)
#----------
  Group.1 outcome1 outcome2 outcome3
1       1 5604.662 2430.258 1646.645
2       3 5174.738 1798.280 1156.325

おそらくそれを包むだけboxplotですか?

        boxplot( ...[ ,2:4], xlab="Sum of minor alleles")

いいえ、私のテストでは、予想どおりボックスプロットでかなり無知であることが示されています。これは表に表示するのと同じくらい効果的ではないでしょうか?

于 2012-08-10T00:34:04.287 に答える