search for: optimres

...2, PACKAGE = "Matrix")))) .Call("mer_coefGets", x, pars, 2, PACKAGE = "Matrix") .Call("mer_gradient", x, 2, PACKAGE = "Matrix") } else NULL optimRes <- nlminb(.Call("mer_coef", x, 2, PACKAGE = "Matrix"), fn, gr, lower = ifelse(constr, 5e-10, -Inf), control = list(iter.max = value$msMaxIter,...

Dear R-er, For a non-linear optimisation, I used optim() with BFGS method but it stopped regularly before to reach a true mimimum. It was not a problem with limit of iterations, just a local minimum. I was able sometimes to reach better minimum using several rounds of optim(). Then I moved to optimx() to do the different optim rounds automatically using "Nelder-Mead" and

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

I'll not copy all the previous material on this thread to avoid overload. The summary is that all the methods Spencer has tried have some issues. The bad news: This is not uncommon with optimization methods, in part because the problems are "hard", in part because getting them implemented and linked to an interfacing approach like R is very tedious and prone to omissions and

---------- Forwarded message ---------- From: niharika singhal <niharikasinghal1990 at gmail.com> Date: Sun, Aug 27, 2017 at 11:57 AM Subject: Re: Find maxima of a function To: "David Winsemius [via R]" <ml+s789695n4745009h56 at n4.nabble.com>, "Ismail SEZEN [via R]" <ml+s789695n4744993h60 at n4.nabble.com>, Ulrik Stervbo <ulrik.stervbo at gmail.com>

> 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

Hello, could anybody direct me where to find code for optim.R? I was able to find the C code at http://docs.rexamine.com/R-devel/optim_8c.html, but the R version would be easier for me to work with and modify. Thank you! [[alternative HTML version deleted]]

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

*** *** WHAT WOULD MAKE R EASIER TO USE FOR YOU? *** The useR! 2017 conference in Brussels, July 4-7, will include a session in which attendees will be invited to offer their thoughts in at least three areas: * A common interface for similar functions. For example, John C. Nash and others have taken the lead in developing packages to provide one standard call that will

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

Hi R-Help, I am using R to do functional outlier detection (using PCA to reduce to 2 dimensions - the functional boxplot methodology used in the Rainbow package), and using Hscv.diag function to calculate the bandwidth matrix where this line of code is run: result <- optim(diag(Hstart), scv.mat.temp, method = "Nelder-Mead", control = list(trace = as.numeric(verbose))) Within the

...= fl)], function(k) 1:((k * (k + 1))/2) <= k))) devLaplace <- function(pars) .Call(glmer_devLaplace, pars, GSpt) rel.tol <- abs(0.01/devLaplace(PQLpars)) if (cv$msVerbose) cat(paste("relative tolerance set to", rel.tol, "\n")) optimRes <- nlminb(PQLpars, devLaplace, lower = ifelse(const, 5e-10, -Inf), control = list(trace = cv$msVerbose, iter.max = cv$msMaxIter, rel.tol = rel.tol)) .Call(glmer_finalize, GSpt) new("glmer", new("lmer", mer, frame = if (model) mf else dat...

Prof Nash, Prof Galanos Is it possible to use a generic code stub in front of packages that use optimx to improve optimx use or curtail it according to the requirements? Best Regards Amit +91 7899381263 ________________________________________________________________________ Please request Skype as available 5th Year FPM (Ph.D.) in Finance and Accounting Area Indian Institute