similar to: normal mixture EM not working?

I fitted a mixture denstiy of two gaussians two my data. I now want to calculated the standard errors of the estimates via the boot.se command of the mixtools package. My question is now, if the output is correct? It seems a bit odd to me, so is this correct what I am doing and can I rely on the values? My data: http://s000.tinyupload.com/?file_id=09285782882980618119 My code:

I fitted a mixture denstiy of two gaussians two my data. I now want to calculated the standard errors of the estimates via the boot.se command of the mixtools package. My question is now, if the output is correct? It seems a bit odd to me, so is this correct what I am doing and can I rely on the values? My data: http://s000.tinyupload.com/?file_id=09285782882980618119 My code:

I tried to post this question two times, each time it seemed to fail, since " An embedded and charset-unspecified text was scrubbed..." So I try it again: I fitted a mixture density of two gaussians two my data. I now want to caluclate the standard errors of the estimates via the boot.se command of the mixtools package. My question is now, if the output is correct? It seems a bit odd

Dear R-Users, I have a large number of data(54000) and the field of data is 50 to 2.0e9. I want to use normalmixEM (package:mixtools) to fit them in finite mixture narmal distributions,but get some mistakes.I don't know which steps make the error. I have used the following functions before >x<-read.table("data") >log.x<-log10(x$V1) >log.x<-sort(log.x)

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

Moved to R-help because there's no obvious financial content. Michael On Sat, Apr 6, 2013 at 10:56 AM, Stat Tistician <statisticiangermany at gmail.com> wrote: > Hi, > I want to implement the EM algorithm manually, with my own loops and so. > Afterwards, I want to compare it to the normalmixEM output of mixtools > package. > > Since the notation is very advanced, I

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

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

Dear all, i am new in R. I am writing a review paper about batteries. However, i am interested in analyzing all the papers by keywords, author, references and year. This could be done by "refviz" a software, which is only running on windows machines and which is not free. So my question to you is, is it somehow possible to write a script that can do all of this work? And if yes, with

(on-list) Hello Cameron On Fri, 6 Aug 2010 19:02:22 +0100 Liviu Andronic <landronimirc at gmail.com> wrote: > On Fri, 6 Aug 2010 11:30:59 -0600 > Cameron Bracken <cameron.bracken at gmail.com> wrote: > > There is no real good way to deal with this in Sweave. Sweave does > > not actually know anything about the graphics it is creating, all > > it knows it that a

Hi, I tried to send several questions to the lists (both normal R and R-Sig-Finance), but everytime I look them up in the archives my messages end up with the following "An embedded and charset-unspecified text was scrubbed... for example see my post here: https://stat.ethz.ch/pipermail/r-sig-finance/2013q2/011496.html This one was a real important for me. Can subscribers still read it?

I fitted a gaussian mixture to my financial data. The data can be found here: http://uploadeasy.net/upload/32xzq.rar I look at the density with plot(density(dat),col="red",lwd=2) this has a skew of library(e1071) skewness(dat) -0.1284311 Now, I fit a gaussian mixture according to: f(l)=πϕ(l;μ1,σ21)+(1−π)ϕ(l;μ2,σ22) with:

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

Please help me in writing the R code for this problem. I've been solving this for 4 days. It was hard for me to solve it. It's a simulation problem in R. The problem is My true model is a normal mixture which is given as 0.5 N(-0.8,1) + 0.5 N(0.8,1). This model has two components. I will get a sample size equal to 100 from this model. I will do this 200 times. That means, I will have

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

Hi, i would to know, if someone have ever write the code to estimate the parameter (mixing proportion, mean, a var/cov matrix) of a mixture of two multivariate normal distribution. I wrote it and it works (it could find mean and mixing proportion, if I fix the var/cov matrix), while if I fix anything, it doesn't work. My suspect is that when the algorithm iterates the var/cov matrix, something

Hi, i would to know, if someone have ever write the code to estimate the parameter (mixing proportion, mean, a var/cov matrix) of a mixture of two multivariate normal distribution. I wrote it and it works (it could find mean and mixing proportion, if I fix the var/cov matrix), while if I fix anything, it doesn't work. My suspect is that when the algorithm iterates the var/cov matrix, something

Dear all How should I attempt debugging this error? > sosInit() Error: invalid connection When this happens R refuses to run anything: > 2+2 Error: invalid connection > sessionInfo () Error in stdout() : invalid connection Debugging it seems impossible. > options(error=recover) > sosInit() Error: invalid connection Error during wrapup: invalid connection The error also

> No HTML this time. Sorry Dear all, I am attempting to model some one-dimensional data using Gaussian mixture model with mclust.? Generally, the data that I have have 3 overlapping populations (with one of them being the majority, and the other two combining to less than 15%) and for some reason, mclust consistently ignores the smaller peaks, giving me strange values for the means

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Greetings! I tried installing CentOS v7 (1406) to an old spare machine, and the video card apparently don't play well with X. It is incorrectly saying that either monitor I connect is only able to support 640x480, when one is 1024x768 and the other supports a higher (but now forgotten) resolution. In no way do I expect a fix to be added just to