similar to: R crashes when loading rgl package before minqa package

Hi all, I wrote the following lines of codes try to do some iterations to find the global optimal values, but the function does not execute properly. Every time codes stop after one iteration right after executing the optim() function. Does anyone could have me to take a look? Thanks. if (count>0){ k=k+0.05; mu0=c(83+k,0,0) Sigma0= diag(0.4,3) initpar=c(.1+10*k,10*k,10*k,10*k) # initial

I've been having some strange problems with R CMD check in the last couple of days, but now believe I have localized the issue. I had been running Ubuntu Hardy on one drive and then upgraded to Jaunty, but put Jaunty on a different drive. I continue to be able to boot Hardy when I wish. I soft-linked my R working area i.e., /home/john/Rstuff ----> /media/lnx804/home/john/Rstuff

Hej, When exporting a image from rgl, the following error is encountered: > rgl.postscript('testing.pdf', fmt="pdf") RGL: ERROR: can't bind glx context to window RGL: ERROR: can't bind glx context to window Warning messages: 1: X11 protocol error: GLXBadContextState 2: X11 protocol error: GLXBadContextState The pdf file is created and is readable, but all the labels

Hi! We?ve recently come across an example where a package (minqa) creates an Rcpp Function object in a static variable. This causes R_ReleaseObject to be called by the destructor at a very late point in time - as part of the system exit function: static Function cf("c"); I?m wondering if that is considered to be ?safe?? Is the R engine supposed to stay in a state where calls to API

Hello: The development version of Ecdat on R-Forge contains a vignette in which optim(?, method=?L-BFGS-B?) stops with an error message while violating the lower bound. To see all the details, try the following: install.packages("Ecdat", repos="http://R-Forge.R-project.org") Then do "help(pac=Ecdat)" -> "User guides, package

Dear r-devel readers, I have a pretty deep problem with package loading and unloading in the development version of the lme4 package <https://github.com/lme4/lme4>; it's not boiled down to a properly minimal example yet (this has been difficult), but I am posting anyway in the hopes that someone has ideas about how to proceed farther, since I'm nearly stumped. Apologies in advance

Hej, When plotting triangles with rgl.triangles and setting the axes afterwards with decorate3d(aspect=TRUE), the axes get the color used for the last triangle plotted. Example: rgl.triangles(c(1,2,3),c(1,2,5),c(1,3,2),col="#55FF55") decorate3d(aspect=TRUE) Using decorate3d(aspect=TRUE,col="#000000") or decorate3d(aspect=TRUE,color="#000000") does not help,

Hi there, Can anyone help with this? I tried to install.packages('lme4') and got lots of warning and error messages, starting with the one in the subject and ending with failed compilation. The same went for the 'minqa' package which is needed by lme4. I have R version 3.0.2 running on Fedora 20. "yum update R" says that's the latest available. Is my system

Thanks for all your great answers. The app I?m working on is indeed an exploratory data analysis tool for gene expression, which requires a bunch of bioconductor packages. I guess for now, my best solution is to divide my app into modules and load/unload packages as the user switch from one module to another. This brought me another question: it seems that unload package with the

Hi R users, I am trying to restrct the range of two of the parameters in a maximization problem. Both parameters should be between -1 and 1. As far as I know, if I choose the estimation method ="L-BFGS-B" under Optim, I can restrict the parameter space. However, the "L-BFGS-B" always require finite values of the loglik function and cannot get around of the problem if an

Hi, Mark et al.: Thanks, Mark. Three comments: 1. Rvmmin was one of the methods I tried after Ravi directed me to optimx. It returned NAs for essentially everything. See my email of this subject stamped 4:43 PM Central time = 21:43 UTC. 2. It would be interesting to know if the current algorithm behind optim and optimx with

Hi i´m starting to work with free software (with free as in freedom) i've been working with GAMS (General Algebraic Modeling System) to solve static and dynamic optimization problems and non-linear large systems of equations i'm wondering if there is any free software and documentation on the subject, wether in R, Maxima, or something else i will appreciate any help and any references

Hi All, I have a situation where I want an 'if' variable to be parameterized. It's entirely possible that the way I'm trying to do this is wrong, especially given the error message I get that indicates I can't do this using an 'if' statement. Essentially, I have data where I think a relationship enters when a variable (here Pwd) is below some value (z). I don't

Hi, Currently I am trying to solve a minimization problem using optim as method Nelder-Mead. However, Neldel-Mead needs many iterations until it finally converges. I have set $control.trace and $control.report such that I can see the value of the function at each iteration. I do see that I set the convergence criteria to strict in the sense that the function value does not change much. However,

Hello,   I am currently working on a modeling project using Fortran to run large repetitive loops (many DO loops). As part of this process I would like to use a model fit in R and currently stored as an R object. This is a rather complex model, a GAMM, and I am curious whether there is a way to call this model from Fortran. I am not sure "call" is correct terminology, but I

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

Hello -- The question I have is about the gmm() function from the 'gmm' package (v. 1.4-5). The manual accompanying the package says that the gmm() function is programmed to use either of four numerical solvers -- optim, optimize, constrOptim, or nlminb -- for the minimization of the GMM objective function. I wonder whether there is a way to pass controls to a solver used while calling

I've done a couple of searches and not found mention of the issue below, though some older scripts mention getting absolute paths for R CMD check. If the issue is "new" or unfamiliar I'll be happy to follow up and document it, but suspect it is in some sense already known and I've missed the right search strategy. The workaround is pretty simple -- copy the files to a new

i create inbound confi my confi is: [incoming] exten=> 18888246463,,1,Dial(SIP/8003,60,rT) exten=> 6463,1,Dial(SIP/8003,60,rT) exten=> 18888246463,,n,Wait(5) exten=> 18888246463,,n,Hangup but y calling and send this error in my CLI: [Apr 15 09:58:48] NOTICE[26985]: chan_sip.c:14383 handle_request_invite: Call from '101396_procall' to extension '8888246463' rejected

Hi all, I am trying to rebuild the results of a study using a different data set. I'm using about 450 observations. The code I've written seems to work well, but I have some troubles minimizing the negative of the LogLikelyhood function using 5 free parameters. As starting values I am using the result of the paper I am rebuiling. The system.time of the calculation of the function is