[R] WinBUGS on survival, simple but confusing question

2011-06-07 Thread tingtingzhan
Hi All, I'm using WinBUGS on a very simple survival model (log-normal with one covariate "Treat"), but I cannot understand the way it handles censored data. I'm posting the R file which generates the data from pre-specified parameters, as well as the .bug file. The question is, if I use NA to de

Re: [R] how to store lme/lmer fit result

2011-03-03 Thread tingtingzhan
Hi All, I'm experiencing difficulties in saving a model fit by gls(). Basically if I just save() a gls object "gls.fit" in my workspace into an .RData file and later reload this .RData file, I get error when running script such as summary(gls.fit)$coefficients Dr. Bates's reply (quoted below)