Re: [R] 回复: cch() and coxph() for case-coh ort

2008-07-16 Thread Jin Wang
How can we implement Prentice jacknife variance estimation in R strictly? Thanks 2008/6/17 Thomas Lumley <[EMAIL PROTECTED]>: > On Mon, 16 Jun 2008, Peter Dalgaard wrote: > > Jin Wang wrote: >> >>> I tried to compare if cch() and coxph() can generate same result

[R] implementation of Prentice method in cch()

2008-07-15 Thread Jin Wang
Case cohort function cch() is in survival package. In cch(), the prentice method is implemented like this: Prentice <- function(tenter, texit, cc, id, X, ntot,robust){ eps <- 0.0001 cens <- as.numeric(cc>0) # Censorship indicators subcoh <- as.numeric(cc<2) # Subcohort indicators

[R] How to run coxph() with time dependent cohort sampling

2008-06-30 Thread Jin Wang
Now that we have case cohort model , we have 1000 people and 50 cases Let the first 10 cases occur at the same time second 10 " third 10 " fourth 10" fifth10" How easy is it to randomly sample 50 di

[R] 回复: cch() and coxph() for case-coh ort

2008-06-16 Thread Jin Wang
eIII 0.597 1.818 1.285 2.571 0.001 stageIV 1.392 4.021 2.670 6.057 0.000 histolUH 1.506 4.507 3.274 6.203 0.000 age 0.043 1.044 0.996 1.095 0.069 2008/6/12, Terry Therneau <[EMAIL PROTECTED]>: > Jin Wang had an error. My original note specified a variable that was 1 for > subje

Re: [R] cch function and time dependent covariates

2008-06-12 Thread Jin Wang
1.allison.2 <- cch(Surv(start, stop, arrest.time) ~ + fin + age + race + wexp + mar + paro + prio + employed, + data=ccoh1.data.fold,subcoh=~subcohort,id=~seqno,cohort.size=19809) Error in cch(Surv(start, stop, arrest.time) ~ fin + age + race + wexp + : Multiple records per id not allow

Re: [R] cch function and time dependent covariates

2008-06-12 Thread Jin Wang
I tried your alternative method on the example in cch() description manual. The example data "nwtco" has not time-dependent covariates yet. I test cch() and coxph() on the same data. But the estimation result is different. I don't know if I did anything wrong. subcoh <- nwtco$in.subcohort selccoh

[R] I wonder if cch function in Survival package can calculate time dependent covariate

2008-06-11 Thread Jin Wang
ity or not? The cch() manual does not say if time dependent covariate is allowed I know coxph() in Survival package can estimate time dependent covariates. Thanks Jin Wang [[alternative HTML version deleted]] __ R-help@r-project.org mailing