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,
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
(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
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
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
5 matches
Mail list logo