# NOT RUN {
data(toydataset)
x <- toydataset$x
Y <- t(toydataset[,2:ncol(toydataset)])
mixhmmr <- emMixHMMR(X = x, Y = Y, K = 3, R = 3, p = 1, verbose = TRUE)
# mixhmmr is a ModelMixHMMR object. It contains some methods such as 'summary' and 'plot'
mixhmmr$summary()
mixhmmr$plot()
# mixhmmr has also two fields, stat and param which are reference classes as well
# Log-likelihood:
mixhmmr$stat$loglik
# Parameters of the polynomial regressions:
mixhmmr$param$beta
# }
Run the code above in your browser using DataLab