library(haven)
library(sjmisc)
library(sjlabelled)
library(tidyverse)
library(factoextra)
library(sjPlot)
library(sjstats)

db <- read_sav("/Users/robsalasco/odrive/Dropbox Personal/Daniel Roberto/Copia de SPSS ICHEM 1708 Casos FINAL CON PANEL y etiquetas PDN 2017 03 23.sav")

#labels
db$CL16 <- sjlabelled::set_label(db$CL16, "Sexo")
db$EDAD <- sjlabelled::set_label(db$EDAD, "Edad")
db$GSE_CALCULADO <- sjlabelled::set_label(db$GSE_CALCULADO, "Grupo Socioeconómico")

db$REF3 <- sjlabelled::set_label(db$REF3, "Tipos de solidaridad")


db$REF3 <- rec(db$REF3, rec = "8,9=NA; else=copy", val.labels = c("Solidaridad directa", 
                                                              "Solidaridad a través del Estado",
                                                              "Solidaridad a través de la comunidad",
                                                              "Individualismo"))

db$V6 <- sjlabelled::set_label(db$V6, "El Estado es...")


db$V6 <- rec(db$V6, rec = "8,9=NA; else=copy", val.labels = c("Estado democratico", 
                                                                  "Estado asistencialista",
                                                                  "Estado oligárquico"))
db$V6A <- sjlabelled::set_label(db$V6A, "El Estado debe ser...")

db$V6A <- rec(db$V6A, rec = "8,9=NA; else=copy", val.labels = c("Estado democratico", 
                                                              "Estado asistencialista",
                                                              "Estado oligárquico"))

db$RE3 <- sjlabelled::set_label(db$RE3, "Percepción democracia")

db$RE3<- rec(db$RE3, rec = "8,9=NA;1=1;2,3=2; else=copy", val.labels = c("Democracia necesaria", 
                                                                "Democracia suficiente"))

db$S7_A <- sjlabelled::set_label(db$S7_A, "Una oportunidad")

db$S7_A<- rec(db$S7_A, rec = "9=NA;1,2=1;3,4=2", val.labels = c("Deacuerdo", 
                                                                "Desacuerdo"))

db$S7_B <- sjlabelled::set_label(db$S7_B, "Una forma de abuso   ")

db$S7_B<- rec(db$S7_B, rec = "9=NA;1,2=1;3,4=2", val.labels = c("Deacuerdo", 
                                                                "Desacuerdo"))

db$CL1 <- sjlabelled::set_label(db$CL1, "Posición Política")

db$CL1<- rec(db$CL1, rec = "99,88=NA;1,2,3,4=1;5,6=2;7,8,9,10=3", val.labels = c("Izquierda","Centro",
                                                                "Derecha"))

db %>% select(REF3,CL16) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Tipos de solidaridad Sexo Total
Hombre Mujer
Solidaridad directa 176
29.1 %
175
26.2 %
351
27.6 %
Solidaridad a través
del Estado
210
34.7 %
235
35.2 %
445
35 %
Solidaridad a través
de la comunidad
88
14.5 %
98
14.7 %
186
14.6 %
Individualismo 131
21.7 %
159
23.8 %
290
22.8 %
Total 605
100 %
667
100 %
1272
100 %
χ2=2.115 · df=3 · Cramer’s V=0.042 · p=0.549
db %>% select(REF3,EDAD) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Tipos de solidaridad Edad Total
18-29 30-45 46-60 60++
Solidaridad directa 89
24.7 %
148
30.5 %
69
26.7 %
45
26.9 %
351
27.6 %
Solidaridad a través
del Estado
126
35 %
176
36.2 %
85
32.9 %
57
34.1 %
444
34.9 %
Solidaridad a través
de la comunidad
56
15.6 %
60
12.3 %
39
15.1 %
31
18.6 %
186
14.6 %
Individualismo 89
24.7 %
102
21 %
65
25.2 %
34
20.4 %
290
22.8 %
Total 360
100 %
486
100 %
258
100 %
167
100 %
1271
100 %
χ2=7.899 · df=9 · Cramer’s V=0.046 · p=0.544
db %>% select(REF3,GSE_CALCULADO) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Tipos de solidaridad Grupo Socioeconómico Total
ABC1 C2 C3 D E
Solidaridad directa 52
35.6 %
71
27.3 %
89
24.5 %
107
25.9 %
31
34.8 %
350
27.5 %
Solidaridad a través
del Estado
48
32.9 %
116
44.6 %
131
36.1 %
131
31.7 %
18
20.2 %
444
34.9 %
Solidaridad a través
de la comunidad
14
9.6 %
25
9.6 %
64
17.6 %
64
15.5 %
19
21.3 %
186
14.6 %
Individualismo 32
21.9 %
48
18.5 %
79
21.8 %
111
26.9 %
21
23.6 %
291
22.9 %
Total 146
100 %
260
100 %
363
100 %
413
100 %
89
100 %
1271
100 %
χ2=28.617 · df=12 · Cramer’s V=0.088 · p=0.004
db %>% select(V6,CL16) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
El Estado es… Sexo Total
Hombre Mujer
Estado democratico 321
43.3 %
339
43.8 %
660
43.6 %
Estado
asistencialista
65
8.8 %
92
11.9 %
157
10.4 %
Estado oligárquico 355
47.9 %
343
44.3 %
698
46.1 %
Total 741
100 %
774
100 %
1515
100 %
χ2=4.202 · df=2 · Cramer’s V=0.054 · p=0.122
db %>% select(V6,EDAD) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
El Estado es… Edad Total
18-29 30-45 46-60 60++
Estado democratico 174
40.8 %
237
42.7 %
154
48 %
96
45.1 %
661
43.6 %
Estado
asistencialista
35
8.2 %
65
11.7 %
32
10 %
25
11.7 %
157
10.4 %
Estado oligárquico 217
50.9 %
253
45.6 %
135
42.1 %
92
43.2 %
697
46 %
Total 426
100 %
555
100 %
321
100 %
213
100 %
1515
100 %
χ2=6.334 · df=6 · Cramer’s V=0.047 · p=0.387
db %>% select(V6,GSE_CALCULADO) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
El Estado es… Grupo Socioeconómico Total
ABC1 C2 C3 D E
Estado democratico 93
51.4 %
122
40.9 %
190
44.3 %
198
41.2 %
57
45.6 %
660
43.6 %
Estado
asistencialista
9
5 %
41
13.8 %
41
9.6 %
55
11.4 %
11
8.8 %
157
10.4 %
Estado oligárquico 79
43.6 %
135
45.3 %
198
46.2 %
228
47.4 %
57
45.6 %
697
46 %
Total 181
100 %
298
100 %
429
100 %
481
100 %
125
100 %
1514
100 %
χ2=18.363 · df=8 · Cramer’s V=0.080 · p=0.019
db %>% select(V6A,CL16) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
El Estado debe
ser…
Sexo Total
Hombre Mujer
Estado democratico 529
71.6 %
556
70.4 %
1085
71 %
Estado
asistencialista
141
19.1 %
181
22.9 %
322
21.1 %
Estado oligárquico 69
9.3 %
53
6.7 %
122
8 %
Total 739
100 %
790
100 %
1529
100 %
χ2=6.888 · df=2 · Cramer’s V=0.069 · p=0.032
db %>% select(V6A,EDAD) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
El Estado debe
ser…
Edad Total
18-29 30-45 46-60 60++
Estado democratico 321
73.8 %
401
71.2 %
225
70.1 %
139
65.9 %
1086
71 %
Estado
asistencialista
80
18.4 %
108
19.2 %
76
23.7 %
58
27.5 %
322
21 %
Estado oligárquico 34
7.8 %
54
9.6 %
20
6.2 %
14
6.6 %
122
8 %
Total 435
100 %
563
100 %
321
100 %
211
100 %
1530
100 %
χ2=11.457 · df=6 · Cramer’s V=0.063 · p=0.075
db %>% select(V6A,GSE_CALCULADO) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
El Estado debe
ser…
Grupo Socioeconómico Total
ABC1 C2 C3 D E
Estado democratico 146
78.5 %
202
68 %
309
71.4 %
340
70.7 %
89
66.4 %
1086
70.9 %
Estado
asistencialista
29
15.6 %
65
21.9 %
84
19.4 %
107
22.2 %
37
27.6 %
322
21 %
Estado oligárquico 11
5.9 %
30
10.1 %
40
9.2 %
34
7.1 %
8
6 %
123
8 %
Total 186
100 %
297
100 %
433
100 %
481
100 %
134
100 %
1531
100 %
χ2=17.366 · df=8 · Cramer’s V=0.077 · p=0.027
db %>% select(V6,CL1) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
El Estado es… Posición Política Total
Izquierda Centro Derecha
Estado democratico 93
31.2 %
173
40.9 %
160
64.3 %
426
43.9 %
Estado
asistencialista
22
7.4 %
43
10.2 %
35
14.1 %
100
10.3 %
Estado oligárquico 183
61.4 %
207
48.9 %
54
21.7 %
444
45.8 %
Total 298
100 %
423
100 %
249
100 %
970
100 %
χ2=71.261 · df=4 · Cramer’s V=0.202 · p=0.000
db %>% select(V6A,CL1) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
El Estado debe
ser…
Posición Política Total
Izquierda Centro Derecha
Estado democratico 224
75.4 %
318
74 %
165
66.3 %
707
72.4 %
Estado
asistencialista
66
22.2 %
91
21.2 %
47
18.9 %
204
20.9 %
Estado oligárquico 7
2.4 %
21
4.9 %
37
14.9 %
65
6.7 %
Total 297
100 %
430
100 %
249
100 %
976
100 %
χ2=18.019 · df=4 · Cramer’s V=0.101 · p=0.001
db %>% select(RE3,CL16) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Percepción
democracia
Sexo Total
Hombre Mujer
Democracia necesaria 545
71.3 %
523
63.2 %
1068
67.1 %
Democracia
suficiente
219
28.7 %
305
36.8 %
524
32.9 %
Total 764
100 %
828
100 %
1592
100 %
χ2=7.124 · df=1 · φ=0.068 · p=0.008
db %>% select(RE3,EDAD) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Percepción
democracia
Edad Total
18-29 30-45 46-60 60++
Democracia necesaria 275
63.1 %
405
68.8 %
237
70.5 %
152
65.5 %
1069
67.1 %
Democracia
suficiente
161
36.9 %
184
31.2 %
99
29.5 %
80
34.5 %
524
32.9 %
Total 436
100 %
589
100 %
336
100 %
232
100 %
1593
100 %
χ2=5.681 · df=3 · Cramer’s V=0.060 · p=0.128
db %>% select(RE3,GSE_CALCULADO) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Percepción
democracia
Grupo Socioeconómico Total
ABC1 C2 C3 D E
Democracia necesaria 134
74.9 %
216
71.5 %
282
63.7 %
342
66.7 %
95
60.9 %
1069
67.1 %
Democracia
suficiente
45
25.1 %
86
28.5 %
161
36.3 %
171
33.3 %
61
39.1 %
524
32.9 %
Total 179
100 %
302
100 %
443
100 %
513
100 %
156
100 %
1593
100 %
χ2=11.458 · df=4 · Cramer’s V=0.085 · p=0.022
credito <-db %>% select(S7_A,S7_B,CL16,EDAD,GSE_CALCULADO)
credito %>% select(S7_A,CL16) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Una oportunidad Sexo Total
Hombre Mujer
Deacuerdo 595
76.8 %
638
76.4 %
1233
76.6 %
Desacuerdo 180
23.2 %
197
23.6 %
377
23.4 %
Total 775
100 %
835
100 %
1610
100 %
χ2=0.169 · df=1 · φ=0.012 · p=0.681
credito %>% select(S7_A,EDAD) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Una oportunidad Edad Total
18-29 30-45 46-60 60++
Deacuerdo 365
79.2 %
450
77.2 %
255
75.2 %
164
71.9 %
1234
76.6 %
Desacuerdo 96
20.8 %
133
22.8 %
84
24.8 %
64
28.1 %
377
23.4 %
Total 461
100 %
583
100 %
339
100 %
228
100 %
1611
100 %
χ2=3.370 · df=3 · Cramer’s V=0.046 · p=0.338
credito %>% select(S7_A,GSE_CALCULADO) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Una oportunidad Grupo Socioeconómico Total
ABC1 C2 C3 D E
Deacuerdo 135
73.8 %
243
82.1 %
345
76.7 %
391
74.3 %
120
76.4 %
1234
76.6 %
Desacuerdo 48
26.2 %
53
17.9 %
105
23.3 %
135
25.7 %
37
23.6 %
378
23.4 %
Total 183
100 %
296
100 %
450
100 %
526
100 %
157
100 %
1612
100 %
χ2=0.985 · df=4 · Cramer’s V=0.025 · p=0.912
credito %>% select(S7_B,CL16) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Una forma de abuso Sexo Total
Hombre Mujer
Deacuerdo 622
82.2 %
677
82.7 %
1299
82.4 %
Desacuerdo 135
17.8 %
142
17.3 %
277
17.6 %
Total 757
100 %
819
100 %
1576
100 %
χ2=1.026 · df=1 · φ=0.028 · p=0.311
credito %>% select(S7_B,EDAD) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Una forma de abuso Edad Total
18-29 30-45 46-60 60++
Deacuerdo 373
82.3 %
472
83.1 %
277
82.9 %
177
80.5 %
1299
82.5 %
Desacuerdo 80
17.7 %
96
16.9 %
57
17.1 %
43
19.5 %
276
17.5 %
Total 453
100 %
568
100 %
334
100 %
220
100 %
1575
100 %
χ2=1.015 · df=3 · Cramer’s V=0.025 · p=0.798
credito %>% select(S7_B,GSE_CALCULADO) %>% sjtab(fun="xtab",show.col.prc = T,weight.by=db$PONDERADOR)
Una forma de abuso Grupo Socioeconómico Total
ABC1 C2 C3 D E
Deacuerdo 161
89 %
232
81.1 %
350
80.3 %
424
82.2 %
131
84 %
1298
82.4 %
Desacuerdo 20
11 %
54
18.9 %
86
19.7 %
92
17.8 %
25
16 %
277
17.6 %
Total 181
100 %
286
100 %
436
100 %
516
100 %
156
100 %
1575
100 %
χ2=5.386 · df=4 · Cramer’s V=0.059 · p=0.250