[R-sig-eco] lognormal gamlss with censored data
Lee, Laura
laura.lee at ncdenr.gov
Thu Jul 3 21:34:11 CEST 2014
Hello all!
I am fitting censored data with the 'gamlss' and 'survival' packages and assuming a lognormal distribution. Do I need to transform or back-transform the fitted values?
cy <- with(data, ifelse(NUMBERS>=Limit, Limit, NUMBERS))
ci <- with(data, ifelse(NUMBERS>=Limit, 0, 1))
data <- data.frame(data, cy, ci)
rm(cy,ci)
gen.cens("LOGNO",type="right")
cfit <- gamlss(Surv(cy, ci) ~ fyear + newWater, data=data, family=LOGNOrc)
Thanks for your help,
Laura
More information about the R-sig-ecology
mailing list