{\beta}_p)^T\sim N((\beta_1,\cdots,\beta_p)^T,{\mbox{\boldmath $\Sigma$}}^*)$, where the $(k,l)$th element in ${\mbox{\boldmath $\Sigma$}}^*$ is $\displaystyle{\mbox{\boldmath $\Sigma$}}_{kl}^*=\sigma^2s_{kl}/(ns_{kk}s_{ll})$.]{}

For ease of notation, let $Z_1,\cdots,Z_p$ be the standardized random variables of $\widehat{\beta}_1,\cdots,\widehat{\beta}_p$, that is, $$\label{b2}
Z_i=\frac{\widehat{\beta}_i}{{\mbox{SD}}(\widehat{\beta}_i)}=\frac{\widehat{\beta}_i}{\sigma/(\sqrt{n}s_{ii})}, \quad\quad i=1,\cdots,p.$$ In the above, we implicitly assume that $\sigma$ is known and the above standardized random variables are z-test statistics. The estimate of residual variance $\sigma^2$ will be discussed in Section 6 via refitted cross-validation (Fan, Guo & Hao, 2011). Then, conditioning on $\{X_j^i\}$, $$\label{c1}
(Z_1,\cdots,Z_p)^T\sim N((\mu_1,\cdots,\mu_p)^T,{\mbox{\boldmath $\Sigma$}}),$$ where $\mu_i=\sqrt{n}\beta_is_{ii}/\sigma$ and covariance matrix ${\mbox{\boldmath $\Sigma$}}$ has the $(k,l)$th element as $s_{kl}$. Simultaneously testing (2) based on $(\widehat{\beta}_1,\cdots,\widehat{\beta}_p)^T$ is thus equivalent to testing

$$\label{c2}
H_{0j}:\quad\mu_j=0\quad\text{vs}\quad H_{1j}:\quad \mu_j\neq0, \quad\quad j=1,\cdots,p$$

based on $(Z_1,\cdots,Z_p)^T$.

In (4), ${\mbox{\boldmath $\Sigma$}}$ is the population covariance matrix of $(Z_1,\cdots,Z_p)^T$, and is known based on the sample data $\{X_j^i\}$. The covariance matrix ${\mbox{\boldmath $\Sigma$}}$ can have arbitrary dependence structure. We would
{ \beta}_p)^T\sim N((\beta_1,\cdots,\beta_p)^T,{\mbox{\boldmath $ \Sigma$}}^*)$, where the $ (k, l)$th element in $ { \mbox{\boldmath $ \Sigma$}}^*$ is $ \displaystyle{\mbox{\boldmath $ \Sigma$}}_{kl}^*=\sigma^2s_{kl}/(ns_{kk}s_{ll})$. ] { } 

 For ease of notation, let $ Z_1,\cdots, Z_p$ be the exchangeable random variable star of $ \widehat{\beta}_1,\cdots,\widehat{\beta}_p$, that is, $ $ \label{b2 } 
 Z_i=\frac{\widehat{\beta}_i}{{\mbox{SD}}(\widehat{\beta}_i)}=\frac{\widehat{\beta}_i}{\sigma/(\sqrt{n}s_{ii }) }, \quad\quad i=1,\cdots, p.$$ In the above, we implicitly assume that $ \sigma$ is known and the above exchangeable random variables are omega - test statistics. The appraisal of residual variance $ \sigma^2 $ will be discussed in Section 6 via refit crabbed - validation (Fan, Guo & Hao, 2011). Then, conditioning on $ \{X_j^i\}$, $ $ \label{c1 } 
 (Z_1,\cdots, Z_p)^T\sim N((\mu_1,\cdots,\mu_p)^T,{\mbox{\boldmath $ \Sigma$}}),$$ where $ \mu_i=\sqrt{n}\beta_is_{ii}/\sigma$ and covariance matrix $ { \mbox{\boldmath $ \Sigma$}}$ suffer the $ (k, l)$th element as $ s_{kl}$. Simultaneously examination (2) based on $ (\widehat{\beta}_1,\cdots,\widehat{\beta}_p)^T$ is thus equivalent to testing 

 $ $ \label{c2 } 
 H_{0j}:\quad\mu_j=0\quad\text{vs}\quad H_{1j}:\quad \mu_j\neq0, \quad\quad j=1,\cdots, p$$ 

 based on $ (Z_1,\cdots, Z_p)^T$. 

 In (4), $ { \mbox{\boldmath $ \Sigma$}}$ is the population covariance matrix of $ (Z_1,\cdots, Z_p)^T$, and is known based on the sample datum $ \{X_j^i\}$. The covariance matrix $ { \mbox{\boldmath $ \Sigma$}}$ can have arbitrary dependence structure. We would
{\betw}_p)^T\sim N((\beta_1,\cdots,\beta_p)^T,{\mnox{\boldmath $\Sigmc$}}^*)$, where the $(k,m)$th elemdnt in ${\mbox{\boldmath $\Sigma$}}^*$ is $\dusplatstyle{\mbox{\boldmath $\Sigoa$}}_{kl}^*=\sigma^2d_{kl}/(ns_{kk}s_{ol})$.]{}

Foc ease of notation, let $Z_1,\gbots,Z_l$ be chx standardized tandom variatles of $\widehad{\bdtc}_1,\cdots,\widehat{\beta}_p$, that is, $$\label{b2}
Z_i=\fwac{\widenah{\beta}_i}{{\mbox{SD}}(\wiqehau{\beea}_i)}=\fdac{\widehat{\beta}_i}{\sigma/(\sqrt{n}s_{ii})}, \quad\qhad i=1,\cdmts,p.$$ In the anove, we implicitly assume hhat $\sigma$ is known anf the above stageardized ranaom variables are z-tesj statistics. The estimate of reskdual variance $\witma^2$ fill be disrussed in Section 6 via refhtted ctoss-validation (Fai, Gui & Hao, 2011). Then, conditioiing on $\{X_j^i\}$, $$\label{c1}
(Z_1,\cqots,Z_p)^T\sik U((\mu_1,\cdots,\mu_p)^T,{\mbox{\boldmqty $\Sigka$}}),$$ wvere $\nu_i=\rqru{n}\bxta_js_{ii}/\sihma$ and covarjance matriz ${\mbox{\boldmath $\Sigms$}}$ rqs the $(k,l)$th emement af $s_{kl}$. Simultaneously testing (2) based on $(\fidshat{\beta}_1,\cdots,\widehat{\betq}_p)^T$ is thus equivalenj to testigg

$$\label{c2}
H_{0j}:\quad\mu_j=0\quad\text{vs}\quad H_{1j}:\quad \mu_j\neq0, \quag\quad g=1,\cditf,o$$

vaded on $(Z_1,\cdots,Z_p)^T$.

In (4), ${\mbox{\boldmath $\Sigma$}}$ is the lokulstion covariange matrix of $(Z_1,\cdotx,Z_o)^T$, wnd is known cased on fhe sample data $\{X_j^l\}$. The cjvariqnce matryx ${\mnox{\boldmath $\Sigma$}}$ can have qrbitrary deiendwnce structure. We could
{\beta}_p)^T\sim N((\beta_1,\cdots,\beta_p)^T,{\mbox{\boldmath $\Sigma$}}^*)$, where the $(k,l)$th element $\Sigma$}}^*$ $\displaystyle{\mbox{\boldmath $\Sigma$}}_{kl}^*=\sigma^2s_{kl}/(ns_{kk}s_{ll})$.]{} ease of notation, random of $\widehat{\beta}_1,\cdots,\widehat{\beta}_p$, that $$\label{b2} Z_i=\frac{\widehat{\beta}_i}{{\mbox{SD}}(\widehat{\beta}_i)}=\frac{\widehat{\beta}_i}{\sigma/(\sqrt{n}s_{ii})}, \quad\quad In the above, we implicitly assume $\sigma$ is known and the above standardized random variables are z-test statistics. The of residual variance $\sigma^2$ will be discussed in Section 6 via refitted cross-validation Guo Hao, Then, on $\{X_j^i\}$, $$\label{c1} (Z_1,\cdots,Z_p)^T\sim N((\mu_1,\cdots,\mu_p)^T,{\mbox{\boldmath $\Sigma$}}),$$ where $\mu_i=\sqrt{n}\beta_is_{ii}/\sigma$ and covariance matrix ${\mbox{\boldmath $\Sigma$}}$ has the $(k,l)$th as $s_{kl}$. Simultaneously testing (2) based on $(\widehat{\beta}_1,\cdots,\widehat{\beta}_p)^T$ thus equivalent to testing H_{0j}:\quad\mu_j=0\quad\text{vs}\quad H_{1j}:\quad \mu_j\neq0, \quad\quad j=1,\cdots,p$$ on In (4), $\Sigma$}}$ the covariance matrix of and is known based on the sample data $\{X_j^i\}$. The covariance matrix ${\mbox{\boldmath $\Sigma$}}$ can have arbitrary structure. We
{\beta}_p)^T\sim N((\beta_1,\cdots,\beta_p)^T,{\Mbox{\boldmaTh $\SigMa$}}^*)$, wHerE tHe $(k,l)$Th elEment in ${\mbox{\bolDMath $\sigma$}}^*$ is $\displaystyle{\mboX{\boldMaTH $\SigMA$}}_{kL}^*=\sigmA^2s_{kl}/(ns_{kK}S_{lL})$.]{}

fOr eAsE oF noTaTIoN, let $Z_1,\CdoTs,Z_p$ be tHe standardIzeD rAndom variablES oF $\widehat{\beTa}_1,\cDots,\widehat{\bEta}_P$, that iS, $$\lAbeL{B2}
Z_i=\frAc{\wIdehaT{\beta}_i}{{\MBox{SD}}(\wIdehat{\betA}_i)}=\FRac{\widEHat{\beta}_I}{\SIgMa/(\sqRt{n}s_{ii})}, \quad\quad i=1,\cdOTs,P.$$ in the above, we imPlicitLy ASsUME thAt $\sIgma$ is knowN aNd the ABove staNDaRDIZed RAndom variableS are z-test stATisTics. ThE eStiMAte of rEsiduAl VAriAnce $\sigma^2$ wiLl be Discussed In SectIOn 6 via reFItted crOss-valIdaTioN (Fan, gUo & haO, 2011). ThEn, COndITiOniNG on $\{x_j^i\}$, $$\label{C1}
(Z_1,\CdOts,Z_p)^t\sim n((\MU_1,\CDots,\Mu_p)^t,{\mboX{\boldMath $\Sigma$}}),$$ wherE $\mu_I=\sqrT{N}\beTa_is_{iI}/\sigmA$ and CoVariaNce matRix ${\mbOx{\Boldmath $\Sigma$}}$ haS the $(K,l)$th elemeNt aS $s_{Kl}$. SImUltanEOusly tEstIng (2) Based on $(\Widehat{\BEta}_1,\CdOTS,\WiDehat{\beta}_p)^T$ is thus eQuIVAlEnt to tesTing

$$\laBEl{C2}
H_{0J}:\Quad\mu_j=0\qUaD\teXt{vs}\QUAd H_{1j}:\qUad \mU_J\nEq0, \quad\quAd j=1,\cdoTS,p$$

BaSed on $(Z_1,\cDoTs,Z_p)^T$.

IN (4), ${\mBox{\BolDmath $\sIgma$}}$ Is the pOpulatioN covaRIance matrix of $(Z_1,\CDots,Z_p)^T$, and is kNOwN BAsED on tHe sAmple data $\{X_j^I\}$. The COvarIancE MaTriX ${\Mbox{\bOldmaTh $\sIgMA$}}$ can have arbitrary dePeNdence StrucTure. We would
{\beta}_p)^T\sim N((\beta_ 1,\cdots,\ beta_ p)^ T,{ \m box{ \bol dmath $\Sigma$ } }^*) $, where the $(k,l)$th  elem en t  in$ {\ mbox{ \boldma t h$ \ Sig ma $} }^* $i s$\dis pla ystyle{ \mbox{\bol dma th  $\Sigma$}}_ { kl }^*=\sigma ^2s _{kl}/(ns_{k k}s _{ll}) $. ]{} 
Foreas e ofnotati o n, let  $Z_1,\cd ot s ,Z_p$b e the s t a nd ardi zed random variab l es of $\widehat{\ beta}_ 1, \ cd o t s,\ wid ehat{\beta }_ p$, t h at is,$ $\ l a b el{ b 2}
Z_i=\frac{ \widehat{\b e ta} _i}{{\ mb ox{ S D}}(\w ideha t{ \ bet a}_i)}=\fra c{\w idehat{\b eta}_i } {\sigma / (\sqrt{ n}s_{i i}) },\qua d \q ua d i =1 , \cd o ts ,p. $ $ I n the ab ov e,  we i mpli c i t l y as sum e th at $\ sigma$ is kno wnandt heabove  stan dard iz ed ra ndom v ariab le s are z-test st atis tics. The  es ti mat eof re s idualvar ian ce $\si gma^2$w ill  b e d is cussed in Section6v i arefitted  cross - va li d ation (F an , G uo & H ao, 2 011) .  T hen, con dition i ng  o n $\{X_ j^ i\}$,$$ \la bel {c1}( Z_1, \cdots ,Z_p)^T\ sim N ( (\mu_1,\cdots, \ mu_p)^T,{\mbo x {\ b o ld m ath$\S igma$}}),$$  whe r e $\ mu_i = \s qrt { n}\be ta_is _{ i i} / \sigma$ and covaria nc e matr ix ${ \mbox{\boldma th $\Sigma $ } } $ has th e $( k ,l ) $th element as  $s_{ kl}$. Simu l taneousl y tes ting (2)  based on $ (\wideha t{\ bet a}_ 1,\ c d ot s,\widehat{\b e t a}_p )^ T$ is t hus  equiva len t t o t est in g

$$\lab el{c2}
H _{ 0j }: \q uad \mu_j = 0\quad\t ex t{v s} \qu ad H_ { 1j}:\q uad \ mu_j \n eq 0 , \ quad\qu a dj = 1,\c do ts ,p$$ 

b as ed on  $(Z _ 1,\ cdots,Z _p)^T$.
In( 4),${ \m box{\bo ldmath $\Sigm a$ }}$ is the  p opu lation c ovarianc e matrix of $(Z_1,\cdot s ,Z_p)^T $,and i s kn own based  on  the s amp l e data  $\{X_ j^i\} $.  Th e covar i a nc e m at rix ${\mbo x { \bo ldmat h$\Si gma$}}$  can have arbitrar y  de pendence stru ctu re.W e  w oul d 
{\beta}_p)^T\sim N((\beta_1,\cdots,\beta_p)^T,{\mbox{\boldmath_$\Sigma$}}^*)$, where_the $(k,l)$th element in_${\mbox{\boldmath $\Sigma$}}^*$_is_$\displaystyle{\mbox{\boldmath $\Sigma$}}_{kl}^*=\sigma^2s_{kl}/(ns_{kk}s_{ll})$.]{}

For_ease_of notation, let_$Z_1,\cdots,Z_p$ be the_standardized random variables of_$\widehat{\beta}_1,\cdots,\widehat{\beta}_p$, that is,_$$\label{b2}
Z_i=\frac{\widehat{\beta}_i}{{\mbox{SD}}(\widehat{\beta}_i)}=\frac{\widehat{\beta}_i}{\sigma/(\sqrt{n}s_{ii})},_\quad\quad i=1,\cdots,p.$$ In the above, we implicitly assume that $\sigma$ is known and the_above_standardized random_variables_are_z-test statistics. The estimate of_residual variance $\sigma^2$ will be_discussed in_Section 6 via refitted cross-validation (Fan, Guo &_Hao,_2011). Then, conditioning_on $\{X_j^i\}$, $$\label{c1}
(Z_1,\cdots,Z_p)^T\sim N((\mu_1,\cdots,\mu_p)^T,{\mbox{\boldmath $\Sigma$}}),$$ where $\mu_i=\sqrt{n}\beta_is_{ii}/\sigma$ and covariance_matrix ${\mbox{\boldmath $\Sigma$}}$ has the $(k,l)$th_element as $s_{kl}$._Simultaneously_testing_(2) based on $(\widehat{\beta}_1,\cdots,\widehat{\beta}_p)^T$_is thus equivalent to testing

$$\label{c2}
H_{0j}:\quad\mu_j=0\quad\text{vs}\quad H_{1j}:\quad_\mu_j\neq0, \quad\quad j=1,\cdots,p$$

based on $(Z_1,\cdots,Z_p)^T$.

In (4),_${\mbox{\boldmath $\Sigma$}}$ is the population covariance matrix_of $(Z_1,\cdots,Z_p)^T$, and is known based_on the sample data $\{X_j^i\}$._The covariance_matrix ${\mbox{\boldmath $\Sigma$}}$ can have_arbitrary dependence structure._We would