* Encoding: UTF-8.
recode area_sanitaria (1=1) (2=2) (3=3) (4=4) (5=5) (6=6) (7=8) (8=7) into area_sanitaria_toledo.
value lables area_sanitaria_toledo 1 'Albacete' 2 'La Mancha-Centro' 3 'Ciudad Real' 4' Cuenca' 5 'Guadalajara' 6 'Talavera de la Reina' 7 'Puertollano' 8'Toledo'.


*compute ln_plazas=ln(plazas_totales).

desc plazas_totales ln_plazas.
* Tablas 1 y 2**********************************************************************************************************************.
USE ALL.
COMPUTE filter_$=(excludecases = 0).
VARIABLE LABELS filter_$ 'excludecases = 0 (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.

DATASET ACTIVATE ConjuntoDatos1.
FREQUENCIES VARIABLES=area_sanitaria
  /ORDER=ANALYSIS.

MEANS TABLES=PLAZAS_TOTALES BY area_sanitaria
  /CELLS=MEAN COUNT STDDEV.

CROSSTABS
  /TABLES=area_sanitaria BY titgest
  /FORMAT=AVALUE TABLES
  /CELLS=COUNT
  /COUNT ROUND CELL.

******Tabla 2*************************************************************************************************************************************.
 

means defunciones BY area_sanitaria
 /CELLS=MEAN COUNT STDDEV SUM
  /STATISTICS ANOVA.
means mortalidad BY area_sanitaria
 /CELLS=MEAN COUNT STDDEV
  /STATISTICS ANOVA.

DATASET ACTIVATE ConjuntoDatos1.
*Nonparametric Tests: Independent Samples. 
MEANS TABLES=inc_covid1865_1 inc_covid_1865_2 BY area_sanitaria
  /CELLS=MEAN COUNT MEDIAN RANGE.


NPTESTS 
  /INDEPENDENT TEST (inc_covid1865_1 inc_covid_1865_2) GROUP (area_sanitaria) 
    KRUSKAL_WALLIS(COMPARE=PAIRWISE) MEDIAN(TESTVALUE=SAMPLE COMPARE=PAIRWISE) 
  /MISSING SCOPE=ANALYSIS USERMISSING=EXCLUDE
  /CRITERIA ALPHA=0.05  CILEVEL=95.

USE ALL.
COMPUTE filter_$=(mvCOVIDdeaths = 0).
VARIABLE LABELS filter_$ 'mvCOVIDdeaths = 0 (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.


means defunciones_COVID  BY area_sanitaria
 /CELLS=MEAN COUNT STDDEV SUM
  /STATISTICS ANOVA.
means mortalidad_COVID  BY area_sanitaria
 /CELLS=MEAN COUNT STDDEV
  /STATISTICS ANOVA.

******Tabla 3 *****************************************************************************************************.************.

USE ALL.
COMPUTE filter_$=(excludecases = 0).
VARIABLE LABELS filter_$ 'excludecases = 0 (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.

MEANS TABLES=PLAZAS_TOTALES BY titgest
  /CELLS=MEAN COUNT STDDEV MEDIAN SUM.

* Modelos lineales generalizados.
GENLIN defunciones BY area_sanitaria_toledo (ORDER=ASCENDING) 
  /MODEL area_sanitaria_toledo  INTERCEPT=YES 
    SCALEWEIGHT=ln_plazas
 DISTRIBUTION=POISSON LINK=LOG
  /CRITERIA METHOD=FISHER(1) SCALE=1 COVB=MODEL MAXITERATIONS=100 MAXSTEPHALVING=5 
    PCONVERGE=1E-006(ABSOLUTE) SINGULAR=1E-012 ANALYSISTYPE=3(WALD) CILEVEL=95 CITYPE=WALD 
    LIKELIHOOD=FULL
  /EMMEANS TABLES=area_sanitaria_toledo SCALE=ORIGINAL COMPARE=area_sanitaria_toledo 
    PADJUST=BONFERRONI
    /MISSING CLASSMISSING=EXCLUDE
  /PRINT CPS DESCRIPTIVES MODELINFO FIT SUMMARY SOLUTION (EXPONENTIATED).

GENLIN defunciones BY area_sanitaria_toledo titgest (ORDER=ASCENDING) WITH inc_covid1865_1_ln 
    inc_covid1865_1_ln2
  /MODEL area_sanitaria_toledo titgest inc_covid1865_1_ln inc_covid1865_1_ln2 INTERCEPT=YES 
    SCALEWEIGHT=ln_plazas
 DISTRIBUTION=POISSON LINK=LOG
  /CRITERIA METHOD=FISHER(1) SCALE=1 COVB=MODEL MAXITERATIONS=100 MAXSTEPHALVING=5 
    PCONVERGE=1E-006(ABSOLUTE) SINGULAR=1E-012 ANALYSISTYPE=3(WALD) CILEVEL=95 CITYPE=WALD 
    LIKELIHOOD=FULL
  /EMMEANS TABLES=area_sanitaria_toledo SCALE=ORIGINAL COMPARE=area_sanitaria_toledo 
    PADJUST=BONFERRONI
  /EMMEANS TABLES=titgest SCALE=ORIGINAL COMPARE=titgest PADJUST=BONFERRONI
  /MISSING CLASSMISSING=EXCLUDE
  /PRINT CPS DESCRIPTIVES MODELINFO FIT SUMMARY SOLUTION (EXPONENTIATED).

USE ALL.
COMPUTE filter_$=(mvCOVIDdeaths = 0).
VARIABLE LABELS filter_$ 'mvCOVIDdeaths = 0 (FILTER)'.
VALUE LABELS filter_$ 0 'Not Selected' 1 'Selected'.
FORMATS filter_$ (f1.0).
FILTER BY filter_$.
EXECUTE.

GENLIN defunciones_COVID BY area_sanitaria_toledo (ORDER=ASCENDING) 
  /MODEL area_sanitaria_toledo  INTERCEPT=YES 
    SCALEWEIGHT=ln_plazas
 DISTRIBUTION=POISSON LINK=LOG
  /CRITERIA METHOD=FISHER(1) SCALE=1 COVB=MODEL MAXITERATIONS=100 MAXSTEPHALVING=5 
    PCONVERGE=1E-006(ABSOLUTE) SINGULAR=1E-012 ANALYSISTYPE=3(WALD) CILEVEL=95 CITYPE=WALD 
    LIKELIHOOD=FULL
  /EMMEANS TABLES=area_sanitaria_toledo SCALE=ORIGINAL COMPARE=area_sanitaria_toledo 
    PADJUST=BONFERRONI
    /MISSING CLASSMISSING=EXCLUDE
  /PRINT CPS DESCRIPTIVES MODELINFO FIT SUMMARY SOLUTION (EXPONENTIATED).

GENLIN defunciones_COVID BY area_sanitaria_toledo titgest (ORDER=ASCENDING) WITH inc_covid1865_1_ln 
    inc_covid1865_1_ln2
  /MODEL area_sanitaria_toledo titgest inc_covid1865_1_ln inc_covid1865_1_ln2 INTERCEPT=YES 
    SCALEWEIGHT=ln_plazas
 DISTRIBUTION=POISSON LINK=LOG
  /CRITERIA METHOD=FISHER(1) SCALE=1 COVB=MODEL MAXITERATIONS=100 MAXSTEPHALVING=5 
    PCONVERGE=1E-006(ABSOLUTE) SINGULAR=1E-012 ANALYSISTYPE=3(WALD) CILEVEL=95 CITYPE=WALD 
    LIKELIHOOD=FULL
  /EMMEANS TABLES=area_sanitaria_toledo SCALE=ORIGINAL COMPARE=area_sanitaria_toledo 
    PADJUST=BONFERRONI
  /EMMEANS TABLES=titgest SCALE=ORIGINAL COMPARE=titgest PADJUST=BONFERRONI
  /MISSING CLASSMISSING=EXCLUDE
  /PRINT CPS DESCRIPTIVES MODELINFO FIT SUMMARY SOLUTION (EXPONENTIATED).