Hola Ramón llegó todo muy bien y ya pude trabajar con los cuadros que te envio para que veas como van quedando revisá los títulos y las unidades que puse en cada uno. Están faltando los que te puse los títulos en rojo que son :
1-prevalencia de lesiones según consumo de tabaco ý género para muestra Int-Mont
load("~/Dropbox/odontologia/relevamiento/casnati/muestra_global/muestra_global.RData")
library(survey)
##
## Attaching package: 'survey'
##
## The following object(s) are masked from 'package:graphics':
##
## dotchart
tabla12 <- svyby(~Grupos.lesiones, ~sexo, disenio.post, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla12[, 2:7] * 100, 1)
## Grupos.lesiones1-cancerizables Grupos.lesiones2-candidosis
## F 0.3 3.4
## M 0.8 3.1
## Grupos.lesiones3-traumatica Grupos.lesiones4-proliferativa-T. benigno
## F 3.2 3.1
## M 2.7 1.4
## Grupos.lesiones5-infecciosa-inmunitaria Grupos.lesiones6-Sin lesiones
## F 0.9 89.1
## M 0.5 91.4
tabla13 <- svyby(~Grupos.lesiones, ~tabaco, disenio.post, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla13[, 2:7] * 100, 1)
## Grupos.lesiones1-cancerizables Grupos.lesiones2-candidosis
## 1 0.5 3.5
## 2 0.0 0.9
## 3 1.8 2.4
## 4 0.4 3.0
## Grupos.lesiones3-traumatica Grupos.lesiones4-proliferativa-T. benigno
## 1 3.2 2.3
## 2 1.9 0.0
## 3 3.0 3.8
## 4 2.0 1.7
## Grupos.lesiones5-infecciosa-inmunitaria Grupos.lesiones6-Sin lesiones
## 1 0.9 89.6
## 2 0.0 97.3
## 3 0.0 89.0
## 4 0.2 92.7
2- Prevalencia de lesiones según consumo de alcohol y género para muestra Int-Mont
tabla14 <- svyby(~Grupos.lesiones, ~sexo * alcohol, disenio.post, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla14[, 3:8] * 100, 1)
## Grupos.lesiones1-cancerizables Grupos.lesiones2-candidosis
## F.1 0.3 4.9
## M.1 1.2 4.8
## F.2 0.3 2.4
## M.2 0.5 1.6
## F.3 0.0 2.0
## M.3 1.5 5.5
## F.4 0.0 1.7
## M.4 0.0 2.2
## Grupos.lesiones3-traumatica Grupos.lesiones4-proliferativa-T. benigno
## F.1 3.5 4.6
## M.1 3.2 2.7
## F.2 2.4 2.3
## M.2 2.7 1.9
## F.3 8.1 0.8
## M.3 1.8 0.0
## F.4 0.0 2.8
## M.4 4.8 0.0
## Grupos.lesiones5-infecciosa-inmunitaria Grupos.lesiones6-Sin lesiones
## F.1 0.5 86.1
## M.1 0.0 88.1
## F.2 0.4 92.3
## M.2 1.1 92.4
## F.3 5.8 83.2
## M.3 0.0 91.3
## F.4 0.0 95.5
## M.4 0.0 93.0
3- Prevalencia de lesiones según perfil de consumo de frutas y verduras para muestra Int Mont
tabla15 <- svyby(~frutas_verduras, ~Grupos.lesiones, disenio.post, svymean,
na.rm = TRUE, deff = TRUE)
round(tabla15[, 2:4] * 1, 1)
## frutas_verduras se DEff.frutas_verduras
## 1-cancerizables 2.3 0.3 1.1
## 2-candidosis 1.9 0.2 1.1
## 3-traumatica 2.6 0.4 1.6
## 4-proliferativa-T. benigno 3.0 0.3 1.8
## 5-infecciosa-inmunitaria 2.0 0.1 0.5
## 6-Sin lesiones 2.3 0.1 1.5
round(confint(tabla15), 1)
## 2.5 % 97.5 %
## 1-cancerizables 1.8 2.8
## 2-candidosis 1.6 2.2
## 3-traumatica 1.9 3.4
## 4-proliferativa-T. benigno 2.4 3.6
## 5-infecciosa-inmunitaria 1.7 2.2
## 6-Sin lesiones 2.2 2.4
plot(table(round(disenio.post$variables$frutas_verduras, 2)), ylab = "Frecuencia",
xlab = "Porciones diarias")
svyhist(~frutas_verduras, disenio.post, main = "Consumo de Fritas y verduras",
col = "yellow", cex.main = 1, xlab = "Porciones diarias", ylab = "Frecuencia")
abline(v = svymean(~frutas_verduras, disenio.post, na.rm = TRUE), col = 2)
abline(v = 5, col = 4)
svyboxplot(frutas_verduras ~ Grupos.lesiones, disenio.post, all.outliers = TRUE,
horizontal = TRUE, main = "Consumo de Fritas y verduras", col = "yellow",
cex.main = 1, xlab = "Porciones diarias", ylab = "Lesiones")
confint(tabla15)
## 2.5 % 97.5 %
## 1-cancerizables 1.766 2.816
## 2-candidosis 1.572 2.197
## 3-traumatica 1.858 3.414
## 4-proliferativa-T. benigno 2.402 3.637
## 5-infecciosa-inmunitaria 1.722 2.199
## 6-Sin lesiones 2.166 2.375
4- Prevalencia de lesiones según INSE para muestra Int-Mont
tabla16 <- svyby(~Grupos.lesiones, ~INSE1, disenio.post, svymean, na.rm = TRUE,
deff = TRUE)
round(tabla16[, 3:8] * 100, 1)
## Grupos.lesiones2-candidosis Grupos.lesiones3-traumatica
## 1-BAJO 4.8 3.5
## 2-MEDIO 2.2 2.7
## 3-ALTO 0.0 0.5
## Grupos.lesiones4-proliferativa-T. benigno
## 1-BAJO 2.9
## 2-MEDIO 1.9
## 3-ALTO 0.0
## Grupos.lesiones5-infecciosa-inmunitaria
## 1-BAJO 1.3
## 2-MEDIO 0.3
## 3-ALTO 0.0
## Grupos.lesiones6-Sin lesiones se.Grupos.lesiones1-cancerizables
## 1-BAJO 87.0 0.3
## 2-MEDIO 92.4 0.2
## 3-ALTO 99.5 0.0
round(confint(tabla16) * 100, 1)
## 2.5 % 97.5 %
## 1-BAJO:Grupos.lesiones1-cancerizables -0.1 1.3
## 2-MEDIO:Grupos.lesiones1-cancerizables 0.1 1.0
## 3-ALTO:Grupos.lesiones1-cancerizables 0.0 0.0
## 1-BAJO:Grupos.lesiones2-candidosis 2.6 6.9
## 2-MEDIO:Grupos.lesiones2-candidosis 0.9 3.5
## 3-ALTO:Grupos.lesiones2-candidosis 0.0 0.0
## 1-BAJO:Grupos.lesiones3-traumatica 1.6 5.4
## 2-MEDIO:Grupos.lesiones3-traumatica 1.3 4.1
## 3-ALTO:Grupos.lesiones3-traumatica -0.5 1.6
## 1-BAJO:Grupos.lesiones4-proliferativa-T. benigno 1.0 4.8
## 2-MEDIO:Grupos.lesiones4-proliferativa-T. benigno 0.4 3.3
## 3-ALTO:Grupos.lesiones4-proliferativa-T. benigno 0.0 0.0
## 1-BAJO:Grupos.lesiones5-infecciosa-inmunitaria -0.2 2.8
## 2-MEDIO:Grupos.lesiones5-infecciosa-inmunitaria 0.0 0.6
## 3-ALTO:Grupos.lesiones5-infecciosa-inmunitaria 0.0 0.0
## 1-BAJO:Grupos.lesiones6-Sin lesiones 83.4 90.5
## 2-MEDIO:Grupos.lesiones6-Sin lesiones 90.0 94.7
## 3-ALTO:Grupos.lesiones6-Sin lesiones 98.4 100.5
5- La prevalencia de las lesiones según la topografía me la enviaste solo con los valores expandidos y recuerdas que tiene que ir en frecuencia relativa porcentuada y tiene que ser para la muestra global y la Int- Mont
tabla11 <- svyby(~muestra, ~Topografia.rec, disenio.post, svytotal, na.rm = TRUE,
deff = TRUE)
round(tabla11[, 2:3] * 1, 0)
## muestraInt muestraMon
## Bordes 10895 0
## Labios 10472 10098
## Lengua 1686 947
## Mucosa 4090 6858
## Paladar 16975 8250
## Rebordes 3116 11519
## Sin trastornos 356998 435338
round(confint(tabla11) * 1, 1)
## 2.5 % 97.5 %
## Bordes:muestraInt 4604.1 17185
## Labios:muestraInt 1910.1 19035
## Lengua:muestraInt 2.8 3370
## Mucosa:muestraInt 1027.5 7153
## Paladar:muestraInt 10496.7 23453
## Rebordes:muestraInt -149.4 6382
## Sin trastornos:muestraInt 344582.3 369414
## Bordes:muestraMon 0.0 0
## Labios:muestraMon 4093.8 16102
## Lengua:muestraMon -910.9 2805
## Mucosa:muestraMon 1179.3 12536
## Paladar:muestraMon 2275.5 14224
## Rebordes:muestraMon 3926.0 19112
## Sin trastornos:muestraMon 422853.3 447823
tabla16 <- svyby(~muestra, ~Trastornos.rec, disenio.post, svytotal, na.rm = TRUE,
deff = FALSE)
round(ftable(tabla16), 0)
## muestraInt muestraMon
## Trastornos.rec
## absceso svytotal 743 4070
## SE 442 2880
## candidosis svytotal 20558 7747
## SE 4242 2942
## hiperplasia fibrom svytotal 2800 4211
## SE 1277 2003
## hiperplasia parapr svytotal 6340 297
## SE 3157 297
## leucoplasia svytotal 2834 1654
## SE 1108 1234
## otras svytotal 6158 3249
## SE 3105 2049
## queratosis fricc svytotal 960 3836
## SE 732 1996
## Sin trastornos svytotal 356998 435338
## SE 6335 6370
## ulcera svytotal 6842 12606
## SE 2146 3510
addmargins(table(disenio.post$variables$Trastornos.rec, disenio.post$variables$muestra))
##
## Int Mon Sum
## absceso 3 2 5
## candidosis 40 11 51
## hiperplasia fibrom 6 5 11
## hiperplasia parapr 12 1 13
## leucoplasia 8 2 10
## otras 13 3 16
## queratosis fricc 2 5 7
## Sin trastornos 820 520 1340
## ulcera 18 14 32
## Sum 922 563 1485
Enfermedad peri
tabla17 <- svyby(~Prev.lesio, ~perio, disenio.post, svymean, na.rm = TRUE, deff = FALSE)
round(ftable(tabla17) * 100, 1)
## Prev.lesiono Prev.lesiosi
## perio
## 1-sano svymean 94.0 6.0
## SE 1.1 1.1
## 2-leve svymean 97.6 2.4
## SE 1.4 1.4
## 3-moderada svymean 89.8 10.2
## SE 3.9 3.9
## 4-grave svymean 95.0 5.0
## SE 2.5 2.5
## 5-desdentada svymean 76.0 24.0
## SE 3.2 3.2
svychisq(~Prev.lesio + perio, disenio.post, statistic = "Wald")
##
## Design-based Wald test of association
##
## data: svychisq(~Prev.lesio + perio, disenio.post, statistic = "Wald")
## F = 8.42, ndf = 4, ddf = 1475, p-value = 1.025e-06