similar to: on the optim function

Dear list, The task view on optimization does not reference a package for non linear constrained optimization problems. Stefan Theussl told me to look at the Rsolnp package, but unfortunately it is not very clear what method is R ported. (The authors ported the matlab code of Yinyu Ye http://www.stanford.edu/~yyye/ <http://www.stanford.edu/%7Eyyye/>) Currently I'm looking for an

Hi all, I'm currently working on the fitdistrplus package (that basically fit distributions). There is something I do not understand about the generic function summary. In the current version on CRAN, there is no NAMESPACE saying S3method(summary, fitdist) . However if we use summary on an object send by fitdist function it works fine... According to R-lang, we have " The most

Hi all, I tried to a nice legend with small boxes filled with the colors used for the plots. But it does nor work, boxes are always filled with black. An example is here plot(1:4,1:4) lines(1:4,4:1, col="blue") legend("top",leg=c("a","b"),col=c("black","blue"), fill=TRUE) How could I specify the colors? the argument col.box is the

Hi, I have a problem when checking the package 'probdistr' (on probability distributions). I got this warning * checking index information ... WARNING Vignettes with missing or empty \VignetteIndexEntry: [1] "probdistr-chi" "probdistr-contextra" "probdistr-discrete" [4] "probdistr-discrextra" "probdistr-exp"

Hi all, I use the symbol functions to draw circles. But the argument lty does not work on windows but works correctly on my macbook: x <- -4:4 y <- -4:4 plot(x,y,type="n") symbols(0, 0, add = TRUE, circles = 1, inches = 1, fg="black", lty="solid") symbols(0, 0, add = TRUE, circles = 1, inches = 2, fg="black", lty="dashed") The

Hi useRs, I'm trying to fit a basic AR(2) model with the 'ar' function. And when I try to check the value of the coefficients, I could not find the same value as the 'ar' function. Here is my example: myserie <- c(212, 205, 210, 213, 217, 222, 216, 218, 220, 212, 215, 236) #plot(myserie, type="l") myserieminus0 <- tail(myserie, -2) myserieminus1 <-

Hi useRs, I'm trying to fit a basic AR(2) model with the 'ar' function. And when I try to check the value of the coefficients, I could not find the same value as the 'ar' function. Here is my example: myserie <- c(212, 205, 210, 213, 217, 222, 216, 218, 220, 212, 215, 236) #plot(myserie, type="l") myserieminus0 <- tail(myserie, -2) myserieminus1 <-

Dear list, Don't you think it could be useful to have in R base a function rowCumSums, that compute cumulative sums for each row of a matrix? My implementation of rowCumSums is rowCumSums <- function(x) t(mapply(function(row)cumsum(x[row,]), 1:NROW(x))) I'm sure it can be improved to have other arguments like na.rm or dims. Is there any hope to have this function in R?

Dear both, I have found that names are not treated in the same way in optim() depending on the optimization method (argument method). The example below shows the difference between the Brent method and the L-BFGS-B method. f <- function(x){ y <- x^2;names(y) <-"f(x)";y} optim(10, f, method="Brent", lower=-1, upper=10)$value optim(10, f, method="L-BFGS-B",

Dear list, I would like to suggest a small update the ?Distributions man page of the stats package. The current version contains the following line. The CRAN package \pkg{SuppDists} for additional distributions. I think it would be better to put in this man page a link to the CRAN task view on Distributions http://cran.r-project.org/web/views/Distributions.html. It is not true that the

Dear useRs, I was wondering how to deal with data where thousands separator is a space or a comma? with google, I thought I found the answer but the thread ( http://tolstoy.newcastle.edu.au/R/help/06/04/25524.html) just focuses on decimal separator. Is there any function to numbers written like this in a csv file? Thanks in advance Christophe -- Christophe DUTANG Ph. D. student at ISFA

Hi R users, I am using the optim funciton to maximize a log likelihood function. My code is as follows: p<-optim(c(-0.2392925,0.4653128,-0.8332286, 0.0657, -0.0031, -0.00245, 3.366, 0.5885, -0.00008, 0.0786,-0.00292,-0.00081, 3.266, -0.3632, -0.000049, 0.1856, 0.00394, -0.00193, -0.889, 0.5379, -0.000063, 0.213, 0.00338, -0.00026, -0.8912, -0.3023, -0.000056), f,

Hi, I would like to draw arrows in a classic 2D plot. Which package should I use? is there R base functions that do job? On google, I could not find any useful discussion about this topic, except a link to the function 'grid.arrows' of the grid package. My problem is I would like to draw arrows at the edge of circles drawn by the 'symbols' function. Maybe there is already a

Dear R users When I use OPTIM with BFGS, I've got a significant result without an error message. However, when I use OPTIMX with BFGS( or spg), I've got the following an error message. ---------------------------------------------------------------------------------------------------- > optimx(par=theta0, fn=obj.fy, gr=gr.fy, method="BFGS", >

In the documentation for 'optim' it gives the following function: fr <- function(x) { ## Rosenbrock Banana function x1 <- x[1] x2 <- x[2] 100 * (x2 - x1 * x1)^2 + (1 - x1)^2 } optim(c(-1.2,1), fr) When I run this code I get: $par [1] 1.000260 1.000506 I am sure I am missing something but why isn't 1,1 a better answer? If I plug 1,1 in the function it seems

Dear all, I have just found a 'good' tutorial R for datamining. I think it should be on the contributed docs. http://cran.r-project.org/other-docs.html Here is the link http://www.liaad.up.pt/~ltorgo/DataMiningWithR/ What do you think? Kind regards Christophe -- Christophe Dutang Ph. D. student at ISFA, Lyon, France website: http://dutangc.free.fr

Hi all, I need to transform classic 32bit Fortran code to 64bit Fortran code, see the discussion [R-SIG-Mac] rnorm.halton. But I'm clearly a beginner in Fortran... Does someone already do this for his package? From here, http://techpubs.sgi.com/library/tpl/cgi-bin/getdoc.cgi?coll=linux&db=bks&fname=/SGI_Developer/Porting_Guide/ch03.html , I identify the following changes

Hi everyone, I am currently using the function optim() to maximize/minimize functions and I would like to see more output of the optimization procedure, in particular the numerical gradient of the parameter vector during each iteration. The documentation of optim() describes that the trace parameter should allow one to trace the progress of the optimization. I use the following command:

Dear R users I am trying to use OPTIMX(OPTIM) for nonlinear optimization. There is no error in my code but the results are so weird (see below). When I ran via OPTIM, the results are that Initial values are that theta0 = 0.6 1.6 0.6 1.6 0.7. (In fact true vales are 0.5,1.0,0.8,1.2, 0.6.) -------------------------------------------------------------------------------------------- >

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20100208/52a6d080/attachment.pl>