m) 
&= \frac{e^{-\psi (\cW_{r_0})} \psi(\cW_{r_0})^m}{m!}
\prod_{i=1}^m \frac{{{\mathrm d}}\psi(r_{(i)}, {{\boldsymbol w}}_{(i)})}{\psi(\cW_{r_0})}\nonumber\\
&\propto \prod_{i=1}^m h( {{\boldsymbol w}}_{(i)}, {{\boldsymbol \theta}}),\end{aligned}$$ where $h$ is a parametric angular density function [e.g. @engelke2012; @beirlant2006 pp. 170–171]. In the next section the angular density models described in Section \[s:par\_models\] are fitted to the data by the maximization of the likelihood. For brevity the asymmetric logistic model is not considered since it has too many parameters. The likelihood is proportional to the product of angular densities, therefore the maximizer of is obtained equivalently by maximizing the log-likelihood $$\label{eq:aploglik}
\ell({{\boldsymbol \theta}})=\sum_{i=1}^{m}\log h({{\boldsymbol w}}_{(i)},{{\boldsymbol \theta}}).$$ Denote by $\widehat{{{\boldsymbol \theta}}}$ the maximizer of $\ell$ and by $\ell'({{\boldsymbol \theta}}) = \nabla_{{{\boldsymbol \theta}}}\,\ell({{\boldsymbol \theta}})$ the score function. Since provides an approximation of the true likelihood, then from the theory on model misspecification [e.g. @davison03 pp. 147–148] it follows that $$\label{eq:godambe}
\sqrt{n}(\widehat{{{\boldsymbol \theta}}}-{{\boldsymbol \theta}})\stackrel{d}{\rightarrow} \mathcal{N}_p({{\boldsymbol 0}},J({{\boldsymbol \theta}})^{-1}\,K({{\boldsymbol \theta}})\,J({{\boldsymbol \theta}})^{-1}),
\quad n\rightarrow\infty,$$ where $\mathcal{N}_p({{\boldsymbol \mu}},\Sigma)$ is the $p$-dimensional normal distribution with mean ${{\boldsymbol \mu}}$ and covariance $\Sigma$, ${{\boldsymbol \theta}}$ is the true parameter and $$J({{\boldsymbol \theta}})=-\E\lbrace\nabla_{{{\boldsymbol \theta}}} \ell'({{\boldsymbol \theta}})\rbrace,\qquad
K({{\boldsymbol \theta}})=\operatorname{Var}_{{{\boldsymbol \theta}}} \lbrace \ell' ( {{\boldsymbol \theta}}) \
m) 
 & = \frac{e^{-\psi (\cW_{r_0 }) } \psi(\cW_{r_0})^m}{m! } 
 \prod_{i=1}^m \frac{{{\mathrm d}}\psi(r_{(i) }, { { \boldsymbol w}}_{(i)})}{\psi(\cW_{r_0})}\nonumber\\ 
 & \propto \prod_{i=1}^m h ({ { \boldsymbol w}}_{(i) }, { { \boldsymbol \theta}}),\end{aligned}$$ where $ h$ is a parametric angular density function [ e.g. @engelke2012; @beirlant2006 pp. 170–171 ]. In the next incision the angular concentration models described in incision \[s: par\_models\ ] are fit to the data by the maximization of the likelihood. For brevity the asymmetric logistic mannequin is not think since it have too many parameters. The likelihood is proportional to the intersection of angular densities, therefore the maximizer of is obtained equivalently by maximize the log - likelihood $ $ \label{eq: aploglik } 
 \ell({{\boldsymbol \theta}})=\sum_{i=1}^{m}\log h({{\boldsymbol w}}_{(i)},{{\boldsymbol \theta}}).$$ Denote by $ \widehat{{{\boldsymbol \theta}}}$ the maximizer of $ \ell$ and by $ \ell'({{\boldsymbol \theta } }) = \nabla_{{{\boldsymbol \theta}}}\,\ell({{\boldsymbol \theta}})$ the score function. Since provides an approximation of the truthful likelihood, then from the theory on model misspecification [ e.g. @davison03 pp. 147–148 ] it be that $ $ \label{eq: godambe } 
 \sqrt{n}(\widehat{{{\boldsymbol \theta}}}-{{\boldsymbol \theta}})\stackrel{d}{\rightarrow } \mathcal{N}_p({{\boldsymbol 0}},J({{\boldsymbol \theta}})^{-1}\,K({{\boldsymbol \theta}})\,J({{\boldsymbol \theta}})^{-1 }), 
 \quad n\rightarrow\infty,$$ where $ \mathcal{N}_p({{\boldsymbol \mu}},\Sigma)$ is the $ p$-dimensional normal distribution with mean $ { { \boldsymbol \mu}}$ and covariance $ \Sigma$, $ { { \boldsymbol \theta}}$ is the true parameter and $ $ J({{\boldsymbol \theta}})=-\E\lbrace\nabla_{{{\boldsymbol \theta } } } \ell'({{\boldsymbol \theta}})\rbrace,\qquad 
 K({{\boldsymbol \theta}})=\operatorname{Var}_{{{\boldsymbol \theta } } } \lbrace \ell' ({ { \boldsymbol \theta } }) \
m) 
&= \vrac{e^{-\psi (\cW_{r_0})} \psi(\cW_{r_0})^m}{m!}
\proa_{i=1}^m \frac{{{\mathrm b}}\psi(r_{(i)}, {{\uoldsymgol w}}_{(i)})}{\psk(\cW_{r_0})}\nonumber\\
&\propto \prod_{i=1}^m h( {{\uoldwymboo w}}_{(i)}, {{\boldsymbol \theta}}),\evd{aligned}$$ where $h$ is e parametric angnmar density fhkctiou [x.g. @engelke2012; @beirkant2006 pp. 170–171]. It the next secdiun the angular density models describqd in Srchion \[s:par\_modelf\] art fyttes to the data by the maximization kf the nikelihood. Fot brevity the asymmetric llgishic model is not clnsidered sunce ut has too mxny parameuexs. The likemihood is proportional to the pfoducc of angulat sejvities, therxfore nhe maximizer of is ottained equivalently ny maeimizing the log-likelihoov $$\label{eq:aploglik}
\ell({{\foldsymbon \cheta}})=\sum_{i=1}^{m}\log h({{\boldsynbil w}}_{(i)},{{\toldvymbuo \tfetz}}).$$ Venkte by $\wivehat{{{\boldsyjbol \theta}}}$ rhe maximizer of $\elk$ wbd by $\ell'({{\boldaymbol \treta}}) = \nabla_{{{\boldsymbol \theta}}}\,\ell({{\boldsymbon \tgeta}})$ the score function. Since provides an apkroximatiog of the true likelihood, then from the theory on kodel oiswpccifkxahion [e.g. @davison03 pp. 147–148] it follows that $$\label{eq:goqzmne}
\xqrt{n}(\widehat{{{\bojdsymbol \thrtw}}}-{{\bpjdsymbol \thetx}})\stackxsl{s}{\rightarrow} \mathcap{N}_p({{\boldfymboo 0}},J({{\boldsyibol \theta}})^{-1}\,K({{\boldsymbol \theta}})\,J({{\bolesymbol \thetc}})^{-1}),
\quqd n\rightarrow\inftv,$$ where $\matheal{N}_p({{\bpldsykbol \mu}},\Sigma)$ is the $p$-dioensjonal normap distribhgion with mean ${{\buldxykbol \mu}}$ and covariance $\Sigia$, ${{\boldsynbol \theta}}$ ks tne truq parameteg and $$J({{\boldsymbol \theta}})=-\E\pbracg\nabla_{{{\toldsymbol \theta}}} \ell'({{\boldsymbol \theta}})\rbracx,\aquad
K({{\boldsykbml \nheta}})=\operctornake{Var}_{{{\boldsymbjl \theta}}} \lbracg \ell' ( {{\bojdsymcol \theta}}) \
m) &= \frac{e^{-\psi (\cW_{r_0})} \psi(\cW_{r_0})^m}{m!} \prod_{i=1}^m \frac{{{\mathrm w}}_{(i)})}{\psi(\cW_{r_0})}\nonumber\\ \prod_{i=1}^m h( w}}_{(i)}, {{\boldsymbol \theta}}),\end{aligned}$$ angular function [e.g. @engelke2012; pp. 170–171]. In next section the angular density models in Section \[s:par\_models\] are fitted to the data by the maximization of the For brevity the asymmetric logistic model is not considered since it has too parameters. likelihood proportional the product of angular densities, therefore the maximizer of is obtained equivalently by maximizing the log-likelihood \ell({{\boldsymbol \theta}})=\sum_{i=1}^{m}\log h({{\boldsymbol w}}_{(i)},{{\boldsymbol \theta}}).$$ Denote by $\widehat{{{\boldsymbol the maximizer of $\ell$ by $\ell'({{\boldsymbol \theta}}) = \nabla_{{{\boldsymbol \theta}})$ score function. provides approximation the true likelihood, from the theory on model misspecification [e.g. @davison03 pp. 147–148] it follows that $$\label{eq:godambe} \sqrt{n}(\widehat{{{\boldsymbol \theta}}}-{{\boldsymbol \theta}})\stackrel{d}{\rightarrow} 0}},J({{\boldsymbol \theta}})^{-1}\,K({{\boldsymbol \quad n\rightarrow\infty,$$ $\mathcal{N}_p({{\boldsymbol is $p$-dimensional normal distribution ${{\boldsymbol \mu}}$ and covariance $\Sigma$, ${{\boldsymbol true parameter and $$J({{\boldsymbol \theta}})=-\E\lbrace\nabla_{{{\boldsymbol \theta}}} \ell'({{\boldsymbol \theta}})\rbrace,\qquad \theta}})=\operatorname{Var}_{{{\boldsymbol \theta}}} \ell' ( {{\boldsymbol \theta}}) \
m) 
&= \frac{e^{-\psi (\cW_{r_0})} \psi(\cW_{r_0})^m}{m!}
\prod_{i=1}^M \frac{{{\mathrM d}}\psi(R_{(i)}, {{\bOldSyMbol W}}_{(i)})}{\psI(\cW_{r_0})}\nonumber\\
&\prOPto \pRod_{i=1}^m h( {{\boldsymbol w}}_{(i)}, {{\boldSymboL \tHEta}}),\eND{aLigneD}$$ where $h$ IS a PARamEtRiC anGuLAr DensiTy fUnction [E.g. @engelke2012; @bEirLaNt2006 pp. 170–171]. In the nexT SeCtion the anGulAr density modEls DescriBeD in sEctioN \[s:pAr\_modEls\] are FItted tO the data bY tHE maximIZation oF THe LikeLihood. For brevity tHE aSYmmetric logistIc modeL iS NoT COnsIdeRed since it HaS too mANy paramETeRS. tHe lIKelihood is proPortional to THe pRoduct Of AngULar denSitieS, tHEreFore the maxiMizeR of is obtaIned eqUIvalentLY by maxiMizing The Log-LikeLIhOoD $$\laBeL{Eq:aPLoGliK}
\Ell({{\BoldsymbOl \ThEta}})=\suM_{i=1}^{m}\lOG H({{\BOldsYmbOl w}}_{(i)},{{\BoldsYmbol \theta}}).$$ DenOte By $\wiDEhaT{{{\boldSymboL \theTa}}}$ The maXimizeR of $\elL$ aNd by $\ell'({{\boldsymbOl \thEta}}) = \nabla_{{{\bOldSyMboL \tHeta}}}\,\eLL({{\boldsYmbOl \tHeta}})$ the Score fuNCtiOn. sINCe Provides an approximAtION oF the true LikeliHOoD, tHEn from thE tHeoRy on MODel miSspeCIfIcation [e.G. @davisON03 pP. 147–148] iT followS tHat $$\labEl{Eq:gOdaMbe}
\sqRT{n}(\wiDehat{{{\bOldsymboL \thetA}}}-{{\Boldsymbol \thetA}})\Stackrel{d}{\righTArROW} \mAThcaL{N}_p({{\Boldsymbol 0}},J({{\BoldSYmboL \theTA}})^{-1}\,K({{\BolDSymboL \thetA}})\,J({{\BOlDSymbol \theta}})^{-1}),
\quad n\rigHtArrow\iNfty,$$ wHere $\mathcal{N}_p({{\Boldsymbol \MU}},\sIgma)$ is thE $p$-diMEnSIonal normal disTribuTion with meAN ${{\boldsymBol \mu}}$ And covarIance $\SigmA$, ${{\BOldsymboL \thEta}}$ Is tHe tRUE pArameter and $$J({{\bOLDsymBoL \theta}})=-\E\LbrAce\nablA_{{{\boLdsYmbOl \tHeTa}}} \ell'({{\boldSymbol \thEtA}})\rBrAcE,\qqUad
K({{\bOLdsymbol \ThEta}})=\OpEraTornaME{Var}_{{{\boLdsymBol \tHeTa}}} \LBraCe \ell' ( {{\boLDsYMBol \tHeTa}}) \
m) 
&= \frac{e^{-\psi (\cW _{r_0})} \ psi(\ cW_ {r_ 0} )^m} {m!} 
\prod_{i=1}^m \fra c{{{\mathrm d}}\psi(r_ {(i)} ,{ {\bo l ds ymbol  w}}_{( i )} ) } {\p si (\ cW_ {r _ 0} )}\no num ber\\
& \propto \p rod _{ i=1}^m h( {{ \ bo ldsymbol w }}_ {(i)}, {{\bo lds ymbol\t het a }}),\ end {alig ned}$$ where$h$ is apa r ametri c  angula r de nsit y function [e.g.@ en g elke2012; @bei rlant2 00 6  p p .  17 0–1 71]. In th enexts ectiont he a n gul a r density mod els describ e d i n Sect io n \ [ s:par\ _mode ls \ ] a re fitted t o th e data by  the m a ximizat i on of t he lik eli hoo d. F o rbr evi ty the as ymm e tri c logist ic  m odelis n o t c onsi der ed s inceit has too ma nypara m ete rs. T he li keli ho od is  propo rtion al  to the product  ofangular d ens it ies ,there f ore th e m axi mizer o f is ob t ain ed e q ui valently by maximi zi n g  t he log-l ikelih o od  $ $ \label{e q: apl ogli k } 
\ell ({{\ b ol dsymbol\theta } }) =\ sum_{i= 1} ^{m}\l og  h( {{\ bolds y mbol  w}}_{ (i)},{{\ bolds y mbol \theta}}) . $$ Denote by$ \w i d eh a t{{{ \bo ldsymbol \t heta } }}$them ax imi z er of  $\el l$ an d  by $\ell'({{\bolds ym bol \t heta} }) = \nabla_{ {{\boldsym b o l  \theta} }}\, \ el l ({{\boldsymbol  \the ta}})$ the score fu nctio n. Since  provides a n approx ima tio n o f t h e  t rue likelihoo d ,  the nfrom th e t heory o n m ode l m iss pe cificatio n [e.g.@d av is on 03pp. 1 4 7–148] i tfol lo wsthat$ $\labe l{eq: goda mb e} \sq rt{n}(\ w id e h at{{ {\ bo ldsy mbo l\thet a}}} - {{\ boldsym bol \thet a}} ) \sta ck re l{d}{\r ightarrow} \m at hcal{N}_p( {{ \bo ldsymb o l  0}},J({ {\boldsymbol \theta}})^ { -1}\,K( {{\ bolds ymbo l \theta} })\ ,J({{\ bol d symbol  \thet a}})^ {- 1}) , \quad n \r igh ta rrow\infty , $ $ w here$\ math cal{N}_ p({{\boldsymbol \m u }}, \Sigma)$ is t he$p$- d i me nsi o na l  no rm a l d i s tribution withmean ${{\b ol d sy mbol \mu}} $  an dcovaria nce $\S igma$ ,  ${{\bo ldsymbol\theta}}$  i s th e tru e paramete r and $$ J({{\bold s ymbol \t heta} })= -\E\lb ra ce\ nabla _{{{\b o lds ymbol  \thet a} }} \el l'({{ \b oldsymbo l \theta}})\rbrace,\qqu ad
K({ {\bol dsy mbol \the ta} } )=\ operatorn ame{ Var}_{{{\b old sym bol \ the t a}}}\lbr a ce  \e l l' ({{\b o ldsymbol\ th eta } } )\
m) 
&=_\frac{e^{-\psi (\cW_{r_0})}_\psi(\cW_{r_0})^m}{m!}
\prod_{i=1}^m \frac{{{\mathrm d}}\psi(r_{(i)}, {{\boldsymbol_w}}_{(i)})}{\psi(\cW_{r_0})}\nonumber\\
&\propto \prod_{i=1}^m_h(_{{\boldsymbol w}}_{(i)},_{{\boldsymbol_\theta}}),\end{aligned}$$ where $h$_is a parametric_angular density function [e.g._@engelke2012; @beirlant2006 pp._170–171]._In the next section the angular density models described in Section \[s:par\_models\] are fitted_to_the data_by_the_maximization of the likelihood. For_brevity the asymmetric logistic model_is not_considered since it has too many parameters. The_likelihood_is proportional to_the product of angular densities, therefore the maximizer of_is obtained equivalently by maximizing the_log-likelihood $$\label{eq:aploglik}
\ell({{\boldsymbol \theta}})=\sum_{i=1}^{m}\log_h({{\boldsymbol_w}}_{(i)},{{\boldsymbol_\theta}}).$$ Denote by $\widehat{{{\boldsymbol_\theta}}}$ the maximizer of $\ell$ and_by $\ell'({{\boldsymbol \theta}}) = \nabla_{{{\boldsymbol \theta}}}\,\ell({{\boldsymbol_\theta}})$ the score function. Since provides an_approximation of the true likelihood, then_from the theory on model_misspecification [e.g._@davison03 pp. 147–148] it follows_that $$\label{eq:godambe}
\sqrt{n}(\widehat{{{\boldsymbol \theta}}}-{{\boldsymbol_\theta}})\stackrel{d}{\rightarrow} \mathcal{N}_p({{\boldsymbol_0}},J({{\boldsymbol \theta}})^{-1}\,K({{\boldsymbol \theta}})\,J({{\boldsymbol_\theta}})^{-1}),
\quad n\rightarrow\infty,$$ where $\mathcal{N}_p({{\boldsymbol \mu}},\Sigma)$ is_the $p$-dimensional normal_distribution with mean ${{\boldsymbol \mu}}$ and_covariance_$\Sigma$, ${{\boldsymbol \theta}}$_is_the_true parameter_and $$J({{\boldsymbol \theta}})=-\E\lbrace\nabla_{{{\boldsymbol_\theta}}}_\ell'({{\boldsymbol \theta}})\rbrace,\qquad
K({{\boldsymbol_\theta}})=\operatorname{Var}_{{{\boldsymbol_\theta}}} \lbrace \ell' ( {{\boldsymbol \theta}})_\