Displaying 1 result from an estimated 1 matches for "probsb2".
Did you mean:
probs2
2009 Sep 21
2
Four concurrent Markov chains
...om(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[stateb2,1:6];
nextstateb2<-rmultinom(1,1,probsb2);
seqb2<-which.max(nextstateb2);
j<-i+1;
observationb2[j,2]<-seqb2;i<-j
#and so on for birds 3 and 4.
Any help you could give me would be much appreciated.
Lucy
--
View this message in c...