rm(list = ls())
library(EMCluster)
# source("./R/libContourPlotNew.r")
load("./data/simu.pp.rda")
load("./data/simu.data.rda")
class.true <- da$id
for(k0 in 2:7){
Pi <- da.s.all[[k0]]$model$pi
Mu <- da.s.all[[k0]]$model$Mu
S <- LTSigma2variance(da.s.all[[k0]]$model$LTSigma)
class <- da.s.all[[k0]]$model$class
da <- da.s.all[[k0]]$x
pdf(paste("./plot/figure4_ppcont_k=", k0, ".pdf", sep = ""),
width = 6, height = 6)
plotppcontour(da, Pi, Mu, S, class, class.true)
dev.off()
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.