[R] addreg help

From: Björn Holmberg <fvoman_at_hotmail.com>
Date: Mon, 25 Feb 2008 09:06:26 +0000

Hello!
Im having a bit of a problem creating martingale-residuals for my aalen additive model. Is it even possible in R to do that? My code for the additive model is the following: ruff <-addreg(Surv(E$stop,E$delta) ~E$Zh +E$Zr +E$PSA +E$Treatment) summary(ruff)plot(ruff)  

If it had been a cox model i would have just written  

stage.ph <- coxph( Surv(E$stop,E$delta) ~E$Zh +E$Zr +E$Treatment, method="breslow", na.action=na.exclude) resid <- resid(stage.ph, type="martingale")str(resid)summary(resid) plot(E$PSA, E$resid, xlab="stage",ylab="Martingale Residuals", ylim=c(-4, 1.0))lines(lowess(E$PSA[!is.na(E$PSA) & !is.na(E$resid)], E$resid[!is.na(E$PSA) & !is.na(E$resid)])) So, how do I do this for the additive model?  

/Eric


        [[alternative HTML version deleted]]



R-help_at_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. Received on Mon 25 Feb 2008 - 09:42:05 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Mon 25 Feb 2008 - 10:30:16 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive