results <- BMRMM(foxp2sm, num.cov = 2, simsize = 50,
cov.labels = list(c("F", "W"), c("U", "L", "A")),
duration.distr = list('mixgamma',shape=rep(1,3),rate=rep(1,3)))
fit.summary <- summary(results)
plot(fit.summary, 'trans.probs.mean')
plot(fit.summary, 'dur.mix.probs')
Run the code above in your browser using DataLab