Re: [R] case-cohort

From: Terry Therneau <therneau_at_mayo.edu>
Date: Thu, 12 Jun 2008 10:27:07 -0500 (CDT)


Jin Wang had an error. My original note specified a variable that was 1 for subjects NOT in the subcohort, so the correct coxph call is

  coxph(Surv(edrel, rel) ~ stage + histol + age +

           offset(-100*(subcohort==0)) + cluster(seqno), data =ccoh.data)            

This gives the same coefficients as the cch example, along with the infinitesimal jackknife or "robust" variance estimate.

        Terry Therneau



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 Thu 12 Jun 2008 - 16:05:06 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 16 Jun 2008 - 20:32:15 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.

list of date sections of archive