R/matrixBC.BIC.R

Defines functions matrixBC.BIC

Documented in matrixBC.BIC

matrixBC.BIC <-
function(x,k,r,lambda,alpha=0.2,beta=0.2,nstart=20,Sigma.init=NULL,Delta.init=NULL){
	x<-x-mean(x)
	BIC<-rep(NA,length(lambda))
    nonzero<-rep(NA,length(lambda))

	for(i in 1:length(lambda)){
		mclustering<-matrixBC(x,k,r,lambda=lambda[i],nstart=nstart,alpha=alpha,beta=beta,Sigma.init=Sigma.init,Delta.init=Delta.init)
		BIC[i]<-CalculateBICMatrix(x,mclustering)
		nonzero[i]<-sum(mclustering$Mus!=0)

	}
	return(list(lambda=lambda[which(BIC==min(BIC))[1]],BIC=BIC,nonzeromus=nonzero))
}

Try the sparseBC package in your browser

Any scripts or data that you put into this service are public.

sparseBC documentation built on May 30, 2017, 3:05 a.m.