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
_______________________________________________
R-sig-ecology mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-sig-ecology
