Page 291 - Invited Paper Session (IPS) - Volume 2
P. 291

IPS246 Tiziana Laureti et al.
               design,  the  weighting  procedure,  the  imputation  schemes  and  the
               nonlinear form of survey estimators should be reflected in the calculation
               of standard errors and confidence intervals (Goedemé, 2013).
              a.  Estimating standard errors: An Empirical analysis for Italian regions
                     In Italy, EU-SILC is conducted by ISTAT to produce estimates of the
                  Italian population living conditions at national and also at regional level
                  (NUTS-2).  In  the  design  of  the  EU-SILC  survey,  regions  are  planned
                  domains  for  which  estimates  are  published  but  without  confidence
                  intervals (ISTAT, 2018). The regional samples are based on a stratified
                  two-stage  sample  design:  in  each  province,  municipalities  are  the
                  Primary  Sampling  Units  (PSUs),  while  households  are  the  Secondary
                  Sampling  Units  (SSUs).  The  PSUs  are  stratified  according  to
                  administrative regions and population size while SSUs are selected by
                  using  systematic  sampling.  We  calculate  confidence  intervals  for  the
                  AROP  for  Italian  regions  in  2017  (NUTS-2)  using  data  from  IT-SILC
                  provided by ISTAT for carrying out research projects in collaboration
                  with  the  Dagum  /Tuscan  Universities  Research  Centre  on  Advanced
                  Statistics for the Equitable and Sustainable Development.
                     AROP is a complex statistic since it is based on a poverty threshold
                  computed from the median of the income distribution, that is:  =
                  ( < 0.6 ∙  0.5 ) ∙ 100. Therefore, the at-risk-of-poverty threshold (ARPT)
                  needs to be estimated first, which is set at 60% of the national median
                                                 ̂
                  equivalized disposable income,  = 0.6 ∙  0.5 . Then the AROP rate is
                  defined  as  the  proportion  of  persons  with  an  equivalized  disposable
                                           ̂
                                                       ̂
                  income below the ARPT:  =  ∑ ∈<    ∙ 100
                                                      ∑  
                     Consequently, there exist two main sources of variability: one is due
                  to the estimated threshold and the other one comes from the estimated
                  proportion given the estimated threshold (e.g., Berger and Skinner 2003;
                  Verma et al 2012). We used a generalized linearization method based
                  on the concept of influence Function (Deville, 1999) which allows us to
                  deal with nonlinear statistics for which the Taylor method cannot be
                  used. This approach does not involve more calculations.
                     We linked household and personal data (H, R, P D files). Household
                  income is equivalised using the modified OECD equivalence scale by
                  assuming that the living standard of all household members is the same.
                  Confidence intervals are estimated using linearization on the basis of
                  complete sample design information regarding “Primary strata”, “PSUs”,
                  “SSU”. The estimates are obtained using R package laeken and the DASP
                  module. Since a region coincides with a sample “design domain” and the
                  regional  AROP  depends  only  on  units  within  region,  variance




                                                               278 | I S I   W S C   2 0 1 9
   286   287   288   289   290   291   292   293   294   295   296