fitgmm: gmmfit

In kyoustat/T4Gauss: Tools and Methods for Gaussian Distribution

Description

gmmfit

Usage

fitgmm(data, k = 2, maxiter = 10, method = c("full", "diag"))

Examples

x1 = rmvnorm(50, mean=rep(-1,4))
x2 = rmvnorm(50, mean=rep(+1,4))
xx = rbind(x1, x2)

cl1 = fitgmm(xx, k=1)
cl2 = fitgmm(xx, k=2)
cl3 = fitgmm(xx, k=3)

cls = array(0,c(10,3))
for (i in 1:10){
  cls[i,] = as.vector(fitgmm(xx, k=i, method="diag")$criteria)
  print(paste("iteration ",i," complete..",sep=""))
}

yy = c(rnorm(30,mean=-1),rnorm(30,mean=1))
cls = array(0,c(10,3))
for (i in 1:10){
  cls[i,] = as.vector(fitgmm(yy, k=i, method="diag")$criteria)
  print(paste("iteration ",i," complete..",sep=""))
}

kyoustat/T4Gauss documentation built on April 9, 2020, 10:47 a.m.