similar to: Mixture of normals

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091125/6a78f78b/attachment-0001.pl>

Hello, I am trying to find out if R can do the following: I have a mixture of normals say f = 0.2*Normal(2, 5) + 0.8*Normal(3,2) How do I find the difference in the densities at any particular point of f and at Normal(2,5)? -- Thanks, Jim. [[alternative HTML version deleted]]

Hi, I have a problem of estimating a mixture of two normal distributions. I need to find the starting points automatically, since this is a part of a larger piece of image processing code. I found the mix2normal1 function in VGAM package that mentions a method of finding starting values for mu1 and mu2 but refers the reader to a book by Everitt and Hand. Unfortunately, I do not have an easy

Dear R-help, I am doing clustering via finite mixture model. Please suggest some packages in R to find clusters via finite mixture model with continuous variables. And also I wish to verify the distributional properties of the mixture distributions by fitting the model with lognormal, gamma, exponentials etc,. Thanks in advance,  warm regards,Ms.Karunambigai M PhD Scholar Dept. of Biostatistics

Dear all, I am attempting to model some one-dimensional data using a mixture model of non-central Student's t distributions. However, I haven't been able to find any R package that provides this functionality. Could there be a way to "manipulate" the EM algorithms from the mixdist or mixtools package to fit the model, or do you have any other suggestions? If anyone could help

Hello I am trying to find an automated way of fitting a mixture of normal and log-normal distributions to data which is clearly bimodal. Here's a simulated example: x.1<-rnorm(6000, 2.4, 0.6)x.2<-rlnorm(10000, 1.3,0.1)X<-c(x.1, x.2) hist(X,100,freq=FALSE, ylim=c(0,1.5))lines(density(x.1), lty=2, lwd=2)lines(density(x.2), lty=2, lwd=2)lines(density(X), lty=4) Currently i am using

Dear list, Can anyone provide a package or code for estimating the parameters of a mixture of c (c >=2) univariate normal distributions? I've tried the algorithm provided by Venables & Ripley (1999) p 263, for the mixture of two normal, but I don't find the "ms" function in R. I've used nls instead, but I'm not sure if it works the same. The data I have is very

Hello, I would like to fit a mixture of two Weibull distributions to my data, estimate the model parameters, and compare the fit of the model to that of a single Weibull distribution. I have used the mix() function in the 'mixdist' package to fit the mixed distribution, and have got the parameter estimates, however, I have not been able to get the log-likelihood for the fit of this model

Hi, I want to fit 3 beta distributions to my data which ranges between 0 and 1. What are the functions that I can easily call and specify that 3 beta distributions should be fitted? I have already looked at normalmixEM and fitdistr but they dont seem to be applicable (normalmixEM is only for fitting normal dist and fitdistr will only fit 1 distribution, not 3). Is that right? Also, my data has 26

Hi, I am currently working on fitting a mixture density to financial data. I have the following data: http://s000.tinyupload.com/?file_id=00083355432555420222 I want to fit a mixture density of two normal distributions. I have the formula: f(l)=πϕ(l;μ1,σ21)+(1−π)ϕ(l;μ2,σ22) my R code is: normalmix<-normalmixEM(dat,k=2,fast=TRUE) pi<-normalmix$lambda[1] mu1<-normalmix$mu[1]

Hello all, I would like to fit a mixture model whose components are normal distributions confined in a closed interval. Since there are already several packages for EM, I would like to extend one of these instead of writing a new script from scratch. What would be the best way to customize such an existing package for doing that? Which one is the best in terms of extensibility? (e.g mixtools,

I want to create the mixture formulation of a discrete uniform ie, say f(x) = 1/10, for i = 1,2,3,4,5,6,7,8,9 and 10 and another discrete distribution which has the same values of x, but he probabilities can vary. Can this be done on any package in R? an if so, can the package estimate the 'probabilities' of each of the x value as well as the mixing proportion if I have the data?

Dear all Please, is there any function or package for dealing with multimodal distributions? I try to fit multimodal distribution or more precisely to find out mixture of normal distributions which can lead to my actual data. I use optim to find (in that case) two parameters but what I want is to let the function find out arbitrary number of normal distributions underlaying my actual data

All: I am trying to fit a mixture of 2 normals with > 110 million observations. I am running R 2.5.1 on a box with 1gb RAM running 32-bit windows and I continue to run out of memory. Does anyone have any suggestions. Thanks so much, Tim [[alternative HTML version deleted]]

hi i'm looking for a gibbs sampling algorithm for R for the case of mixture of K normals, and in particular for the case of bivariate normals. i'd be grateful if anyone could send its own R-routine, at least for the univariate case. thank you in advance matteo

Dear Ingmar & Dave, Thanks a lot for your help and sorry for the late reply. Finally, I've found a way to separate the mixture of distributions (empirically). But the gamlss package looks great, I'm sure it will help me during my further studies. Kind regards, Susanne On 15 Jun 2009, at 20:09, Ingmar Visser wrote: > Dear Susanne & Dave, > > The gamlss package family

Hi, Is there any package that I can use to simulate the Gaussian Mixture Model , which is a mixture modeling method that is widely used in statistical learning theory. I know there is a mclust, however, I think it is a little bit different from my problem. Thanks very much.. regards. -------------------------- Peng Jiang ?? Ph.D. Candidate Antai College of Economics &

Hi, I would like to know if it is possible to have a "R code" to estimate the parameters of a mixture of bivariate (or multivariate) normals via EM Algorithm. I tried to write it, but in the estimation of the matrix of variance and covariance, i have some problems. I generate two bidimensional vectors both from different distribution with their own vector means and variance and

Dear R helper, I hope that u can help me to sort out my problem because I sent an E-mail last night to R-list but I have not receive any help and at the same time I think this problem is not so hard. I have used the following functions before > K<-10 > prime<-c(2,3,5,7,11,13,17) > UN<-seq(1:K)%*%t(sqrt(prime)) > U1<-UN-as.integer(UN) > U<-matrix(qnorm(U1),K,7)

I couldn't find a direct answer in CRAN to this question, so I'm asking with some trepidation. I have a multivariate dataset (data.frame) with columns that can be expressed as a set of mixed normals (at least I think) and need to impute values that have constraints (truncated mixture of normals where the values cannot be below zero). If there isn't a package that can do this, is there