Claude Messiaen - Urc Necker
2005-Jul-13 10:10 UTC
[R] fitting Weibull distribution on observed percentiles
Hi , R Users
I'm trying to fit a Weibull ditribution on observed percentiles using nls
but it doesn't work. Here is the code I use: is there something wrong ?
# p corresponds to percentiles
# and q to the observed values
# the datas are from the livebirth in france in 1998 distribution
ined1998 <- data.frame( p = c(0.01638,0.49629,0.99284) ,
q = c( 18,27,41))
a0<- 3
b0<- ined1998$q[2]/gamma(1+1/a0)
dist_w1998 <- nls(ined1998$q ~ qweibull(ined1998$p , a0 ,
b0) , start = c(a0,b0) , data=ined1998)
I look forward your reply
C.M
URC NECKER
[[alternative HTML version deleted]]