similar to: Fwd: Find maxima of a function

Hi, Thanks for your mail, and time It is not working for some arguments, when mean value is like >6. case mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0. 1134463,0.2074027) rv <-UnivarMixingDistribution(Norm(486.4255, 53.24133), Norm(664.0713, 3.674773), Norm(669.0484, 4.101381),

Dear Niharika, As I said before, the problem is basically an optimization issue. You should isolate the problematic part from the rest of your study. Sometimes, more information does not help to solution. All the answers from us (Ulrik, David, me) are more or less are correct to find a maximum point. Newton?s method is also correct. But after answers, you only say, it didn?t find the right

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

> On 26 Aug 2017, at 16:39, niharika singhal <niharikasinghal1990 at gmail.com> wrote: > > Hi, > > Thanks for your mail, and time > > It is not working for some arguments, when mean value is like >6. > > > case > > mc0 <- c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0. > 1134463,0.2074027) > > rv

Please keep the list in cc. Sorry, it didn't work as expected. Maybe someone else have an appropriate solution. Best, Ulrik On Sa., 26. Aug. 2017, 12:57 niharika singhal <niharikasinghal1990 at gmail.com> wrote: > Hi > > Thanks for you mail, > I really appreciate your time on my problem > > I have posted this problem on > > >

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, 0.2428986) 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

Hello, I am facing a problem with optimization in R from 2-3 weeks. I have some Gaussian mixtures parameters and I want to find the maximum in that *Parameters are in the form * mean1 mean2 mean3 sigma1 sigma2 sigma3 c1 c2 c3 506.8644 672.8448 829.902 61.02859 9.149168 74.84682 0.1241933 0.6329082 0.2428986 I have used optima and optimx to find the

I''m looking to write a javascript heavy clientside program with a something serverside backend that connects to the free maxima math program. I have extensive knowledge of ruby on rails, so I would prefer to call Maxima with ruby, but I don''t know if this is even possible. Its fairly easy to call Maxima (with a lisp implementation) using ANSI C, it is a little less easy to

If a run a LOESS model and then produce a smoothed surface: Is there any way to determine the coordinates of the local maxima on the surface? [[alternative HTML version deleted]]

Dear all, I wonder which R algorithm could perform a search of local maxima in an spatial grid, in other words, having the coordinates of a map (x,y,z... up to 6 coordinates) and then the altitude/height at each point (h) (in total 7 numerical variables) I would like to localise the peaks (local maxima) of that topological surface. Which R algorithm can perform this? Thanks a lot, Josep Maria,

Dear all, suppose I have a bi-variate function f(x,y), I want to find the maxima. I define x and y vector, and get matrix z=f(x,y). how can I get which (x0,y0) makes z become the maxima? I can do two loops to get the x0 and y0, but I think there may exist a function to do this.

Hi. Suppose I have a vector that I partition into disjoint, contiguous subvectors. For example, let v = c(1,4,2,6,7,5), partition it into three subvectors, v1 = v[1:3], v2 = v[4], v3 = v[5:6]. I want to find the maximum element of each subvector. In this example, max(v1) is 4, max(v2) is 6, max(v3) is 7. If I knew that the successive subvector maxima would never decrease, as in the example,

Hello I have a vector v=c(0.08844446,0.1744455,0.1379778,0.1209769,0.1573065,0.1134463,0.2074027) when i do sum(v) or 0.08844446+0.1744455+0.1379778+0.1209769+0.1573065+0.1134463+0.2074027 i am getting output as 1 But if i add them manually i get 1.00000026 I do not want to round of my value since it effect my code further Can anyone suggest how can i avoid this. Thanks & Regards Niharika

I have data in the form of (x,y) pairs and would like to find local minima and maxima (typically the zeros of the 2nd derivative) of the y values. I looked at numericDeriv, but I don't have an "expression" per se. I looked at optim, also, but it looks like it will find only one "global" max or min. I can code up my own piecewise derivatives, but wondered if there is

Is there any R package that can help me with digging out the maxima of a 1-D trajectory ? I have 975 1-D curves. They are only known as time series. That is a set of points ordered with respect to time. Some curves exhibit one only peak. Others have two peaks of different height. We wish to find the number of peaks and their position along the time axis. Apparently it's a trivial problem

Hello, I'm aware that recently there has been discussion of R joining forces with an algebra system in some way. Here is my $0.02. I suggest linking R with Maxima (http://maxima.sf.net). Maxima is written in Lisp. Maxima objects are pretty simple. Maxima works on expressions, and almost every expression is a Lisp list of the form ((op) arg1 arg2 arg3 ...) which associates an operator and

Thank your reply. Pankaj. Actually, I have done it very similar to yours. But I think for my demand, it is better to implement in Front End. Maybe I will re-implement it later in clang. ---------- Forwarded message ---------- From: Pankaj Gode [via LLVM] <ml-node+s1065342n59345h22 at n5.nabble.com> Date: Mon, Jul 15, 2013 at 2:35 PM Subject: Re: Regarding scope information for variable

Hi I have a large dataset which follows a multimodal distribution. And I would like to identify the maxima. As the data is obtained from a stochastic simulation, not all maxima in the data are "real maxima of the dirstribution" but rather small random fluctuations. Unfortunately, it is not possible for me to run more simulations to smooth the obtained distribution. What I am doing

Hello, I have initial parameters for HMM model and I want to use depmixS4 package. The parameters are in the form intial_prob_matrix=matrix(c(0.07614213, 0.45177665, 0.47208122), nrow=1, ncol=3, byrow = TRUE) transition_matrix=matrix(c(0.46666667,0.46666667,0.06666667, 0.06741573,0.5617978,0.37078652, 0.02173913,0.3478261,0.63043478), nrow = 3, ncol =

Dear userRs, "playing around" with combinations of replicate() and random number generating functions inside a self-defined "wrapper" function I encounterd a puzzling behaviour. The following are intentionally simple (and rather nonsense-) examples to isolate the relevant aspects. Please, note the seemingly "inconsistent" behaviour for the second call of