\HeaderA{cch}{Fits proportional hazards regression model to case-cohort data}{cch}
\keyword{survival}{cch}
\begin{Description}\relax
Returns estimates and standard errors from relative risk
regression fit to data from case-cohort studies. A choice is available
among the Prentice, Self-Prentice and Lin-Ying methods for unstratified
data. For stratified data the choice is between Borgan I, a
generalization of the Self-Prentice estimator for unstratified
case-cohort data, and Borgan II, a generalization of the Lin-Ying
estimator.
\end{Description}
\begin{Usage}
\begin{verbatim}
cch(formula, data = sys.parent(), subcoh, id, stratum=NULL, cohort.size,
    method =c("Prentice","SelfPrentice","LinYing","I.Borgan","II.Borgan"),
    robust=FALSE)
\end{verbatim}
\end{Usage}
\begin{Arguments}
\begin{ldescription}
\item[\code{formula}] A formula object that must have a \code{\LinkA{Surv}{Surv}} object as the response. 
The Surv object must be of type \code{"right"}, or of type \code{"counting"}.

\item[\code{subcoh}] Vector of indicatorsfor subjects sampled as part of the
sub-cohort. Code \code{1} or \code{TRUE} for members of the
sub-cohort, \code{0} or \code{FALSE} for others. If \code{data} is a
data frame then \code{subcoh} may be a one-sided formula.

\item[\code{id}] Vector of unique identifiers, or formula specifying such a vector.

\item[\code{stratum}] A vector of stratum indicators or a formula specifying
such a vector
\item[\code{cohort.size}] Vector with size of each stratum original cohort from which subcohort was sampled

\item[\code{data}] An optional data frame in which to interpret the variables 
occurring in the formula. 

\item[\code{method}] Three procedures are available. The default method is "Prentice", with 
options for "SelfPrentice" or "LinYing". 
\item[\code{robust}] For \code{"LinYing"} only, if \code{robust=TRUE}, use design-based standard errors even for
phase I
\end{ldescription}
\end{Arguments}
\begin{Details}\relax
Implements methods for case-cohort data analysis described by Therneau and
Li (1999). The three methods differ in the choice of "risk sets" used to
compare the covariate values of the failure with those of others at risk at
the time of failure. "Prentice" uses the sub-cohort members "at risk" plus
the failure if that occurs outside the sub-cohort and is score unbiased.
"SelfPren" (Self-Prentice) uses just the sub-cohort members "at risk". These
two have the same asymptotic variance-covariance matrix. "LinYing" (Lin-Ying)
uses the all members of the sub-cohort and all failures outside the sub-cohort
who are "at risk". The methods also differ in the weights given to different
score contributions.

The \code{data} argument must not have missing values for any variables
in the model.  There must not be any censored observations outside the subcohort.
\end{Details}
\begin{Value}
An object of class "cch"  incorporating a list of estimated regression coefficients and two estimates of their 
asymptotic variance-covariance matrix.
\begin{ldescription}
\item[\code{coef}] regression coefficients.


\item[\code{naive.var}] Self-Prentice model based variance-covariance matrix.


\item[\code{var}] Lin-Ying empirical variance-covariance matrix. 


\end{ldescription}
\end{Value}
\begin{Author}\relax
Norman Breslow, modified by Thomas Lumley
\end{Author}
\begin{References}\relax
Prentice, RL (1986). A case-cohort design for epidemiologic cohort studies and
disease prevention trials. Biometrika 73: 1--11.


Self, S and Prentice, RL (1988). Asymptotic distribution theory and efficiency
results for case-cohort studies. Annals of Statistics 16: 64--81.


Lin, DY and Ying, Z (1993). Cox regression with incomplete covariate measurements.
Journal of the American Statistical Association 88: 1341--1349.


Barlow, WE (1994). Robust variance estimation for the case-cohort design. Biometrics
50: 1064--1072


Therneau, TM and Li, H (1999). Computing the Cox model for case-cohort designs.
Lifetime Data Analysis 5: 99--112.


Borgan, \eqn{\O}{O}, Langholz, B, Samuelsen, SO, Goldstein, L and Pogoda, J (2000)
Exposure stratified case-cohort designs. Lifetime Data Analysis 6, 39-58.
\end{References}
\begin{SeeAlso}\relax
\code{twophase} and \code{svycoxph} in the "survey" package for
more general two-phase designs. \url{http://faculty.washington.edu/tlumley/survey/}
\end{SeeAlso}
\begin{Examples}
\begin{ExampleCode}
## The complete Wilms Tumor Data 
## (Breslow and Chatterjee, Applied Statistics, 1999)
## subcohort selected by simple random sampling.
##

subcoh <- nwtco$in.subcohort
selccoh <- with(nwtco, rel==1|subcoh==1)
ccoh.data <- nwtco[selccoh,]
ccoh.data$subcohort <- subcoh[selccoh]
## central-lab histology 
ccoh.data$histol <- factor(ccoh.data$histol,labels=c("FH","UH"))
## tumour stage
ccoh.data$stage <- factor(ccoh.data$stage,labels=c("I","II","III","IV"))
ccoh.data$age <- ccoh.data$age/12 # Age in years

##
## Standard case-cohort analysis: simple random subcohort 
##

fit.ccP <- cch(Surv(edrel, rel) ~ stage + histol + age, data =ccoh.data,
   subcoh = ~subcohort, id=~seqno, cohort.size=4028)

fit.ccP

fit.ccSP <- cch(Surv(edrel, rel) ~ stage + histol + age, data =ccoh.data,
   subcoh = ~subcohort, id=~seqno, cohort.size=4028, method="SelfPren")

summary(fit.ccSP)

##
## (post-)stratified on instit
##
stratsizes<-table(nwtco$instit)
fit.BI<- cch(Surv(edrel, rel) ~ stage + histol + age, data =ccoh.data,
   subcoh = ~subcohort, id=~seqno, stratum=~instit, cohort.size=stratsizes,
   method="I.Borgan")

summary(fit.BI)
\end{ExampleCode}
\end{Examples}

