我有一个数组可以接受建模包的输出:
M <- array(list(NULL), c(trials,3))
试验是一个可产生大约50组数据的数字。 从采样循环,我插入输出的特定方面。建模包的输出看起来有点像这样:
Mt$effects
c_name effect Other
1 DPC_I 0.0818277549 0
2 DPR_I 0.0150814475 0
3 DPA_I 0.0405341027 0
4 DR_I 0.1255416311 0
5 (etc.)
我通过循环将其插入到我的数组中
For(x in 1:trials) {
Mt<-run_model(params)
M[[x,3]] <- Mt$effects
}
该对象现在看起来如下
M[,3]
[[1]]
c_name effect Other
1 DPC_I 0.0818277549 0
2 DPR_I 0.0150814475 0
3 DPA_I 0.0405341027 0
4 DR_I 0.1255416311 0
5 (etc.)
[[2]]
c_name effect Other
1 DPC_I 0.0717384637 0
2 DPR_I 0.0190812375 0
3 DPA_I 0.0856456427 0
4 DR_I 0.2330002551 0
5 (etc.)
[[3]]
And so on (up to 50 elements).
我想要做的是在这50次试运行中计算每个c_name分组的效果的平均值(和sd),但我无法将数据提取到单个数据帧中(例如)所以我可以对它们进行ddply总结。
我已尝试过rbind,cbind,unlist的各种组合,但我无法理解如何正确地将这些数据从顺序元素中提取出来。我还注意到对.names的任何引用都会导致NULL。 任何解决方案都将非常受欢迎!