Hi, I once found this somewhere on stackoverflow:
values <- rnorm(20, mean = c(2.15,2.0,2.9), sd = c(0.1,0.1,0.1)) v_dens <- density(values) v_dens_y <- v_dens$y r <- rle(v_dens_y) # These functions ignore the extremes if they are the first or last point maxima_index <- which(rep(x = diff(sign(diff(c(-Inf, r$values, -Inf)))) == -2, times = r$lengths)) minima_index <- which(rep(x = diff(sign(diff(c(-Inf, r$values, -Inf)))) == 2, times = r$lengths)) plot(v_dens_y) HTH Ulrik On Sat, 26 Aug 2017 at 11:49 niharika singhal <niharikasinghal1...@gmail.com> wrote: > I have a Gaussian mixture model with some parameters > > mean=(506.8644,672.8448,829.902) > > sigma=(61.02859,9.149168,74.84682) > > c=(0.1241933, 0.6329082 <06329%20082>, 0.2428986 <02428%20986>) > > And the plot look something like below.[image: enter image description > here] > <https://i.stack.imgur.com/4uUQ9.png> > > Also, if I change my parameters to > > mean=(2.15,2.0,2.9) > > sigma=(0.1,0.1,0.1) > > c=(1/3,1/3,1/3) > > Then plot would change to[image: enter image description here] > <https://i.stack.imgur.com/kESYX.png> > > Is there any way to find the maxima. I have tried Newton's method but it > gave me the wrong output. > > Like in general some common solution, which would work on all the cases, is > needed.Can someone suggest me how can I achieve this > > Thanks in advance > > Niharika Singhal > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.