提供个简单的方法。
灰常简单。
a<-c(-1,-1,-1,1,1,1,1,1,1,-1,1,1,-1,1,1,1)
table(a)
tmp<-embed(a,2)
tmp<-aggregate(matrix(1,length(a)-1),list(tmp[,1],tmp[,2]),sum)
tmp<-matrix(tmp$V1,2,dimnames=list(tmp$Group.1[1:2],tmp$Group.1[1:2]))
#mk状态转移矩阵
mk<-t(apply(tmp,2,function(v) v/sum(v)))
#求姐最终状态
res<-solve(t(mk)-diag(2)+1,c(1,1))