Question

我如何制作结果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

Answer 1

我不确定boxplot问题，因为我发现这个构造大大超过了使用而且没有努力学习它，但是聚合步骤可能会是这样的：

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

没有，我的测试显示我在预测的箱形图中相当无能为力。这不会像在桌子上有效地显示一样吗？

用于箱图的R编程

1 个答案: