Re: [R] problem with fitdistr ?

2004-02-18 Thread Prof Brian Ripley
It works under the versions 7.2-0 dated Jan 22 or later: that on CRAN is dated Jan 14 and predates 7.1-14. Since R-devel is `under development', the pieces are not at all times in sync. On 18 Feb 2004, Peter Dalgaard wrote: Prof Brian Ripley [EMAIL PROTECTED] writes: Which version of

Re: [R] problem with fitdistr ?

2004-02-18 Thread Ernesto Jardim
On Tue, 2004-02-17 at 18:59, Prof Brian Ripley wrote: Which version of MASS (not that you gave me any credit)? This works in the current 7.1-14. On Tue, 17 Feb 2004, Ernesto Jardim wrote: I'm trying fitdistr but I'm getting some errors fitdistr(rnorm(100),Normal) Error in

Re: [R] problem with fitdistr ?

2004-02-18 Thread Prof Brian Ripley
On Wed, 18 Feb 2004, Ernesto Jardim wrote: On Tue, 2004-02-17 at 18:59, Prof Brian Ripley wrote: Which version of MASS (not that you gave me any credit)? This works in the current 7.1-14. On Tue, 17 Feb 2004, Ernesto Jardim wrote: I'm trying fitdistr but I'm getting some errors

Re: [R] problem with fitdistr ?

2004-02-18 Thread Peter Dalgaard
Prof Brian Ripley [EMAIL PROTECTED] writes: It works under the versions 7.2-0 dated Jan 22 or later: that on CRAN is dated Jan 14 and predates 7.1-14. Since R-devel is `under development', the pieces are not at all times in sync. OK. I'm still picking up the older version from CRAN on

Re: [R] problem with fitdistr ?

2004-02-18 Thread Prof Brian Ripley
On 18 Feb 2004, Peter Dalgaard wrote: Prof Brian Ripley [EMAIL PROTECTED] writes: It works under the versions 7.2-0 dated Jan 22 or later: that on CRAN is dated Jan 14 and predates 7.1-14. Since R-devel is `under development', the pieces are not at all times in sync. OK. I'm

[R] problem with fitdistr ?

2004-02-17 Thread Ernesto Jardim
Hi, I'm trying fitdistr but I'm getting some errors fitdistr(rnorm(100),Normal) Error in fitdistr(rnorm(100), Normal) : 'start' must be a named list fitdistr(rnorm(100),Normal,start=list(mean=0,sd=1)) Error in fitdistr(rnorm(100), Normal, start = list(mean = 0, sd = 1)) : supplying

Re: [R] problem with fitdistr ?

2004-02-17 Thread Spencer Graves
Which version of R and the MASS library? It worked for me just now with R 1.8.1: library(MASS) fitdistr(rnorm(100),Normal,start=list(mean=0,sd=1)) mean sd 0.03760724 0.97053743 (0.09705374) (0.06862975) hope this helps. spencer graves Ernesto Jardim wrote:

Re: [R] problem with fitdistr ?

2004-02-17 Thread Prof Brian Ripley
Which version of MASS (not that you gave me any credit)? This works in the current 7.1-14. On Tue, 17 Feb 2004, Ernesto Jardim wrote: I'm trying fitdistr but I'm getting some errors fitdistr(rnorm(100),Normal) Error in fitdistr(rnorm(100), Normal) : 'start' must be a named list

Re: [R] problem with fitdistr ?

2004-02-17 Thread Peter Dalgaard
Prof Brian Ripley [EMAIL PROTECTED] writes: Which version of MASS (not that you gave me any credit)? This works in the current 7.1-14. Odd things are happening for me with r-devel, though: library(MASS) fitdistr(rnorm(100),Normal) Error in dn[[2]] : subscript out of bounds

Re: [R] Problem with fitdistr for beta

2003-07-04 Thread Prof Brian Ripley
rbeta(100,0.1,0.1) is generating samples which contain 1, an impossible value for a beta and hence the sample has an infinite log-likelihood. It is clearly documented on the help page that the range is 0 x 1. However, that is not so surprising as P(X 1-1e-16) is about 1% and hence values will

Re: [R] Problem with fitdistr for beta

2003-07-04 Thread Spencer Graves
My standard work-around for the kind of problem you identified is to shrink the numbers just a little towards 0.5. For example: library(MASS) a - rbeta(100,0.1,0.1) fitdistr(x=a, beta, start=list(shape1=0.1,shape2=0.1)) Error in optim(start, mylogfn, x = x, hessian = TRUE, ...) :

Re: [R] Problem with fitdistr for beta

2003-07-04 Thread Prof Brian Ripley
In this example shrinking by (1 - 2e-16) leads to a significant change in the distribution: see my probability calculation. And you can't shrink by much less. A beta(0.1, 0.1) is barely a continuous distribution. On Fri, 4 Jul 2003, Spencer Graves wrote: My standard work-around for the kind