search for: vmmin

I am calling 'vmmin' several times from a C function (which is called via .C). It works very well, except for memory consumption. The cause is that vmmin allocates memory via R_alloc, and this memory is not freed as vmmin exits. Instead all the allocated memory is freed on return of the .C call. In one applicati...

Hi, I'm getting an error message using polr(): Error in optim(start, fmin, gmin, method = "BFGS", hessian = Hess, ...) : initial value in 'vmmin' is not finite The outcome variable is ordinal and factored, and the independant variable is continuous. I've checked the source code for both polr() and optim() and can't find any variable called "vmmin" as referenced in the error message. Any suggestions welcome. Thank...

Dear listers, Am currently using MCMC approaches to estimate some parameters of my model. One parameter has to be updated using a tuned gamma distribution. So at each iteration I estimate the mean and variance of the density of the gamma approximation using "vmmin" (i also supply the gradient argument). For moderate replications the procedure works, but if I increase them R crashes. If instead of the tuned gamma density i use the adaptive rejection sampling procedure by Gilks and Wild (downloaded from their website) the procedure returns the results...

Error in optim(method = "BFGS", c(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, : initial value in 'vmmin' is not finite I am running a logit model with latent class segments. I successfully got estimates for 2 segments. However when I tried to increase the no. of segments, I got this error message at the end. I checked my code again but can't find anything wrong. Is this error message related...

I am writing code in C and would like to call R's functions samin and vmmin (optimization routines: simulated annealing and BFGS) I do not understand how to create and pass in the function (as well as the extra arguments it needs) I am optimizing. I have read the R Extensions manual but it is still unclear to me. Could you give me some pointers and/or direct me to some e...

Dear r helpers, I'm trying to run the models of the BIOMOD package and got this message: * Error in nnet.default (x, y, w,...): **initial value in 'vmmin' is not finite* What does this mean? and how can I correct it? Thank you *Andreia* [[alternative HTML version deleted]]

Dear r-helpers, I'm experiencing some problems in fitting a maximum likelihood binomial model to some of my data. The error is in optim, which founds: Error in optim(par = c(0.2, 0.5), fn = function (p) : initial value in 'vmmin' is not finite Yes, I know it's a common problem, and I've carefully searched and readed all the posts about the issue, But, I can't find a solution. It looks like it will be a small piece of my code, but I'm unable to find it. Maybe anyone could help me... Here is my data and...

Dear r helpers I run the following code for nested logit and got a message that Error in optim(c(0, 0, 0, 0, 0.1, -2, -0.2), fr, hessian = TRUE, method = "BFGS") : initial value in 'vmmin' is not finite What does this mean? and how can I correct it? Thank you June > yogurt = read.table("yogurtnp.csv", header=F,sep=",")> attach(yogurt)> dim(yogurt)[1] 12784 25> choice = yogurt[,2:5]> price=yogurt[,14:17]> feature=yogurt[,6:9]> n = n...

...it running for smaller adjacency matrices (<1,500) it works just fine but when my matrices are bigger 4000+. I get the error: > lnam1_01.adj<- lnam(data01$adopt,x01,ec2001.csr) Error in optim(rho, n2ll.rho, method = method, control = control, beta = parm$beta, : initial value in 'vmmin' is not finite I have looked at the lnam code and cant even figure out what vmmin is. Is there anyway around this? Am I doing something wrong? What makes me think that its about the size of the adjacency matrix is that I can run the same command on similar objects that are just smaller and i...

...;- polr(Npf ~ o_stddev + o_skewness + o_kurtosis + o_acl_1e + dispersal, rlc, Hess=TRUE) And that worked; I got a good model fit. However, a variety of other things that I've tried give me this error: Error in optim(s0, fmin, gmin, method = "BFGS", ...) : initial value in 'vmmin' is not finite This occurs, for example, when I try to use the method="probit" option of polr(). It also occurs when I try a regression involving interactions, such as: pomod <- polr(Npf ~ o_stddev * o_skewness * o_kurtosis * o_acl_1e * dispersal, rlc, Hess=TRUE) I have goo...

...maxRoutine(fn = logLik, grad = grad, hess = hess, start = start, : NA in the initial gradient > p = summary(maxLik(fr,start=c(0,0,0,1,0,-25,-0.2),method="BFGS")) Error in optim(start, func, gr = gradient, control = control, method = "BFGS", : initial value in 'vmmin' is not finite The codes are as follows yogurt = read.table("yogurtnp.csv", header=F,sep=",") attach(yogurt) dim(yogurt) choice = yogurt[,2:5] price=yogurt[,14:17] feature=yogurt[,6:9] n = nrow(yogurt) constant = rep(1,n) yop=cbind(constant,feature[,1],price[,1]) dan=cbind...

[This email is either empty or too large to be displayed at this time]

Good-day Sir, I am R.Language users but am try to? estimate parameter of beta distribution particular dataset but give this error, which is not clear to me: (Initial value in "vmmin" is not finite) beta.fit <- fitdistr(data,densfun=dbeta,shape1=value , shape2=value) kindly assist. expecting your reply:

...s (I think in particular Brian) Especially for larger optimization problems, it would be nice to have an entry point for C/C++ code to the R optimizers (the ones which are called when using optim()), where the client just has to provide the functions fminfn() and fmingr() and calls directly, e.g., vmmin() (all from $RHOME/src/main/optim.c). Are there any plans for providing such an entry point? best Adrian -- Dr. Adrian Trapletti, Olsen & Associates Ltd. Seefeldstrasse 233, CH-8008 Zürich, Switzerland Phone: +41 (1) 386 48 47 Fax: +41 (1) 422 22 82 E-mail: adrian@olsen.ch WWW: http:...

...ge to the given set of data, it requires initial parameters of the distribution, to start off with. Hence, i tried this >x1fitbeta<-fitdistr(x1,densfun=dbeta, start=list(shape1=2,shape2=3)) Error in optim(x = c(1.89074018737135, 1.52649293971978, 2.19950245230280, : initial value in 'vmmin' is not finite I tried with "f" and "chi-square" what i did with "t". Please find below the output. > x1fitt<-fitdistr(x1,"t") Error in fitdistr(x1, "t") : optimization failed In addition: Warning messages: 1: In log(s) : NaNs produced...

...I got outputs even it sounds a problem for delata. Actually, there is a mistake made for the parameters (p & alpha1) set up. If I correct them as: p<-parameters[1] alpha1<-parameters[2] beta1<-parameters[3] delta1<-parameters[4] Then error messga there: initial value in 'vmmin' is not finite In addition: There were 38 warnings (use warnings() to see them). Could you give some advice please? Thanks a lot! Xin Shi [[alternative HTML version deleted]]

Hello: I am sure this question was dealt with several years ago. Is the function vmmin() available from Rmath Standalone? If not is it possible to call optim() or nlm() from Rmath in C. Thank you. Mervyn

...use the function optmin optmin(start,fn,gr,method="L-BFGS-B", lower, upper, hess,...) But then I get the Error in optmin ...: L-BFGS-B needs finite values of fn Then I used optmin(start,fn,gr,method="BFGS", hess, ...) But then I get the Error in optmin ...: initial value in vmmin is not finite I know the final parameter estimates from S-Plus which I use as starting values in R. The upper and lower bounds are close around the final estimates. So there is not much to maximize. What can I do? Thank you for help, Peter [[alternative HTML version deleted]]

Hello, I get a segfault when running glmmboot in my own package glmmML. Has happened many time before, but this time I get no hint of where in my C functions the error might be. I give the output below. Can this be an R bug? I suspect it has to do with repeated calls to 'vmmin' like this: for (...){ vmax = vmaxget(); vmmin(*p, b, &Fmin, bfun, bfun_gr, *maxit, *trace, mask, abstol, reltol, nREPORT, ext, &fncount, &grcount, &fail); vmaxset(vmax); ... } but I am only...

...to do betareg on my dataset. Dependent variable is not normally distributed and is proportion (of condom use (0,1)). But I'm having problems: gyl<-betareg(cond ~ alcoh + drug, data=results) Error in optim(par = start, fn = loglikfun, gr = gradfun, method = method, : initial value in 'vmmin' is not finite Why is R returning me error in optim()? What should I do? In advance, thank you for advice! -- ************************** Vlatka Matkovic Puljic +32/ 485/ 453340 [[alternative HTML version deleted]]