* Encoding: UTF-8.
* Fichero EPA4T2014-Segmentación.sav'.
* Los datos se refieren a la submuestra de asalariados: EMPRESA=1 OR EMPRESA=2.

FREQUENCIES CONTRATO JORNADA DURACION ANTIGUEDAD OCUPACION ACTIVIDAD EMPRESA EDUCACION
                        EDAD SEXO NACIONALIDAD.

*** Análisis de Correspondencias Múltiples.

* Obtención de todos los valores propios: m-p=40-8=32.
MULTIPLE CORRESPONDENCE VARIABLES=CONTRATO JORNADA DURACION ANTIGUEDAD OCUPACION ACTIVIDAD EMPRESA 
    EDUCACION EDAD SEXO NACIONALIDAD
  /ANALYSIS=CONTRATO(WEIGHT=1) JORNADA(WEIGHT=1) DURACION(WEIGHT=1) ANTIGUEDAD(WEIGHT=1) 
    OCUPACION(WEIGHT=1) ACTIVIDAD(WEIGHT=1) EMPRESA(WEIGHT=1) EDUCACION(WEIGHT=1) EDAD SEXO NACIONALIDAD    
  /MISSING=CONTRATO(PASSIVE,MODEIMPU) JORNADA(PASSIVE,MODEIMPU) DURACION(PASSIVE,MODEIMPU) 
    ANTIGUEDAD(PASSIVE,MODEIMPU) OCUPACION(PASSIVE,MODEIMPU) ACTIVIDAD(PASSIVE,MODEIMPU) 
    EMPRESA(PASSIVE,MODEIMPU) EDUCACION(PASSIVE,MODEIMPU) EDAD(PASSIVE,MODEIMPU) SEXO(PASSIVE,MODEIMPU) 
    NACIONALIDAD(PASSIVE,MODEIMPU) 
  /SUPPLEMENTARY=VARIABLE(EDAD SEXO NACIONALIDAD)
  /DIMENSION=32
  /NORMALIZATION=SYMMETRICAL
  /MAXITER=100
  /CRITITER=.00001
  /PRINT=NONE 
  /PLOT=NONE.

* ACM con dos factores.
MULTIPLE CORRESPONDENCE VARIABLES=CONTRATO JORNADA DURACION ANTIGUEDAD OCUPACION ACTIVIDAD EMPRESA 
    EDUCACION EDAD SEXO NACIONALIDAD
  /ANALYSIS=CONTRATO(WEIGHT=1) JORNADA(WEIGHT=1) DURACION(WEIGHT=1) ANTIGUEDAD(WEIGHT=1) 
    OCUPACION(WEIGHT=1) ACTIVIDAD(WEIGHT=1) EMPRESA(WEIGHT=1) EDUCACION(WEIGHT=1) EDAD SEXO NACIONALIDAD    
  /MISSING=CONTRATO(PASSIVE,MODEIMPU) JORNADA(PASSIVE,MODEIMPU) DURACION(PASSIVE,MODEIMPU) 
    ANTIGUEDAD(PASSIVE,MODEIMPU) OCUPACION(PASSIVE,MODEIMPU) ACTIVIDAD(PASSIVE,MODEIMPU) 
    EMPRESA(PASSIVE,MODEIMPU) EDUCACION(PASSIVE,MODEIMPU) EDAD(PASSIVE,MODEIMPU) SEXO(PASSIVE,MODEIMPU) 
    NACIONALIDAD(PASSIVE,MODEIMPU) 
  /SUPPLEMENTARY=VARIABLE(EDAD SEXO NACIONALIDAD)
  /DIMENSION=2
  /NORMALIZATION=SYMMETRICAL
  /MAXITER=100
  /CRITITER=.00001
  /PRINT=CORR DISCRIM OCORR QUANT(CONTRATO JORNADA DURACION ANTIGUEDAD OCUPACION ACTIVIDAD EMPRESA
    EDUCACION EDAD SEXO NACIONALIDAD)
  /PLOT=BIPLOT(20) OBJECT(20) JOINTCAT(CONTRATO JORNADA DURACION ANTIGUEDAD OCUPACION ACTIVIDAD
    EMPRESA EDUCACION EDAD SEXO NACIONALIDAD) (20) DISCRIM (20)
  /SAVE=OBJECT (Factor).

* Muestra del 15% al no poder procesar todos los casos..
USE ALL.
COMPUTE filter_$=(uniform(1)<=.15).
VARIABLE LABELS filter_$ 'Aproximadamente 30% de los casos (SAMPLE)'.
FORMATS filter_$ (f1.0).
FILTER  BY filter_$.
EXECUTE.

*** Análisis de clasificación: Método Ward.
CLUSTER   Factor1_1 Factor2_1
  /METHOD WARD
  /MEASURE=SEUCLID
  /PRINT SCHEDULE
  /PLOT NONE
  /SAVE CLUSTER(2,10).

FREQUENCIES CLU3_1.

* Generador de gráficos.
GGRAPH
  /GRAPHDATASET NAME="graphdataset" VARIABLES=Factor1_1 Factor2_1 CLU3_1 MISSING=LISTWISE 
    REPORTMISSING=NO
  /GRAPHSPEC SOURCE=INLINE
  /FITLINE TOTAL=NO SUBGROUP=NO.
BEGIN GPL
  SOURCE: s=userSource(id("graphdataset"))
  DATA: Factor1_1=col(source(s), name("Factor1_1"))
  DATA: Factor2_1=col(source(s), name("Factor2_1"))
  DATA: CLU3_1=col(source(s), name("CLU3_1"), unit.category())
  GUIDE: axis(dim(1), label("Dimensión de puntuaciones de objeto 1"))
  GUIDE: axis(dim(2), label("Dimensión de puntuaciones de objeto 2"))
  GUIDE: legend(aesthetic(aesthetic.color.interior), label("Ward Method                          ",
    "   "))
  GUIDE: text.title(label("Dispersión agrupada de Dimensión de puntuaciones de objeto 2 por ",
    "Dimensión de puntuaciones de objeto 1 por Ward Method"))
  ELEMENT: point(position(Factor1_1*Factor2_1), color.interior(CLU3_1))
END GPL.

IF CLU3_1=1 CLU3_1Media1=0.0886.
IF CLU3_1=1 CLU3_1Media2=0.5294.
IF CLU3_1=2 CLU3_1Media1=1.0441.
IF CLU3_1=2 CLU3_1Media2=-0.6494.	
IF CLU3_1=3 CLU3_1Media1=-1.0153.	
IF CLU3_1=3 CLU3_1Media2=-0.7616.
EXECUTE.

* Generador de gráficos.
GGRAPH
  /GRAPHDATASET NAME="graphdataset" VARIABLES=CLU3_1Media1 CLU3_1Media2 CLU3_1 MISSING=LISTWISE 
    REPORTMISSING=NO
  /GRAPHSPEC SOURCE=INLINE
  /FITLINE TOTAL=NO SUBGROUP=NO.
BEGIN GPL
  SOURCE: s=userSource(id("graphdataset"))
  DATA: CLU3_1Media1=col(source(s), name("CLU3_1Media1"))
  DATA: CLU3_1Media2=col(source(s), name("CLU3_1Media2"))
  DATA: CLU3_1=col(source(s), name("CLU3_1"), unit.category())
  GUIDE: axis(dim(1), label("CLU3_1Media1"))
  GUIDE: axis(dim(2), label("CLU3_1Media2"))
  GUIDE: legend(aesthetic(aesthetic.color.interior), label("Ward Method                          ",
    "   "))
  GUIDE: text.title(label("Dispersión agrupada de CLU3_1Media2 por CLU3_1Media1 por Ward Method"))
  ELEMENT: point(position(CLU3_1Media1*CLU3_1Media2), color.interior(CLU3_1))
END GPL.

CROSSTABS
  /TABLES=CONTRATO JORNADA DURACION ANTIGUEDAD OCUPACION ACTIVIDAD EMPRESA EDUCACION EDAD SEXO NACIONALIDAD 
                 BY CLU3_1
  /FORMAT=AVALUE TABLES
  /CELLS=COUNT ROW COLUMN 
  /COUNT ROUND CELL.

USE ALL.
***  Análisis de clasificación: Método de centros móviles.
QUICK CLUSTER Factor1_1 Factor2_1
  /MISSING=LISTWISE
  /CRITERIA=CLUSTER(3) MXITER(10) CONVERGE(0)
  /METHOD=KMEANS(NOUPDATE)
  /SAVE CLUSTER DISTANCE
  /PRINT INITIAL ANOVA.

***  Análisis de clasificación: Método bietápico.
TWOSTEP CLUSTER
  /CONTINUOUS VARIABLES=Factor1_1 Factor2_1
  /DISTANCE LIKELIHOOD
  /NUMCLUSTERS AUTO 15 BIC
  /NOSTANDARDIZE VARIABLES=Factor1_1 Factor2_1
  /HANDLENOISE 0
  /MEMALLOCATE 64
  /CRITERIA INITHRESHOLD(0) MXBRANCH(8) MXLEVEL(3)
  /VIEWMODEL DISPLAY=YES
  /SAVE VARIABLE=TSC_5865.