results <- BMRMM(foxp2sm, num.cov = 2, simsize = 80,
duration.distr = list('mixgamma',shape=rep(1,3),rate=rep(1,3)))
diag.BMRMM(results)
diag.BMRMM(results, cov.combs = list(c(1,1),c(1,2)),
transitions = list(c(1,1)), components = c(3))
Run the code above in your browser using DataLab