[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