Displaying 1 result from an estimated 1 matches for "probsb1".
Did you mean:
probs1
2009 Sep 21
2
Four concurrent Markov chains
...tionb1[1,2]<-sample(1:6, 1, prob=w, repl=T)
i<-1;observationb2[1,2]<-sample(1:6, 1, prob=w, repl=T)
i<-1;observationb3[1,2]<-sample(1:6, 1, prob=w, repl=T)
i<-1;observationb4[1,2]<-sample(1:6, 1, prob=w, repl=T)
while(i<n+1){
stateb1<-observationb1[i,2];
probsb1<-transitions[stateb1,1:6];
nextstateb1<-rmultinom(1,1,probsb1);
seqb1<-which.max(nextstateb1);
j<-i+1;
observationb1[j,2]<-seqb1;i<-j
#-------------------------------------if I close the brackets here this
works but I want to inlcude this next bit of c...