Re: [R] Direct adjusted survival?

2008-01-31 Thread Terry Therneau
The lines that I hoped to be the survival probabilities for each edtrt-group adjusted for confounding by log(bili) are nearly identical to the KM-lines, and they certainly don't appear adjusted for the very strong confounding by log(bili). I'm not quite sure what they are, though. Yes,

Re: [R] Direct adjusted survival?

2008-01-31 Thread Peter Jepsen
Dear Terry Therneau: Thank you for replying. Please forgive me for carrying on. Here is the example I gave, now with some output shown: require(survival) data(pbc) coxph(Surv(time,status)~edtrt,data=pbc) --- OUTPUT Call: coxph(formula = Surv(time, status) ~ edtrt, data = pbc) coef

Re: [R] Direct adjusted survival?

2008-01-30 Thread Peter Jepsen
(avg.fits) - list(NULL,var.values) list(time=bfit$time,fits=avg.fits) } *** -Oprindelig meddelelse- Fra: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] På vegne af Thomas Lumley Sendt: 30. januar 2008 00:31 Til: Peter Jepsen Cc: r-help Emne: Re: [R] Direct adjusted survival? On Wed, 30 Jan

[R] Direct adjusted survival?

2008-01-29 Thread Peter Jepsen
Hello, I am trying to find an R function to compute 'direct adjusted survival' with standard errors. A SAS-macro to do this is presented in Zhang X, Loberiza FR, Klein JP, Zhang MJ. A SAS macro for estimation of direct adjusted survival curves based on a stratified Cox regression model. Comput

Re: [R] Direct adjusted survival?

2008-01-29 Thread Thomas Lumley
On Wed, 30 Jan 2008, Peter Jepsen wrote: I am trying to find an R function to compute 'direct adjusted survival' with standard errors. A SAS-macro to do this is presented in Zhang X, Loberiza FR, Klein JP, Zhang MJ. A SAS macro for estimation of direct adjusted survival curves based on a