Displaying 1 result from an estimated 1 matches for "llfunc".
Did you mean:
lfunc
2003 Oct 20
1
Fitting a Weibull/NaNs
...o fit a Weibull distribution to some data via maximum
likelihood estimation. I'm following the procedure described by Doug
Bates in his "Using Open Source Software to Teach Mathematical
Statistics" but I keep getting warnings about NaNs being converted to
maximum positive value:
> llfunc <- function (x) { -sum(dweibull(AM,shape=x[1],scale=x[2], log=TRUE))}
> mle <- nlm(llfunc,c(shape=1.5,scale=40), hessian=TRUE)
Warning messages:
1: NaNs produced in: dweibull(x, shape, scale, log)
2: NA/Inf replaced by maximum positive value
3: NaNs produced in: dweibull(x, shape, scale...