" I am trying to fit a non-linear model for a parasite dataset.
Initially, I ...
The survreg() function will fit a weibull distribution. It has the
further ability to include censored data, a facility you don't need.
library(survival)
fit <- survreg(Surv(coccidiaopg) ~ age + year + sex*season, data=...
One day I should update it so that if there is no censoring then the
Surv() wrapper is not required. Also, see the examples in the ?survreg
documentation for comments on shape and scale. This routine uses a
location-scale parameterization which is different than rweibull.
Terry Therneau
______________________________________________
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.
