On 19/01/14 06:20, peter dalgaard wrote:
Offhand, I'd say that the issue isn't so much the starting values as the risk of
hitting the upper bound on gamma. If you get gamma set to min(x) at any point during
the iteration, the density at the corresponding point is zero, the log-likelihood is
-Inf, and <poof>, the magic smoke escapes.
Perhaps try min(x)-1e-7 or so?
Yep, that did it. Thanks Peter.
(Actually I used min(x) - eps where I set eps equal to
sqrt(.Machine$double.eps) --- which seems to be a kind of standard
construction for a small-but-not-too-small number.)
cheers,
Rolf
______________________________________________
R-help@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.