Displaying 1 result from an estimated 1 matches for "seqb1".
Did you mean:
seq1
2009 Sep 21
2
Four concurrent Markov chains
...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 code for bird 2
stateb2<-observationb2[i,2];
probsb2<-transitions[state...