search for: vinv

...d using EM on mclust assuming the data was lognormally distributed and could only get it to work in "E" mode, i.e. the equal variance mode. I could never get it to work on "V" mode [ the warning message was Warning message: In meV(data, z = z, prior = prior, control = control, Vinv = parameters$Vinv, : sigma-squared falls below threshold ] so I gave up even after trying the two suggestions I encountered on the web (eliminating zeros in the data and changing the version of R and mclust). Does anyone understand what's going on here? My EM script in MATLAB has no probl...

Dear All, I have written the following programs to find a non-singular (10*10) covariance matrix. Here is the program: nitems <- 10 x <- array(rnorm(5*nitems,3,3), c(5,nitems)) sigma <- t(x)%*%x inverse <- try(solve(sigma), TRUE) while(inherits(inverse, "try-error")) { x <- array(rnorm(5*nitems,3,3), c(5,nitems)) sigma <- t(x)%*%x inverse <-

...strained models. Sometimes, constrained models mess things up! EMclusters <- mclustBIC(mds$points, G=Clusterrange, modelNames= c("VII", "VVI", "VVV"), prior=NULL, control=emControl(), initialization=list(hcPairs=NULL, subset=NULL, noise=NULL), Vinv=NULL, warn=FALSE, x=NULL) The input data are in the form of an N X N matrix of pairwise genetic distances between strains. Those distances can either be the total number of differences over X characters, or can be normalized to the fraction of characters that differ by dividing the number of diff...

we have implemented teh following code for determinging the clustering model of a dataset. bicvals <- EMclust( hdata, 7) sumry1 <- summary(bicvals, hdata,7) # summary object for emclust() print(sumry1) This set of code gives the following output classification table: 1 2 3 4 5 6 7 1 1 1 4 1 1 1 which I think means there is 1 gene in the 1st cluster...1 gene in the 2nd cluster ,

Hi, I'm trying to use the library mclust for gaussian mixture on a numeric vector. The function Mclust(data,G=3) is working fine but the fitting is not optimal and is using modelNames="E". When I'm trying Mclust(data,G=3,modelName="V") I have the following message: Error in if (Sumry$G > 1) ans[c(orderedNames, "z")] else ans[orderedNames] : argument is