 function by the inclusion-exclusion principle, then using we have $$\label{eq:est_prob_ruin_set}
\P\{Y_1>y_1,\ldots,Y_d>y_d\}\approx  \sum_{j=1}^d  \frac{1}{y_j} \bar{\Phi}_{d-1} 
\Big\lbrace \Big (\lambda_{k,j} +\frac{\log y_k/y_j }{2\lambda_{k,j}} \Big )_{k \in I_j}; \bar{\Lambda}_j \Big\rbrace,$$ where $\bar{\Phi}_{d-1}$ is the survival function of the multivariate normal distribution [@nikoloulopoulos2009]. Similar to @cooley2010 we define three extreme events: $\{\text{PM}10 > 95, \text{NO} > 270, \text{SO}2 > 95\}$, $\{\text{NO}2 > 110, \text{SO}2 > 95, \text{NO} > 270\}$ and $\{\text{PM}10 > 95, \text{NO} > 270, \text{NO}2 > 110, \text{SO}2 > 95\}$. Then, we compute probability using in place of the parameters their estimates. Table \[table:proba\_excess\] reports the results. For the three events the estimates fall inside the $95\%$ confidence intervals highlighting the ability of the model to estimate such extreme events.

                 Event $1$                    Event $2$                    Event $3$
  -------------- ---------------------------- ---------------------------- -----------------------------
  Excess / $n$   $18 /528 $                   $14 / 562$                   $12 /528$
  Emp. Est.      $0.034\; ( 0.019, 0.050 )$   $0.025\; ( 0.012, 0.038 )$   $0.023\; ( 0.010, 0.035 )$
  Mod. Est.      $0.038$                      $0.030$                      $0.030$

  : Probability estimates of excesses. The first row reports the number of excess and the sample size. The second row reports the empirical estimates and between brackets the $95\%$ confidence intervals obtained with the normal approximation. The third row reports the model estimates.[]{data-label="table:proba_excess"}

$
\begin{array}{
function by the inclusion - exclusion principle, then use we take $ $ \label{eq: est_prob_ruin_set } 
 \P\{Y_1 > y_1,\ldots, Y_d > y_d\}\approx   \sum_{j=1}^d   \frac{1}{y_j } \bar{\Phi}_{d-1 } 
 \Big\lbrace \Big (\lambda_{k, j } + \frac{\log y_k / y_j } { 2\lambda_{k, j } } \Big) _ { k \in I_j }; \bar{\Lambda}_j \Big\rbrace,$$ where $ \bar{\Phi}_{d-1}$ is the survival function of the multivariate normal distribution [ @nikoloulopoulos2009 ]. Similar to @cooley2010 we specify three extreme events: $ \{\text{PM}10 > 95, \text{NO } > 270, \text{SO}2 > 95\}$, $ \{\text{NO}2 > 110, \text{SO}2 > 95, \text{NO } > 270\}$ and $ \{\text{PM}10 > 95, \text{NO } > 270, \text{NO}2 > 110, \text{SO}2 > 95\}$. Then, we compute probability use in position of the parameters their estimate. Table \[table: proba\_excess\ ] reports the resultant role. For the three events the estimates decrease inside the $ 95\%$ confidence intervals highlighting the ability of the model to estimate such extreme events. 

                  Event $ 1 $                     Event $ 2 $                     Event $ 3 $ 
   -------------- ---------------------------- ---------------------------- ----------------------------- 
   Excess / $ n$    $ 18 /528 $                    $ 14 / 562 $                    $ 12 /528 $ 
   Emp. Est.       $ 0.034\; (0.019, 0.050) $    $ 0.025\; (0.012, 0.038) $    $ 0.023\; (0.010, 0.035) $ 
   Mod. Est.       $ 0.038 $                       $ 0.030 $                       $ 0.030 $ 

  : Probability estimate of excesses. The first row report the number of excess and the sample size. The second row report the empirical estimates and between brackets the $ 95\%$ confidence intervals prevail with the normal approximation. The third row reports the model estimates.[]{data - label="table: proba_excess " } 

 $ 
 \begin{array } {
 fujction by the inclusion-ewclusion principle, then using we have $$\label{eq:est_prob_ruin_set}
\P\{Y_1>y_1,\ldovs,Y_d>t_d\}\appeox  \sum_{j=1}^d  \frac{1}{y_j} \bar{\Ohi}_{d-1} 
\Big\lhrace \Bit (\lanvda_{k,j} +\frac{\log y_k/y_j }{2\lambdz_{n,j}} \Bng )_{k \in I_j}; \bar{\Lakbda}_j \Big\rtrace,$$ where $\bas{\Pfi}_{b-1}$ is the survival function of the mujtivaristf normal distrybutpog [@niioloulopoulos2009]. Similar to @cooley2010 we define three extremr events: $\{\text{PM}10 > 95, \text{NO} > 270, \texh{SO}2 > 95\}$, $\{\text{NO}2 > 110, \texh{SO}2 > 95, \text{NI} > 270\}$ wbd $\{\text{PM}10 > 95, \gext{NO} > 270, \ttxc{NO}2 > 110, \text{SK}2 > 95\}$. Then, we compute probability usiny in place if thf parameters theig estimates. Table \[tabne:proba\_rxcess\] reports thx rewults. For the three etents the estimates sall insige the $95\%$ confidence unrervans hhghlkthtkng tie zbilitj oh the model to estimatw such extreme evenus.

                 Evene $1$                    Event $2$                    Event $3$
  -------------- ---------------------------- ---------------------------- -----------------------------
  Excess / $n$   $18 /528 $                   $14 / 562$                   $12 /528$
  Emp. Est.      $0.034\; ( 0.019, 0.050 )$   $0.025\; ( 0.012, 0.038 )$   $0.023\; ( 0.010, 0.035 )$
  Jud. Tsb.      $0.038$                      $0.030$                      $0.030$

  : Pgobability estimabes of excesses. Thr vitft row reportr the uhmger of excess and hhe samkle size. The sesond row reports the empirical wstimates anb bwtween brackets thz $95\%$ confidencz intetvals pbtained with the normau aplroximation. The thirs row reports the mocen estimauds.[]{data-label="table:pwoba_excesw"}

$
\begnn{array}{
function by the inclusion-exclusion principle, then using $$\label{eq:est_prob_ruin_set} \sum_{j=1}^d \frac{1}{y_j} \Big\lbrace \Big (\lambda_{k,j} \in \bar{\Lambda}_j \Big\rbrace,$$ where is the survival of the multivariate normal distribution [@nikoloulopoulos2009]. to @cooley2010 we define three extreme events: $\{\text{PM}10 > 95, \text{NO} > 270, > 95\}$, $\{\text{NO}2 > 110, \text{SO}2 > 95, \text{NO} > 270\}$ and $\{\text{PM}10 95, > \text{NO}2 110, \text{SO}2 > 95\}$. Then, we compute probability using in place of the parameters their estimates. \[table:proba\_excess\] reports the results. For the three events estimates fall inside the confidence intervals highlighting the ability the to estimate extreme Event Event $2$ Event -------------- ---------------------------- ---------------------------- ----------------------------- Excess / $n$ $18 /528 $ $14 / 562$ $12 /528$ Emp. Est. ( 0.019, $0.025\; ( 0.038 $0.023\; 0.010, 0.035 )$ $0.038$ $0.030$ $0.030$ : Probability estimates first row reports the number of excess and sample size. second row reports the empirical estimates between brackets the $95\%$ confidence intervals obtained with normal approximation. The third row reports the model estimates.[]{data-label="table:proba_excess"} $ \begin{array}{
 function by the inclusion-excLusion prinCiple, TheN usInG we hAve $$\lAbel{eq:est_prob_rUIn_seT}
\P\{Y_1>y_1,\ldots,Y_d>y_d\}\approx  \sum_{J=1}^d  \fraC{1}{y_J} \Bar{\PHI}_{d-1} 
\big\lbRace \Big (\LAmBDA_{k,j} +\FrAc{\Log Y_k/Y_J }{2\lAmbda_{K,j}} \BIg )_{k \in I_j}; \Bar{\Lambda}_j \big\RbRace,$$ where $\bar{\pHi}_{D-1}$ is the survIvaL function of tHe mUltivaRiAte NOrmal DisTribuTion [@niKOlouloPoulos2009]. SimIlAR to @cooLEy2010 we defINE tHree Extreme events: $\{\text{pm}10 > 95, \tEXt{NO} > 270, \text{SO}2 > 95\}$, $\{\text{nO}2 > 110, \text{sO}2 > 95, \TExT{no} > 270\}$ anD $\{\teXt{PM}10 > 95, \text{NO} > 270, \TeXt{NO}2 > 110, \tEXt{SO}2 > 95\}$. TheN, We COMPutE Probability usIng in place oF The ParameTeRs tHEir estImateS. TABle \[Table:proba\_eXcesS\] reports tHe resuLTs. For thE Three evEnts thE esTimAtes FAlL iNsiDe THe $95\%$ cONfIdeNCe iNtervals HiGhLightIng tHE ABIlitY of The mOdel tO estimate such ExtReme EVenTs.

                 EveNt $1$                    EveNt $2$                    EvEnT $3$
  -------------- ---------------------------- ---------------------------- -----------------------------
  ExceSs / $n$   $18 /528 $                   $14 / 562$                   $12 /528$
  Emp. est.      $0.034\; ( 0.019, 0.050 )$   $0.025\; ( 0.012, 0.038 )$   $0.023\; ( 0.010, 0.035 )$
  MoD. ESt.      $0.038$                      $0.030$                      $0.030$

  : Probability esTimaTes of exceSseS. THe fIrSt row REports The NumBer of exCess and THe sAmPLE SiZe. The second row repoRtS THe EmpiricaL estimATeS aND between BrAckEts tHE $95\%$ ConfiDencE InTervals oBtaineD WiTh The normAl ApproxImAtiOn. THe thiRD row ReportS the modeL estiMAtes.[]{data-label="tABle:proba_excesS"}

$
\BeGIN{aRRay}{
 function by the inclusion -exclusion  prin cip le,  t henusin g we have $$\l a bel{ eq:est_prob_ruin_set}\P\{Y _1 > y_1, \ ld ots,Y _d>y_d\ } \a p p rox  \s um_ {j = 1} ^d  \ fra c{1}{y_ j} \bar{\P hi} _{ d-1} 
\Big\l b ra ce \Big (\ lam bda_{k,j} +\ fra c{\log  y _k/ y _j }{ 2\l ambda _{k,j} }  \Big)_{k \inI_ j }; \ba r {\Lambd a } _j  \Bi g\rbrace,$$ where $\ b ar{\Phi}_{d-1} $ is t he su r v iva l f unction of  t he mu l tivaria t en o r mal distribution[@nikoloulo p oul os2009 ].  Si m ilar t o @co ol e y20 10 we defin e th ree extre me eve n ts: $\{ \ text{PM }10 >95,  \t ext{ N O}  >  27 0, \te x t{ SO} 2  >95\}$, $ \{ \t ext{N O}2> 1 1 0, \ tex t{SO }2 >95, \text{NO}  >270\ } $ a nd $\ {\tex t{PM }1 0 > 9 5, \te xt{NO }> 270, \text{NO }2 >  110, \te xt{ SO }2>95\}$ .  Then,  we  co mpute p robabil i tyus i n g  i n place of the par am e t er s theirestima t es .T able \[t ab le: prob a \ _exce ss\] re ports th e resu l ts .For the  t hree e ve nts  th e est i mate s fall  insidethe $ 9 5\%$ confidenc e  intervals hi g hl i g ht i ng t heability ofthem odel  toe st ima t e suc h ext re m ee vents.

               Ev ent $ 1$                   E v e n t $2$               Event  $3$  -------- - ----- -- ----- -------- --------- - - -- ----- --- --- --- --- - - -- ------- ----- - - ---- -- ------- --- ------  E xce ss/ $ n$    $18 /5 28 $             $ 1 4 / 562$              $12  /528 $
 Em p. Est .     $ 0. 0 3 4\;(0. 019,  0. 05 0 )$  $0 . 025 \; ( 0. 012, 0.03 8 ) $    $ 0. 02 3\; ( 0 .010, 0.035 ) $  Mod. Est .    $0.0 3 8 $                     $0.030$                $0.030$
 : Prob abi l ity es timate s ofex ces s e s. Th e fi rst  r ow reports t henumbe rof e xcess a nd the sample size .  Th e second rowrep orts t he  em p ir i cal  e s tim a t es and betweenbrackets t he $9 5\%$ confi d enc einterva ls obta inedw ith the  normal a pproximat io n. T h e  th ird row re ports th e model e s timat e s. []{da ta- label= "t abl e:pro ba_exc e ss" }

$\begin {a rray}{ 
 function_by the_inclusion-exclusion principle, then using_we have_$$\label{eq:est_prob_ruin_set}
\P\{Y_1>y_1,\ldots,Y_d>y_d\}\approx_ \sum_{j=1}^d__\frac{1}{y_j} \bar{\Phi}_{d-1} 
\Big\lbrace_\Big (\lambda_{k,j} +\frac{\log_y_k/y_j }{2\lambda_{k,j}} \Big )_{k_\in I_j}; \bar{\Lambda}_j_\Big\rbrace,$$_where $\bar{\Phi}_{d-1}$ is the survival function of the multivariate normal distribution [@nikoloulopoulos2009]. Similar to_@cooley2010_we define_three_extreme_events: $\{\text{PM}10 > 95, \text{NO}_> 270, \text{SO}2 > 95\}$,_$\{\text{NO}2 >_110, \text{SO}2 > 95, \text{NO} > 270\}$ and_$\{\text{PM}10_> 95, \text{NO}_> 270, \text{NO}2 > 110, \text{SO}2 > 95\}$. Then,_we compute probability using in place_of the parameters_their_estimates._Table \[table:proba\_excess\] reports the_results. For the three events the_estimates fall inside the $95\%$ confidence_intervals highlighting the ability of the model_to estimate such extreme events.

 _    _ _    _  _Event $1$_  _     _  _     __Event $2$ ___ _  __ __     __Event $3$
  -------------- ---------------------------- ---------------------------- -----------------------------
_ Excess / $n$__ $18 /528 $_    _     _ _ _   $14 / 562$     _     _   _ __$12 /528$
 _Emp._Est. _    $0.034\; ( 0.019,_0.050 )$_  $0.025\; ( 0.012, 0.038_)$   $0.023\;_(_0.010, 0.035 )$
  Mod. Est._     $0.038$_    ___    _     _  _$0.030$      __     __     __$0.030$

_ : Probability estimates of_excesses. The first row reports the_number of excess and the sample size. The second_row reports the_empirical estimates and between brackets_the_$95\%$_confidence intervals obtained with the normal approximation. The third row_reports the_model estimates.[]{data-label="table:proba_excess"}

$
\begin{array}{