Skip to main content

Return to Question

removed snipped wrapper
Source Link
MrFlick
MrFlick
  • 209.1k
  • 19
  • 300
  • 324

I'm fairly new with R and I can't get my head around a for loop. I've developed a big markov model based on a loop (now I know it's not the most efficient way but no time to change it now). the final outcome of the model is a pair of matrices for each value of one parameter j. A sample of what I'm trying to do is 

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
that is: saving in the third dimension of the array the output of every loop. But as I see it saves only the last loop in the last dimension of the array. The model is much more complex so that I can't really change the basic structure but only find a way of saving its outputs. Any thoughts? 

hope it's clear and any help would be really appreciated!

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

that is: saving in the third dimension of the array the output of every loop. But as I see it saves only the last loop in the last dimension of the array. The model is much more complex so that I can't really change the basic structure but only find a way of saving its outputs. Any thoughts? hope it's clear and any help would be really appreciated!

I'm fairly new with R and I can't get my head around a for loop. I've developed a big markov model based on a loop (now I know it's not the most efficient way but no time to change it now). the final outcome of the model is a pair of matrices for each value of one parameter j. A sample of what I'm trying to do is 

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
that is: saving in the third dimension of the array the output of every loop. But as I see it saves only the last loop in the last dimension of the array. The model is much more complex so that I can't really change the basic structure but only find a way of saving its outputs. Any thoughts? hope it's clear and any help would be really appreciated!

I'm fairly new with R and I can't get my head around a for loop. I've developed a big markov model based on a loop (now I know it's not the most efficient way but no time to change it now). the final outcome of the model is a pair of matrices for each value of one parameter j. A sample of what I'm trying to do is 

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

that is: saving in the third dimension of the array the output of every loop. But as I see it saves only the last loop in the last dimension of the array. The model is much more complex so that I can't really change the basic structure but only find a way of saving its outputs. Any thoughts? hope it's clear and any help would be really appreciated!

Source Link
user2008338
user2008338
  • 21
  • 3

For loop just saves the final value in an array

I'm fairly new with R and I can't get my head around a for loop. I've developed a big markov model based on a loop (now I know it's not the most efficient way but no time to change it now). the final outcome of the model is a pair of matrices for each value of one parameter j. A sample of what I'm trying to do is 

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
that is: saving in the third dimension of the array the output of every loop. But as I see it saves only the last loop in the last dimension of the array. The model is much more complex so that I can't really change the basic structure but only find a way of saving its outputs. Any thoughts? hope it's clear and any help would be really appreciated!