 same as that of $\hat{\beta}_n$. These arguments can be justified similarly as those in @John:Stra:indu:2009.

Variance Estimation {#sect:vari}
===================

The asymptotic variance of $\tilde\beta_n$ is even harder to estimate for case-cohort data than for full cohort data because of the extra complexity caused by the data structure. The terms in the summation in $\tilde{U}_n(\beta)$ are not independent since the sub-cohort is drawn from the full cohort without replacement. We propose four variance estimators; one is fully resampling based while the other three use resampling to a component of the sandwich variance estimator.

Multiplier Bootstrap
--------------------

The multiplier bootstrap estimator of @Jin:Lin:Wei:Ying:rank:2003 is adapted to case-cohort data by inserting proper case-cohort weights, $h_i$’s, in the multiplier bootstrap estimating equations. Let $\eta_i$, $i = 1, \dots, n$, be independent and identically distributed positive random variables with ${\mathrm{E}}(\eta_i) = {\mathrm{Var}}(\eta_i) = 1$. Define $$\tilde{U}_n^{c*}(\beta) = \sum_{i=1}^n \sum_{j=1}^n \eta_i\eta_jh_j\Delta_i (X_i- X_j) \Phi \left[ \frac{e_j(\beta)-e_i(\beta)}{r_{ij}^2} \right].
  \label{eq:UnSmoothC:Res}$$ For a realization of $(\eta_1, \ldots, \eta_n)$, the solution to  provides one draw of $\tilde\beta_n$ from its asymptotic distribution. By repeating this process a large number $B$ times, the variance matrix of $\tilde\beta_n$ can be estimated directly by the sampling variance matrix of the bootstrap sample of $\tilde\beta_n$.

Since the asymptotic variance of $\hat{\beta}_n$ is the same as that of $\tilde{\beta}_n$, the covariance matrix of $\tilde{\beta}_n$ can also be estimated by  through multiplier bootstrap. This is, however, not recommended because it would need to solve a large number $B$ nonsmooth estimating equations. As will be seen in our simulation study
same as that of $ \hat{\beta}_n$. These arguments can be justified similarly as those in @John: Stra: indu:2009. 

 Variance Estimation { # sect: vari } 
 = = = = = = = = = = = = = = = = = = = 

 The asymptotic discrepancy of $ \tilde\beta_n$ is even hard to estimate for case - cohort datum than for full cohort data because of the excess complexity induce by the data social organization. The term in the summation in $ \tilde{U}_n(\beta)$ are not independent since the sub - cohort is drawn from the full age group without replacement. We propose four variance calculator; one is fully resampling based while the other three use resampling to a component of the sandwich discrepancy estimator. 

 Multiplier Bootstrap 
 -------------------- 

 The multiplier bootstrap estimator of @Jin: Lin: Wei: Ying: rank:2003 is adapted to case - cohort datum by inserting proper case - cohort weights, $ h_i$ ’s, in the multiplier bootstrap estimating equations. lease $ \eta_i$, $ i = 1, \dots, n$, be independent and identically distributed positive random variables with $ { \mathrm{E}}(\eta_i) = { \mathrm{Var}}(\eta_i) = 1$. Define $ $ \tilde{U}_n^{c*}(\beta) = \sum_{i=1}^n \sum_{j=1}^n \eta_i\eta_jh_j\Delta_i (X_i- X_j) \Phi \left [ \frac{e_j(\beta)-e_i(\beta)}{r_{ij}^2 } \right ]. 
   \label{eq: UnSmoothC: Res}$$ For a realization of $ (\eta_1, \ldots, \eta_n)$, the solution to    provides one draw of $ \tilde\beta_n$ from its asymptotic distribution. By repeating this process a large act $ B$ times, the variance matrix of $ \tilde\beta_n$ can be calculate immediately by the sample variance matrix of the bootstrap sample of $ \tilde\beta_n$. 

 Since the asymptotic variance of $ \hat{\beta}_n$ is the same as that of $ \tilde{\beta}_n$, the covariance matrix of $ \tilde{\beta}_n$ can also be estimated by    through multiplier bootstrap. This is, however, not recommended because it would need to resolve a large number $ B$ nonsmooth estimating equations. As will be seen in our simulation discipline
 sale as that of $\hat{\beta}_n$. Tmese arguments ccb be jnstifies similafly as those in @John:Stra:indu:2009.

Tariqnce Tftimation {#sect:vari}
===================

Ghe asympnotic varuanct of $\tilde\beta_n$ is even harder fl escinate for case-cphort data than for full cuhlrt data because of the extra complqxity csuded by the datw stgusturs. The terms in the summation in $\timde{U}_n(\beua)$ are not indepencent since the sub-cohort id drwwn from the full fohort withiut wwplacement. Wd propose four variancg estimators; one is fully resampuing yased while tye ljher three uwe refampling to a componant of yhe sandwich vsrience estimator.

Multiplier Uootstrap
--------------------

The multiplyer bootsdrcp estimator of @Jin:Lib:Wwi:Yinc:rann:2003 is qdaotes vo dase-cojorv data by ihserting priper case-cohort weibhew, $h_i$’s, in the jultipjiqr bootstrap estimating equations. Let $\eua_i$, $i = 1, \dots, n$, be independenr and identically disjributed pjsitive random variables with ${\mathrm{E}}(\eta_i) = {\mathrm{Ear}}(\ete_i) = 1$. Bcnine $$\ripde{U}_n^{c*}(\beta) = \sum_{i=1}^n \sum_{j=1}^n \eta_i\eta_jh_j\Delta_i (X_i- X_t) \Pni \left[ \frac{e_j(\beba)-e_i(\beta)}{r_{ij}^2} \right].
  \lwbrj{eq:UnSmoothC:Rgs}$$ For c rsalization of $(\eta_1, \pdots, \eja_n)$, thw solutiog to  provides one draw of $\tilde\veta_n$ from ins awymptotic distribucion. By repecting jhis ptocess a large number $B$ timss, the variwnce matrjb of $\tilde\beta_n$ zan ba estimaudd directly by thq samplinj varnance magrix of thq bootstrao samine of $\tilde\beta_n$.

Sijce tke asfmptotic vwriance of $\hat{\beta}_n$ is the same as that of $\tikda{\bena}_n$, the cjvarisnce matrix os $\tilde{\beta}_n$ ccn also ye estkmated by  nhrough mnltiplier bojtstrap. This hd, however, nov recommegded becquse it dould need to xolve a lcxge numbee $B$ nonsmooth estikatkhg equations. As wiol be seen in oir rimtlwtmon segdy
same as that of $\hat{\beta}_n$. These arguments justified as those @John:Stra:indu:2009. Variance Estimation of is even harder estimate for case-cohort than for full cohort data because the extra complexity caused by the data structure. The terms in the summation $\tilde{U}_n(\beta)$ are not independent since the sub-cohort is drawn from the full cohort replacement. propose variance one is fully resampling based while the other three use resampling to a component of the variance estimator. Multiplier Bootstrap -------------------- The multiplier bootstrap of @Jin:Lin:Wei:Ying:rank:2003 is adapted case-cohort data by inserting proper weights, in the bootstrap equations. $\eta_i$, $i = \dots, n$, be independent and identically distributed positive random variables with ${\mathrm{E}}(\eta_i) = {\mathrm{Var}}(\eta_i) = 1$. Define = \sum_{i=1}^n (X_i- X_j) \left[ \right]. For a realization \ldots, \eta_n)$, the solution to provides $\tilde\beta_n$ from its asymptotic distribution. By repeating this a large $B$ times, the variance matrix of can be estimated directly by the sampling variance of the bootstrap sample of $\tilde\beta_n$. Since the asymptotic variance of $\hat{\beta}_n$ is the same of $\tilde{\beta}_n$, the covariance of $\tilde{\beta}_n$ can be by multiplier This is, not recommended because it would need to solve a large number nonsmooth estimating equations. As will be seen in our simulation
 same as that of $\hat{\beta}_n$. These Arguments cAn be jUstIfiEd SimiLarlY as those in @John:sTra:iNdu:2009.

Variance Estimation {#sEct:vaRi}
===================

tHe asYMpTotic VariancE Of $\TILde\BeTa_N$ is EvEN hArder To eStimate For case-cohOrt DaTa than for fulL CoHort data beCauSe of the extra ComPlexitY cAusED by thE daTa strUcture. tHe termS in the sumMaTIon in $\tILde{U}_n(\beTA)$ ArE not Independent since tHE sUB-cohort is drawn From thE fULl COHorT wiThout replaCeMent. WE Propose FOuR VARiaNCe estimators; oNe is fully reSAmpLing baSeD whILe the oTher tHrEE usE resampling To a cOmponent oF the saNDwich vaRIance esTimatoR.

MuLtiPlieR boOtStrAp
--------------------

tHe mULtIplIEr bOotstrap EsTiMator Of @JiN:lIN:wei:YIng:Rank:2003 Is adaPted to case-cohOrt Data BY inSertiNg proPer cAsE-cohoRt weigHts, $h_i$’S, iN the multiplier bOotsTrap estimAtiNg EquAtIons. LET $\eta_i$, $i = 1, \DotS, n$, bE indepeNdent anD IdeNtICALlY distributed positiVe RANdOm variabLes witH ${\MaThRM{E}}(\eta_i) = {\maThRm{VAr}}(\etA_I) = 1$. definE $$\tilDE{U}_N^{c*}(\beta) = \suM_{i=1}^n \sum_{J=1}^N \eTa_I\eta_jh_j\deLta_i (X_i- x_j) \phi \LefT[ \frac{E_J(\betA)-e_i(\betA)}{r_{ij}^2} \righT].
  \labeL{Eq:UnSmoothC:Res}$$ fOr a realizatioN Of $(\ETA_1, \lDOts, \eTa_n)$, The solution To  prOVideS one DRaW of $\TIlde\bEta_n$ fRoM ItS Asymptotic distributIoN. By repEatinG this process a Large numbeR $b$ TImes, the vAriaNCe MAtrix of $\tilde\beTa_n$ caN be estimatED directlY by thE samplinG variance MATrix of thE boOtsTraP saMPLe Of $\tilde\beta_n$.

SINCe thE aSymptotIc vAriance Of $\hAt{\bEta}_N$ is ThE same as thAt of $\tildE{\bEtA}_n$, ThE coVariaNCe matrix Of $\TilDe{\BetA}_n$ can ALso be eStimaTed bY  tHrOUgh MultiplIEr BOOtstRaP. THis iS, hoWeVer, noT recOMmeNded becAuse it wouLd nEEd to SoLvE a large Number $B$ nonsmoOtH estimatinG eQuaTions. AS WIll be seeN in our simulation study
 same as that of $\hat{\be ta}_n$. Th ese a rgu men ts  can  bejustified simi l arly  as those in @John:Str a:ind u: 2 009. 
V arian ce Esti m at i o n { #s ec t:v ar i }===== === ======= ====

Theasy mp totic varian c eof $\tilde \be ta_n$ is eve n h arderto  es t imate  fo r cas e-coho r t data  than for  f u ll coh o rt data b ec ause  of the extra com p le x ity caused bythe da ta st r u ctu re.  The terms  i n the summati o ni n $\t i lde{U}_n(\bet a)$ are not ind epende nt  si n ce the  sub- co h ort  is drawn f romthe fullcohort without replace ment.Wepro pose fo ur  va ri a nce es tim a tor s; one i sfu lly r esam p l i n g ba sed  whi le th e other three  us e re s amp lingto acomp on ent o f thesandw ic h variance esti mato r.

Multi pli er  Bo ot strap ------ --- --- ------- -

Them ult ip l i e rbootstrap estimato ro f  @ Jin:Lin: Wei:Yi n g: ra n k:2003 i sada pted t o cas e-co h or t data b y inse r ti ng  proper  c ase-co ho rtwei ghts, $h_i $’s, i n the mu ltipl i er bootstrap e s timating equa t io n s .L et $ \et a_i$, $i =1, \ d ots,  n$, be  in d epend ent a nd id e ntically distribute dpositi ve ra ndom variable s with ${\ m a t hrm{E}}( \eta _ i) = {\mathrm{Var }}(\e ta_i) = 1$ .  Define$$\ti lde{U}_n ^{c*}(\be t a ) = \sum _{i =1} ^n\su m _ {j =1}^n \eta_i\ e t a_jh _j \Delta_ i ( X_i- X_ j)\Ph i \ lef t[  \frac{e_ j(\beta) -e _i (\ be ta) }{r_{ i j}^2} \r ig ht] .  \ label { eq:UnS mooth C:Re s} $$ For  a real i za t i on o f$( \eta _1,  \ ldots , \e t a_n )$, the  solution  to  pro vi de s one d raw of $\tild e\ beta_n$ fr om  it s asym p t otic dis tribution. By repeating this pr oce ss alarg e number$B$  times , t h e vari ance m atrix  o f $ \ t ilde\ b e ta _n$  c an be esti m a ted  dire ct ly b y the s ampling variance m a tri x of the boot str ap s a m pl e o f  $ \ til de \ bet a _ n$.

Since theasymptotic  v a ri ance of $\ h at{ \b eta}_n$  is the  same as that  of $\til de{\beta} _n $, t h e  co variance m atrix of  $\tilde{ \ beta} _ n$  canals o be e st ima ted b y  thr o ugh  mult iplier  b ootstr ap. T hi s is, ho wever, not recommendedbecaus e itwou ld need t o s o lve  a largenumb er $B$ non smo oth  esti mat i ng eq uati o ns . A s  will  bes een in ou r  s imu l a ti on study
 same_as that_of $\hat{\beta}_n$. These arguments_can be_justified_similarly as_those_in @John:Stra:indu:2009.

Variance Estimation_{#sect:vari}
===================

The asymptotic variance_of $\tilde\beta_n$ is even_harder to estimate_for_case-cohort data than for full cohort data because of the extra complexity caused by_the_data structure._The_terms_in the summation in $\tilde{U}_n(\beta)$_are not independent since the_sub-cohort is_drawn from the full cohort without replacement. We_propose_four variance estimators;_one is fully resampling based while the other three_use resampling to a component of_the sandwich variance_estimator.

Multiplier_Bootstrap
--------------------

The_multiplier bootstrap estimator of_@Jin:Lin:Wei:Ying:rank:2003 is adapted to case-cohort data_by inserting proper case-cohort weights, $h_i$’s,_in the multiplier bootstrap estimating equations. Let_$\eta_i$, $i = 1, \dots, n$,_be independent and identically distributed_positive random_variables with ${\mathrm{E}}(\eta_i) = {\mathrm{Var}}(\eta_i)_= 1$. Define_$$\tilde{U}_n^{c*}(\beta) =_\sum_{i=1}^n \sum_{j=1}^n \eta_i\eta_jh_j\Delta_i_(X_i- X_j) \Phi \left[ \frac{e_j(\beta)-e_i(\beta)}{r_{ij}^2} \right].
_ \label{eq:UnSmoothC:Res}$$ For_a realization of $(\eta_1, \ldots, \eta_n)$,_the_solution to  provides_one_draw_of $\tilde\beta_n$_from its asymptotic_distribution._By repeating_this_process a large number $B$ times,_the_variance matrix of $\tilde\beta_n$ can be estimated_directly by the sampling_variance_matrix of the bootstrap_sample of $\tilde\beta_n$.

Since the asymptotic_variance of $\hat{\beta}_n$ is the same_as that_of $\tilde{\beta}_n$,_the covariance matrix of $\tilde{\beta}_n$ can also be estimated by  through_multiplier bootstrap. This is, however, not_recommended because it would_need to_solve_a large number_$B$_nonsmooth estimating_equations. As will be seen in our_simulation study