similar to: nlm question

Hi I am trying to use nlm with an additional call to nlm within the function but after the first pass, the parameters to the outer call are being passed to the inner call. The inner call is a very trivial problem. ie: test.outer<-function(param.outer){ slope<-nlm(test.inner,param.inner) ... loglikelihood<-sum(...) return(-loglikelihood) } and nlm(test.outer,param.outer) on the

Dear R users, I am trying to minimise a function using "nlm". I am getting the following error message: "Error: function is too long to keep source" The function is really very long (about 100 A4 pages). Is there anything I could do to solve this problem? At the moment I am using "nlmin" in S-Plus with no problems but I'd prefer to use R. Thank you very

I'm a very recent user of R. I have been adapting my Splus programmes and I found only one (important) problem. There exists no function "nlmin" in R and its substitute, "nlm", does not work well with my kind of problems, sometimes no achieving convergence, other tines "converging" to impossible values. My models are highly nonlinear and are to be estimated by

Having compiled and successfully installed e1071 on a previous machine (and run under 0.62.4) I want to compile it on a new machine to run under 0.63.3 but the loader can't find netpbm - both machined are RH5.2 Linux installations. I have some quite old netpbm's on a CD but none more recent and I certainly didn't install them on my old machine. I can't locate or rpm -q or rpm

Having compiled and successfully installed e1071 on a previous machine (and run under 0.62.4) I want to compile it on a new machine to run under 0.63.3 but the loader can't find netpbm - both machined are RH5.2 Linux installations. I have some quite old netpbm's on a CD but none more recent and I certainly didn't install them on my old machine. I can't locate or rpm -q or rpm

Dear all, I am using the functions 'optim' and 'nlminb'. For both, you can provide a function which computes the gradient of the objective function (to enhance speed and precision). In my case, both the objective function and the gradient take time to be computed and share many common computations (similar matrix, products, etc...). Therefore, I have to compute these

I have pulled Jim Lindsey's repeated library and put it under the libraries directory, as are all the other toolboxes. 2 questions: Running R COMPILE *.c and *.f goes fine except that Mathlib.h cannot be found - it is not in the include library. I cannot find any other Mathlib.h on my system anywhere - presumably this is an R specific module. Was this missing from the 0.62.4 distribution -

Hi I have some data (v. large amount) with a (0,1) response where I want to minimise the errors in the betas rather than SS or deviance. So can anyone point me to a ridge regression function or equivalent for such a logistic regression case? John -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

Hello R-ists Have been lurking on the list for a while, I have noticed recently two problems: 1 There seem to be problems with installing R on SuSe Linux version 5.3. As the list only gets problems and doesn't log the successes, has anyone successfully loaded R on SuSe 5.3? Is this a feature of 0.62.3 only or do the problems apply also to earlier versions. 2 There was a list of

Hello and Happy New R Two points: 1 Noting the existence of 0.63.2 as a tgz file on CRAN, but being careful or lazy depending on how you want to see it, I also note that the binaries for Redhat stop at 0.63.1 on RH 5.1. I recently got the RH 5.2 Power Tools where I was pleased to see R 0.62.4 included, lots of libraries including V&R. This had been compiled into an rpm - does anyone know

Hello and Happy New R Two points: 1 Noting the existence of 0.63.2 as a tgz file on CRAN, but being careful or lazy depending on how you want to see it, I also note that the binaries for Redhat stop at 0.63.1 on RH 5.1. I recently got the RH 5.2 Power Tools where I was pleased to see R 0.62.4 included, lots of libraries including V&R. This had been compiled into an rpm - does anyone know

Hi, recently I saw but did not pay too much attention to a question that concerned regression with positive coefficients. In Splus, thereis the nnls() function that can be used if I am not wrong, but what about R ? Now I have the same problem: doing a logistic regression under constraint that coefs are non negative. What can I do with R? is there a (weighted) nnls() counterpart available? Thanks

Hi R-Gurus, I've been cutting along quite nicely with nlm, until I threw in the following condition in the function that nlm is minimising: if (((term*bexp) < 0.0001)) { #warning(term*bexp, "=term*bexp",psi,"=psi") theta<-2000 } Now when I run this function anywhere else, there is no problem, whether or the if's condition is met. When

Greetings, Before I reinvent the wheel has anyone done a profile function for a fit using optim (or nlm)? (like the buggy profile.ms in S+) . It seems a bit tricky as the function to be minimised has to have arguments corresponding to the variables being fitted - which is one less than the function provided to optim()... I guess you can create another function on the fly somehow. Cheers,

Hi Is there a SuSE binary package for R 0.64? John -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch

This is a trivial question probably buried in page 21 somewhere of An Idiot's Guide to R. In MASS(2) P26 for S-Plus you can write: ordered(income)<-c("Lo","Mid","Hi") but in R 1.5.0 you get the response to such an instruction: Error: couldn't find function "ordered<-" What's the fix to this? I have a series of models where the factors

I'm sorry if this is a basic question, but I'm stumped. I'm just trying to plot the residuals from a linear model against another variable in the data frame. Here are the lines I'm trying to execute: size <- read.table(file="/u67/abasl70/surveys/annenberg/mega/smschl.dat", header=T) sizef <- data.frame(size, row.names=size$unit) attach(sizef) mschmod <- lm

I note that the CRAN at .at land still refers to 0.63.2. Is this just a typo? I don't want to download a whole lot only to find I still have it! Does 0.63.3 live elsewhere? And I note the new function kronecker - very useful but can it be applied to non-numerics as in factors, data frames etc preserving the class, levels and names? John

I note that the CRAN at .at land still refers to 0.63.2. Is this just a typo? I don't want to download a whole lot only to find I still have it! Does 0.63.3 live elsewhere? And I note the new function kronecker - very useful but can it be applied to non-numerics as in factors, data frames etc preserving the class, levels and names? John

Despite searching (no doubt in the wrong place) I can't seem to find simple things like the linear predictor or deviance residuals from the glm() function. Obviously it is possible to code particular instances but would I be wasting my time? ie are ther such functions? \John -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read