data(CouplesCope)
my.states<-StateExpand(CouplesCope, 2:49, 50:97)
my.trans<-StateTrans(my.states, FALSE)
my.logseq<-LogSeq(my.trans)
my.logseq
single.LogSeq(my.logseq, 41) # prints estimates for case 41
Run the code above in your browser using DataLab