Plot a NormMixClus object.
1 2 3 4 5 6 
x 
An object of class 
y_profiles 
y (n x q) matrix of observed profiles for n
observations and q variables to be used for graphing results (optional for

K 
If desired, the specific model to use for plotting. If 
threshold 
Threshold used for maximum conditional probability; only observations with maximum conditional probability greater than this threshold are visualized 
conds 
Condition labels, if desired 
average_over_conds 
If 
graphs 
Graphs to be produced, one (or more) of the following:

order 
If 
profiles_order 
If 
... 
Additional optional plotting arguments 
Andrea Rau, Cathy MaugisRabusseau
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19  ## Simulate toy data, n = 300 observations
set.seed(12345)
countmat < matrix(runif(300*4, min=0, max=500), nrow=300, ncol=4)
countmat < countmat[which(rowSums(countmat) > 0),]
conds < rep(c("A","B","C","D"), each=2)
## Run the Normal mixture model for K = 2,3,4
run_arcsin < coseq(y=countmat, K=2:4, iter=5, transformation="arcsin")
## Plot and summarize results
plot(run_arcsin)
summary(run_arcsin)
## Compare ARI values for all models (no plot generated here)
ARI < compareARI(run_arcsin, plot=FALSE)
## Compare ICL values for models with arcsin and logit transformations
run_logit < coseq(y=countmat, K=2:4, iter=5, transformation="logit")
compareICL(list(run_arcsin, run_logit))

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.