[R] Assistant

Rui Barradas ruipbarradas at sapo.pt
Tue Jan 22 18:17:03 CET 2013


Hello,

You are calling the function in a wrong way. In the case of a beta fit, 
densfun should be the quoted string "beta" and the initial parameter 
values are elements of a named list. Like this:


library(MASS)

x <- rbeta(1000, shape1 = 2, shape2 = 0.5)
fitdistr(x, densfun = "beta", start = list(shape1 = 1, shape2 = 1))


As for your error, I only got it if the data clearly can not fit a beta.

y <- rgamma(1000, shape = 2, rate = 0.5)
fitdistr(y, densfun = "beta", start = list(shape1 = 1, shape2 = 1))
Error in optim(x = c(6.19809706003757, 2.32632108817696, 
3.60844436009277,  :
   initial value in 'vmmin' is not finite


So revise the way you call fitdistr and then, if the error persists, 
revise the parametric distribution to be fitted.


Hope this helps,

Rui Barradas

Em 22-01-2013 16:07, Adelabu Ahmmed escreveu:
> 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:
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>



More information about the R-help mailing list