【发布时间】:2014-09-21 11:23:25
【问题描述】:
我对 R 很陌生,我无法理解 for 循环。我已经开发了一个基于循环的大马尔可夫模型(现在我知道这不是最有效的方法,但现在没有时间改变它)。模型的最终结果是一个参数 j 的每个值的一对矩阵。我正在尝试做的一个示例是
for (j in 1:10){
x<-3:13
y<-rep(x[j],10)
a<-matrix(y,5,2)
y<-rep(2,10)
b<-matrix(y,5,2)
outcome<-array(0,dim=c(5,2,10))
outcome[,,j]<-a*b
}
outcome
即:将每个循环的输出保存在数组的第三维中。但正如我所见,它只保存了数组最后一维中的最后一个循环。该模型要复杂得多,因此我无法真正更改基本结构,只能找到一种保存其输出的方法。有什么想法吗? 希望它很清楚,任何帮助将不胜感激!
【问题讨论】:
-
附言。 Stack sn-ps 尚未针对 R 代码运行。