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