# get data from the package
data(epil1)
fit.BCC <- epil1
# for local cluster
trajplot(fit=fit.BCC,feature.ind=1, which.cluster = "local.cluster",
title= "Local Clustering",xlab="time (months)",
ylab="anxiety",color=c("#00BA38", "#619CFF"))
# for global cluster
trajplot(fit=fit.BCC,feature.ind=1,
which.cluster = "global.cluster",
title="Global Clustering",xlab="time (months)",
ylab="anxiety",color=c("#00BA38", "#619CFF"))
Run the code above in your browser using DataLab