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.