search for: emcluster

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 ,

Dear Knowledgeable R Community Members, Please excuse my ignorance -- I apologize in advance if this is an easy question, but I am a bit stumped and could use a little guidance regarding parameters in 2 functions in the "mclust" package. -------------------- PROBLEM DESCRIPTION -------------------- I have a finite mixture modeling problem -- for example, a 2-component gaussian

Hi All, I am trying to cluster a one-dimensional data (see the file attached) using Mclust() but got an error message like: >Mclust(x) Error in rep(1, n) : Object "n" not found When I do a simulation sometimes it works sometimes doesn't. >Mclust(c(rnorm(50),rnorm(56,-0.5))) Error in rep(1, n) : Object "n" not found >Mclust(c(rnorm(56),rnorm(56,-0.5))) best

Dear group, I want to use the mclust package on large data, and therefore I want to use a subset in the initial clustering phase. From help(mclust): k: If `k' is specified, the hierarchical clustering phase will use a sample of size `k' of the data in the initial hierarchical clustering phase. The default is to use the entire data set. m2 is a

Dear Developers: I have a question about memory management in R 2.2.1 and am wondering if you would be kind enough to help me understand what is going on. (It has been a few years since I have done software development on Windows, so I apologize in advance if these are easy questions.) ------------- MY SYSTEM ------------- I am currently using R (version 2.2.1) on a PC running Windows 2000

Dear All; While using EMclust() in the mclust package, I frequently encountered a program error. A message window popped up with the message " Rgui.exe has generated errors and will be closed by Windows. You will need to restart the program. An error log is be created." > version _ platform i386-pc-mingw32 arch i386 os mingw32

Dear Knowledgeable R Community Members, Please excuse my ignorance, I apologize in advance if this is an easy question, but I am a bit stumped and could use a little guidance. I have a finite mixture modeling problem -- for example, a 2-component gaussian mixture -- where the components have a large overlap, and I am trying to use the "mclust" package to solve this problem. I need

Dear R-users, I have a two - dimensional data set which needs to be clustered into groups: I'm searching for groups of points which show a positive correlation (in a twodimensional plot of the data set), but I do not have any knowledge about how many groups there might be. Do you know of a clustering algorithm in R (or in general) which can use a-priori information about the cluster's

Dear R-List I have a question about how to interpret BIC as a goodness-of-fit statistic. I was trying to use "EMclust" and other "mclust" library and found that BIC was used as a goodness-of-fit statistic. Although I know that smaller BIC indicates a better fit, it is not clear to me how good a fit is by reading a BIC number. Is there a standard way of interpreting a BIC

Hi I have two questions, the first perhaps dumber than the second. Firstly, I have a data set, and when I plot a histogram it looks like a normal distribution. So I want to overlay a bell-shaped normal distribution on top of it, to demonstrate how similar it is to the normal distribution. I have read the help on dnorm(), rnorm(), pnorm() etc but still can't figure out how to plot a normal

Hello. I know my dissimilarity matrix but not my original data. Is there any way i could use the clustering function Mclust or EMclust with this dissimilarity matrix? or at least some equivalent of these functions? As this is model based clustering i dont know if it is actually possible to do it without the original data thanks in advance for your help [[alternate HTML version deleted]]

--Hipsgkxbeg Content-Type: text/plain; charset=us-ascii Content-Description: message body text Content-Transfer-Encoding: 7bit This one was sent privately to me. --Hipsgkxbeg Content-Type: message/rfc822 Content-Description: forwarded message Content-Transfer-Encoding: 7bit Received: from tuvok.kom.tuwien.ac.at (tuvok.kom.tuwien.ac.at [192.35.241.66]) by fangorn.ci.tuwien.ac.at (8.9.3/8.8.5)

Dear Christian, first of all thank you for your answer. I am going to parse through the pages you told me. Meanwhile I'd like to note that probably it is a good idea to put 2-3 lines of R-code demonstrating such a simple needs somnewhere in docs of `cluster' package. E.g. x<-rnorm(500) ... # output means we have rather 1 claster x<-c(rnorm(500), rnorm(500)+5)

...0.1, wascores = TRUE, expand = TRUE, trace = FALSE, plot = FALSE, old.wa = FALSE) ######################### BEGIN EM ANALYSIS ######################### #Use the points determined by MDS to perform EM clustering. #Allow only the unconstrained 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...

Hi, I would like to know if someone knows an extended implementation of k-means in R to find appropriate number of clusters for a given k-dimensional data. Also, I am working on clustering for forecasting, if someone is interested or has knowledge on implementational details please mail me, I would appreciate it. Regards Skanda Kallur "Cogito, ergo sum" (I think, therefore I

I'like to do cluster analysis by using mahalanobis distance. Could you tell me how to do?

I am trying to follow the mclust examples in "MCLUST: Software for Model Based Clustering, Density Estimation and Disriminant Analysis" by Chris Fraley and Adrian Raftery, but I cannot reproduce the density and uncertainty surfaces for the Lansing Woods maples. I am using R 1.8.1 with the code below. The same code works fine in S-Plus 6.2 Am I missing something or is this a bug?

Dear all, I'm trying to obtain a fuzzy clustering with fanny from the cluster package, using a given set of data. That worked just fine. I have another separate sample of data from the same problem. For each case in this new sample I would like to know their membership coefficients with respect to the clustering obtained with the first dataset. In effect I want to have a kind of prediction

I have not followed the history of this thread, but I am quite flummoxed as to why the OP is rewriting code to estimate parameters from an univariate Gaussian mixture model when alternatives such as EMCluster (which generally appears to handle initialization better than MClust) exist. Or perhaps there is more to it in which case I apologize. But I thought that I would make the OP aware of the package (of which, in full disclosure, I am co-author). Best wishes, Ranjan On Sun, 27 Aug 2017 16:01:59 +030...

Hi all, once more about the old subj :-) My data has too much various distribution families and for every particular experiment I need just to decide whether the data is "quite homogeneous" or it has two or more clusters. I've revisited the following libraries: amap, clust, cclust, mclust, multiv, normix, survey. And I didn't find any ready-to-use general