|\theta_{10j}|)=\gamma\neq 0$.

Implementation
--------------

Similar to Leng and Tang (2012), a nonlinear optimization procedure can be employed to maximize $\ell_p(\theta)$ given in Equation (\[PETL13\]). It is quite difficult to implement the nonlinear optimization procedure because of the nonconcave penalty function $p_\gamma(|\theta_j|)$ involved. To address the issue, we consider the following local quadratic approximation to the penalty function (Fan and Li, 2001) at a fixed value $\theta_j^{(m)}$ of $\theta_j$: $p_\gamma(|\theta_j|)\approx p_\gamma(|\theta_j^{(m)}|)+\frac{1}{2}\{p_\gamma '(|\theta_j^{(m)}|)/|\theta_j^{(m)}|\}\{\theta_j^2-(\theta_j^{(m)})^2\}$, where $\theta_j^{(m)}$ is the estimated value of $\theta_j$ at the $m$th step and $\theta_j$ is the $j$th component of $\theta$. Thus, the nonlinear optimization algorithm given in Owen (2001) can be adopted to maximize Equation (\[PETL13\]) based on the above local quadratic approximation of $p_\gamma(|\theta_j|)$. Repeating the nonlinear optimization procedure until convergence yields the PET estimate $\hat\theta$ of $\theta$.

To obtain the PET estimate of $\theta$, it is also necessary to find a data-driven approach to select the penalty parameter $\gamma$. To select an appropriate penalty parameter $\gamma$, we consider the following adjusted aBIC criterion: ${\rm aBIC}(\gamma)=-2\ell(\hat\theta_\gamma)+C_n\frac{\log(n)}{n}{\rm df}_\gamma$, where $\ell(\theta)=\ell(\nu,\theta)$ is given in Equation (\[PETL22\]), $\hat\theta_\gamma$ is the PET estimator of $\theta$ depending on the tuning parameter $\gamma$, df$_\gamma$ is the number of nonzero components in $\theta$ representing the “degrees of freedom" of the estimated unconditional moment models, $C_n$ is a scaling factor diverging to infinity at a slow rate as $p\rightarrow \infty$. When $p$ is fixed, we set $C_n=1$, otherwise we take $C_n={\rm max}\{\
|\theta_{10j}|)=\gamma\neq 0$. 

 Implementation 
 -------------- 

 Similar to Leng and Tang (2012), a nonlinear optimization procedure can be employed to maximize $ \ell_p(\theta)$ given in Equation (\[PETL13\ ]). It is quite unmanageable to enforce the nonlinear optimization procedure because of the nonconcave penalty affair $ p_\gamma(|\theta_j|)$ involve. To address the issue, we regard the trace local quadratic approximation to the punishment affair (Fan and Li, 2001) at a fixed value $ \theta_j^{(m)}$ of $ \theta_j$: $ p_\gamma(|\theta_j|)\approx p_\gamma(|\theta_j^{(m)}|)+\frac{1}{2}\{p_\gamma' (|\theta_j^{(m)}|)/|\theta_j^{(m)}|\}\{\theta_j^2-(\theta_j^{(m)})^2\}$, where $ \theta_j^{(m)}$ is the calculate value of $ \theta_j$ at the $ m$th step and $ \theta_j$ is the $ j$th component of $ \theta$. therefore, the nonlinear optimization algorithm given in Owen (2001) can be adopted to maximize Equation (\[PETL13\ ]) based on the above local quadratic approximation of $ p_\gamma(|\theta_j|)$. Repeating the nonlinear optimization procedure until convergence succumb the PET estimate $ \hat\theta$ of $ \theta$. 

 To obtain the PET estimate of $ \theta$, it is besides necessary to find a data - drive approach path to select the punishment parameter $ \gamma$. To select an appropriate penalty parameter $ \gamma$, we study the following adjusted aBIC criterion: $ { \rm aBIC}(\gamma)=-2\ell(\hat\theta_\gamma)+C_n\frac{\log(n)}{n}{\rm df}_\gamma$, where $ \ell(\theta)=\ell(\nu,\theta)$ is given in Equation (\[PETL22\ ]), $ \hat\theta_\gamma$ is the PET estimator of $ \theta$ depend on the tuning parameter $ \gamma$, df$_\gamma$ is the number of nonzero components in $ \theta$ representing the “ degrees of freedom " of the estimated unconditional moment models, $ C_n$ is a scaling component diverging to infinity at a dull rate as $ p\rightarrow \infty$. When $ p$ is fixed, we set $ C_n=1 $, otherwise we consider $ C_n={\rm max}\{\
|\theha_{10j}|)=\gamma\neq 0$.

Implementatiun
--------------

Similar to Leut and Vang (2012), a nonlinexr optimization procedure cai be empliyed to maximize $\ell_p(\tfeta)$ givej in Equqtioi (\[PETL13\]). It is quivs diffigblt tk impnxment the nonlikear optimisation proceduse bzcause of the nonconcave penalty funstion $p_\balma(|\theta_j|)$ invojved. Eo asdress the issue, we consider the fkllowinj local quadratoc approximation to the pejaltj function (Fan and Li, 2001) at a fuxed calue $\theta_j^{(o)}$ of $\theta_j$: $p_\gamma(|\thetz_j|)\approx p_\gamma(|\theta_j^{(m)}|)+\frac{1}{2}\{p_\gammx '(|\theca_j^{(m)}|)/|\theta_j^{(m)}|\}\{\jkwta_u^2-(\jheta_j^{(m)})^2\}$, wherx $\thetw_j^{(m)}$ is the estimated value pf $\theta_j$ at tme $m$ti stwp and $\theta_j$ is the $o$th component of $\theja$. Thus, tha uonlinear optimizatiob qlgorhthm givdb iv Osei (2001) dan be adkpted to mzximize Equqtion (\[PETL13\]) based on trv above local quadrwtyc approximation of $p_\gamma(|\theta_j|)$. Repeatpng fhe nonlinear optimizatuon procedure until clnvergencq yields the PET estimate $\hat\theta$ of $\theta$.

To obtdin tie PEU cfgumwte of $\theta$, it is also necessary to find a dwfa-crpven approach to felect the lejakjy parameter $\gxmma$. To sslect an appropriahe penajty pqrameter $\damms$, we consider the following adjusted aBPC ceiterion: ${\rm aBIC}(\galma)=-2\ell(\hat\thzta_\gamka)+C_n\ftac{\log(n)}{n}{\rm df}_\gamma$, wherz $\ell(\tgeta)=\ell(\nu,\thfta)$ is gibdn in Equation (\[PDTL22\]), $\hdt\theta_\gamma$ is the PET eseimator oh $\theca$ depenaing on thq tuning pwramebar $\gamma$, df$_\gamma$ id the ngmber of nlnzero components in $\theta$ reprxxenting the “cecrevs of frezdom" on the estimated unconditional moment modeus, $C_n$ is a scalinj factor divqrging to infhjity at a slmw rate ws $p\eighrarrow \kvfty$. When $p$ is fixed, we set $C_n=1$, otyerwise we take $C_n={\vm mab}\{\
|\theta_{10j}|)=\gamma\neq 0$. Implementation -------------- Similar to Leng (2012), nonlinear optimization can be employed Equation It is quite to implement the optimization procedure because of the nonconcave function $p_\gamma(|\theta_j|)$ involved. To address the issue, we consider the following local quadratic to the penalty function (Fan and Li, 2001) at a fixed value $\theta_j^{(m)}$ $\theta_j$: p_\gamma(|\theta_j^{(m)}|)+\frac{1}{2}\{p_\gamma where is the estimated value of $\theta_j$ at the $m$th step and $\theta_j$ is the $j$th component $\theta$. Thus, the nonlinear optimization algorithm given in (2001) can be adopted maximize Equation (\[PETL13\]) based on above quadratic approximation $p_\gamma(|\theta_j|)$. the optimization procedure until yields the PET estimate $\hat\theta$ of $\theta$. To obtain the PET estimate of $\theta$, it is also to find approach to the parameter To select an parameter $\gamma$, we consider the following ${\rm aBIC}(\gamma)=-2\ell(\hat\theta_\gamma)+C_n\frac{\log(n)}{n}{\rm df}_\gamma$, where $\ell(\theta)=\ell(\nu,\theta)$ is given in (\[PETL22\]), $\hat\theta_\gamma$ the PET estimator of $\theta$ depending the tuning parameter $\gamma$, df$_\gamma$ is the number nonzero components in $\theta$ representing the “degrees of freedom" of the estimated unconditional moment models, a scaling factor diverging infinity at a rate $p\rightarrow When is fixed, set $C_n=1$, otherwise we take $C_n={\rm max}\{\
|\theta_{10j}|)=\gamma\neq 0$.

ImplementatIon
--------------

Similar To LenG anD TaNg (2012), A nonLineAr optimization PRoceDure can be employed to maxImize $\ElL_P(\theTA)$ gIven iN EquatiON (\[Petl13\]). It Is QuIte DiFFiCult tO imPlement The nonlineAr oPtImization proCEdUre because Of tHe nonconcave PenAlty fuNcTioN $P_\gammA(|\thEta_j|)$ iNvolveD. to addrEss the issUe, WE consiDEr the foLLOwIng lOcal quadratic apprOXiMAtion to the penaLty funCtIOn (fAN anD Li, 2001) At a fixed vaLuE $\thetA_J^{(m)}$ of $\theTA_j$: $P_\GAMma(|\THeta_j|)\approx p_\gAmma(|\theta_j^{(m)}|)+\FRac{1}{2}\{P_\gamma '(|\ThEta_J^{(M)}|)/|\theta_J^{(m)}|\}\{\theTa_J^2-(\TheTa_j^{(m)})^2\}$, where $\thEta_j^{(M)}$ is the estImated VAlue of $\tHEta_j$ at tHe $m$th sTep And $\ThetA_J$ iS tHe $j$Th COmpONeNt oF $\TheTa$. Thus, thE nOnLineaR optIMIZAtioN alGoriThm giVen in Owen (2001) can bE adOpteD To mAximiZe EquAtioN (\[PeTL13\]) baSed on tHe aboVe Local quadratic aPproXimation oF $p_\gAmMa(|\tHeTa_j|)$. RePEating The NonLinear oPtimizaTIon PrOCEDuRe until convergence YiELDs The PET esTimate $\HAt\ThETa$ of $\thetA$.

TO obTain THE PET eStimATe Of $\theta$, iT is alsO NeCeSsary to FiNd a datA-dRivEn aPproaCH to sElect tHe penaltY paraMEter $\gamma$. To selECt an appropriaTE pENAlTY parAmeTer $\gamma$, we cOnsiDEr thE folLOwIng ADjustEd aBIc cRItERion: ${\rm aBIC}(\gamma)=-2\ell(\hAt\Theta_\gAmma)+C_N\frac{\log(n)}{n}{\rm dF}_\gamma$, wherE $\ELL(\theta)=\elL(\nu,\tHEtA)$ Is given in EquatIon (\[PEtL22\]), $\hat\theta_\GAmma$ is thE PET eStimator Of $\theta$ dePENding on tHe tUniNg pAraMETeR $\gamma$, df$_\gamma$ IS The nUmBer of noNzeRo compoNenTs iN $\thEta$ RePresentinG the “degrEeS oF fReEdoM" of thE EstimateD uNcoNdItiOnal mOMent moDels, $C_N$ is a ScAlINg fActor diVErGINg to InFiNity At a SlOw ratE as $p\RIghTarrow \iNfty$. When $p$ Is fIXed, wE sEt $c_n=1$, otherWise we take $C_n={\rM mAx}\{\
|\theta_{10j}|)=\gamma\neq  0$.

Impl ement ati on-- ---- ---- ----

Similart o Le ng and Tang (2012), anonli ne a r op t im izati on proc e du r e  ca nbe  em pl o ye d tomax imize $ \ell_p(\th eta )$  given in Eq u at ion (\[PET L13 \]). It is q uit e diff ic ult to im ple mentthe no n linear  optimiza ti o n proc e dure be c a us e of  the nonconcave p e na l ty function $p _\gamm a( | \t h e ta_ j|) $ involved .To ad d ress th e  i s s u e,w e consider th e following loc al qua dr ati c  appro ximat io n  to  the penalt y fu nction (F an and Li, 200 1 ) at afixedval ue$\th e ta _j ^{( m) } $ o f  $ \th e ta_ j$: $p_\ ga mm a(|\t heta _ j | ) \app rox  p_\ gamma (|\theta_j^{( m)} |)+\ f rac {1}{2 }\{p_ \gam ma  '(|\ theta_ j^{(m )} |)/|\theta_j^{( m)}| \}\{\thet a_j ^2 -(\ th eta_j ^ {(m)}) ^2\ }$,  where$\theta _ j^{ (m ) } $  i s the estimated va lu e of  $\theta _j$ at th e$ m$th ste pand  $\t h e ta_j$  ist he  $j$th c ompone n tof  $\thet a$ . Thus ,the  no nline a r op timiza tion alg orith m  given in Owen (2001) can be ad o p te d  tomax imize Equat ion( \[PE TL13 \ ])  ba s ed on  theab o ve local quadratic app ro ximati on of  $p_\gamma(|\ theta_j|)$ . R epeating  the no n linear optimiz ation  procedure until co nverg ence yie lds the P E T  estimat e $ \ha t\t het a $  o f $\theta$.
T o  obt ai n the P ETestimat e o f $ \th eta $,  it is al so neces sa ry  t ofin d a d a ta-drive napp ro ach  to s e lect t he pe nalt ypa r ame ter $\g a mm a $ . To  s el ectanap propr iate pen alty pa rameter $ \ga m ma$,  w econside r the followi ng  adjustedaB ICcriter i o n: ${\rm  aBIC}(\gamma)=-2\ell(\ h at\thet a_\ gamma )+C_ n\frac{\l og( n)}{n} {\r m  df}_\ gamma$ , whe re  $\ e l l(\th e t a) =\e ll (\nu,\thet a ) $ i s giv en  inEquatio n (\[PETL22\]), $\ h at\ theta_\gamma$  is  the P ET  es t im a tor  o f  $\ t h eta$ dependingon the tun in g  p arameter $ \ gam ma $, df$_ \gamma$  is t h e numbe r of nonz ero compo ne ntsi n  $\ theta$ rep resentin g the “de g reeso ffreed om"  of th eest imate d unco n dit ional  momen tmodels , $C_ n$  is a sc aling factor divergingto inf inity  at  a slow r ate as$p\righta rrow  \infty$.Whe n $ p$ is  fi x ed, w e se t  $ C_n = 1$, o ther w ise we ta k e$C_ n = {\ rm max}\{\
|\theta_{10j}|)=\gamma\neq 0$.

Implementation
--------------

Similar_to Leng_and Tang (2012), a_nonlinear optimization_procedure_can be_employed_to maximize $\ell_p(\theta)$_given in Equation_(\[PETL13\]). It is quite_difficult to implement_the_nonlinear optimization procedure because of the nonconcave penalty function $p_\gamma(|\theta_j|)$ involved. To address the_issue,_we consider_the_following_local quadratic approximation to the_penalty function (Fan and Li,_2001) at_a fixed value $\theta_j^{(m)}$ of $\theta_j$: $p_\gamma(|\theta_j|)\approx p_\gamma(|\theta_j^{(m)}|)+\frac{1}{2}\{p_\gamma_'(|\theta_j^{(m)}|)/|\theta_j^{(m)}|\}\{\theta_j^2-(\theta_j^{(m)})^2\}$,_where $\theta_j^{(m)}$ is_the estimated value of $\theta_j$ at the $m$th step_and $\theta_j$ is the $j$th component_of $\theta$. Thus,_the_nonlinear_optimization algorithm given in_Owen (2001) can be adopted to_maximize Equation (\[PETL13\]) based on the_above local quadratic approximation of $p_\gamma(|\theta_j|)$. Repeating_the nonlinear optimization procedure until convergence_yields the PET estimate $\hat\theta$_of $\theta$.

To_obtain the PET estimate of_$\theta$, it is_also necessary_to find a_data-driven approach to select the penalty_parameter $\gamma$. To_select an appropriate penalty parameter $\gamma$,_we_consider the following_adjusted_aBIC_criterion: ${\rm_aBIC}(\gamma)=-2\ell(\hat\theta_\gamma)+C_n\frac{\log(n)}{n}{\rm df}_\gamma$, where_$\ell(\theta)=\ell(\nu,\theta)$_is given_in_Equation (\[PETL22\]), $\hat\theta_\gamma$ is the PET_estimator_of $\theta$ depending on the tuning parameter_$\gamma$, df$_\gamma$ is the_number_of nonzero components in_$\theta$ representing the “degrees of_freedom" of the estimated unconditional moment_models, $C_n$_is a_scaling factor diverging to infinity at a slow rate as $p\rightarrow_\infty$. When $p$ is fixed, we_set $C_n=1$, otherwise we_take $C_n={\rm_max}\{\