Sensibilidad

Veamos un poco de descriptiva de la base de datos:

      table1( ~ .|grupo,data=df
      , render.continuous=c(.="Mean (SD) ", " Mediana [Q1,Q3] "= "Q2 [Q1, Q3]",
                          "CV%, [Rango] "="CV%, [MIN, MAX]"))

## Warning in table1.formula(~. | grupo, data = df, render.continuous = c(.
## = "Mean (SD) ", : Terms to the right of '|' in formula 'x' define table
## columns and are expected to be factors with meaningful labels.

	0 (n=41)	1 (n=39)	Overall (n=80)
grupo
Mean (SD)	0.00 (0.00)	1.00 (0.00)	0.488 (0.503)
Mediana [Q1,Q3]	0.00 [0.00, 0.00]	1.00 [1.00, 1.00]	0.00 [0.00, 1.00]
CV%, [Rango]	NA	0%, [1.00, 1.00]	103.2%, [0.00, 1.00]
edad
Mean (SD)	27.9 (4.24)	41.8 (10.5)	34.7 (10.5)
Mediana [Q1,Q3]	27.0 [25.0, 29.1]	40.5 [34.3, 50.3]	31.2 [26.3, 40.7]
CV%, [Rango]	15.2%, [23.0, 43.3]	25.2%, [20.4, 63.8]	30.4%, [20.4, 63.8]
ESTATURA
Mean (SD)	1.65 (0.0491)	1.65 (0.0610)	1.65 (0.0549)
Mediana [Q1,Q3]	1.64 [1.61, 1.69]	1.64 [1.61, 1.69]	1.64 [1.61, 1.69]
CV%, [Rango]	3.0%, [1.57, 1.77]	3.7%, [1.54, 1.78]	3.3%, [1.54, 1.78]
PESO
Mean (SD)	60.1 (7.01)	59.6 (6.96)	59.8 (6.94)
Mediana [Q1,Q3]	59.0 [55.0, 62.0]	60.0 [56.0, 64.0]	59.0 [55.0, 64.0]
CV%, [Rango]	11.7%, [49.0, 75.0]	11.7%, [41.0, 73.5]	11.6%, [41.0, 75.0]
IMC
Mean (SD)	36.4 (4.06)	36.1 (4.17)	36.3 (4.09)
Mediana [Q1,Q3]	35.4 [33.7, 39.2]	36.0 [34.5, 39.3]	35.8 [33.9, 39.2]
CV%, [Rango]	11.2%, [29.3, 46.5]	11.6%, [26.1, 46.5]	11.3%, [26.1, 46.5]
PIE
Mean (SD)	0.135 (0.187)	0.0369 (0.0280)	0.0871 (0.143)
Mediana [Q1,Q3]	0.0700 [0.0400, 0.160]	0.0400 [0.0200, 0.0400]	0.0400 [0.0200, 0.0700]
CV%, [Rango]	138.3%, [0.0200, 1.00]	75.8%, [0.00800, 0.160]	164.0%, [0.00800, 1.00]
MANO
Mean (SD)	0.0279 (0.0259)	0.0448 (0.158)	0.0361 (0.111)
Mediana [Q1,Q3]	0.0200 [0.0200, 0.0200]	0.0200 [0.00800, 0.0200]	0.0200 [0.00800, 0.0200]
CV%, [Rango]	92.8%, [0.00800, 0.160]	352.1%, [0.00800, 1.00]	308.0%, [0.00800, 1.00]
ANOD
Mean (SD)	0.0248 (0.0211)	0.0343 (0.0353)	0.0294 (0.0291)
Mediana [Q1,Q3]	0.0200 [0.00800, 0.0200]	0.0200 [0.00800, 0.0400]	0.0200 [0.00800, 0.0400]
CV%, [Rango]	85.1%, [0.00600, 0.0700]	103.1%, [0.00800, 0.160]	99.0%, [0.00600, 0.160]
ANAI
Mean (SD)	0.0247 (0.0219)	0.0234 (0.0278)	0.0241 (0.0248)
Mediana [Q1,Q3]	0.0200 [0.00800, 0.0400]	0.0200 [0.00800, 0.0200]	0.0200 [0.00800, 0.0400]
CV%, [Rango]	88.8%, [0.00800, 0.0700]	118.7%, [0.00800, 0.160]	103.1%, [0.00800, 0.160]
VAGD
Mean (SD)	0.0212 (0.0268)	0.0171 (0.0139)	0.0192 (0.0215)
Mediana [Q1,Q3]	0.00800 [0.00800, 0.0200]	0.00800 [0.00800, 0.0200]	0.00800 [0.00800, 0.0200]
CV%, [Rango]	126.7%, [0.00800, 0.160]	81.4%, [0.00800, 0.0700]	112.0%, [0.00800, 0.160]
VAGI
Mean (SD)	0.0170 (0.0147)	0.0169 (0.0133)	0.0169 (0.0139)
Mediana [Q1,Q3]	0.00800 [0.00800, 0.0200]	0.00800 [0.00800, 0.0200]	0.00800 [0.00800, 0.0200]
CV%, [Rango]	86.3%, [0.00800, 0.0700]	79.0%, [0.00800, 0.0700]	82.3%, [0.00800, 0.0700]
CLIT
Mean (SD)	0.0221 (0.0185)	0.0174 (0.0119)	0.0199 (0.0157)
Mediana [Q1,Q3]	0.0200 [0.00800, 0.0200]	0.00800 [0.00800, 0.0200]	0.0200 [0.00800, 0.0200]
CV%, [Rango]	83.6%, [0.00800, 0.0700]	68.5%, [0.00800, 0.0400]	79.3%, [0.00800, 0.0700]
PARTOS
Mean (SD)	NA	0.564 (0.502)	0.564 (0.502)
Mediana [Q1,Q3]	NA	1.00 [0.00, 1.00]	1.00 [0.00, 1.00]
CV%, [Rango]	NA	89.1%, [0.00, 1.00]	89.1%, [0.00, 1.00]
Missing	41 (100%)	0 (0%)	41 (51.2%)
TIPO
Mean (SD)	NA	1.03 (1.11)	1.03 (1.11)
Mediana [Q1,Q3]	NA	1.00 [0.00, 2.00]	1.00 [0.00, 2.00]
CV%, [Rango]	NA	108.4%, [0.00, 3.00]	108.4%, [0.00, 3.00]
Missing	41 (100%)	0 (0%)	41 (51.2%)
EPISIOTOMIA
Mean (SD)	NA	0.333 (0.478)	0.333 (0.478)
Mediana [Q1,Q3]	NA	0.00 [0.00, 1.00]	0.00 [0.00, 1.00]
CV%, [Rango]	NA	143.3%, [0.00, 1.00]	143.3%, [0.00, 1.00]
Missing	41 (100%)	0 (0%)	41 (51.2%)
CIRUGIA
0	0 (0%)	15 (38.5%)	15 (18.8%)
CESAREA-LEGRADO	0 (0%)	1 (2.6%)	1 (1.2%)
CESAREA-LIBERACION PIRAMIDAL	0 (0%)	1 (2.6%)	1 (1.2%)
CIRUGIA CANCER RECTO	0 (0%)	1 (2.6%)	1 (1.2%)
ESTRABISMOS	0 (0%)	1 (2.6%)	1 (1.2%)
FISURA ANAL,HEMORROIDES	0 (0%)	1 (2.6%)	1 (1.2%)
HEMORROIDECTOMIA-ESFINTEROTOMIA-QISTE UTERO-CESAREA	0 (0%)	1 (2.6%)	1 (1.2%)
HERNIA INGUINAL	0 (0%)	1 (2.6%)	1 (1.2%)
HISTERECTOMIA	0 (0%)	2 (5.1%)	2 (2.5%)
LAPAROSCOPIA	0 (0%)	2 (5.1%)	2 (2.5%)
LEGRADO-LAPAROSCOPIA	0 (0%)	2 (5.1%)	2 (2.5%)
LIPOSUCION-ABDOMINAL	0 (0%)	1 (2.6%)	1 (1.2%)
MALLAS IU	0 (0%)	1 (2.6%)	1 (1.2%)
PECHO-IVE	0 (0%)	1 (2.6%)	1 (1.2%)
PUDENDOS-LAPAROSCOPIA	0 (0%)	1 (2.6%)	1 (1.2%)
QUISTE-RECTAL	0 (0%)	1 (2.6%)	1 (1.2%)
RTU-RESECCION URETRAL	0 (0%)	1 (2.6%)	1 (1.2%)
SUTURAPERINEAL	0 (0%)	1 (2.6%)	1 (1.2%)
TABIQUE-VAGINAL	0 (0%)	1 (2.6%)	1 (1.2%)
URETER	0 (0%)	1 (2.6%)	1 (1.2%)
Missing	41 (100%)	2 (5.1%)	43 (53.8%)
REGLA
Mean (SD)	NA	0.270 (0.450)	0.270 (0.450)
Mediana [Q1,Q3]	NA	0.00 [0.00, 1.00]	0.00 [0.00, 1.00]
CV%, [Rango]	NA	166.6%, [0.00, 1.00]	166.6%, [0.00, 1.00]
Missing	41 (100%)	2 (5.1%)	43 (53.8%)
ENFERMEDAD_ACTUAL
COXIGODINEA	0 (0%)	1 (2.6%)	1 (1.2%)
DIPAREUNIA	0 (0%)	1 (2.6%)	1 (1.2%)
DISPAREUNIA	0 (0%)	11 (28.2%)	11 (13.8%)
DISPAREUNIA-MIGRAÑAS	0 (0%)	2 (5.1%)	2 (2.5%)
DISPAREUNIA SUPERFICIAL	0 (0%)	1 (2.6%)	1 (1.2%)
DISPREUNIA	0 (0%)	1 (2.6%)	1 (1.2%)
DOLOR	0 (0%)	6 (15.4%)	6 (7.5%)
DOLOR ANAL	0 (0%)	2 (5.1%)	2 (2.5%)
DOLOR CICATRIZ	0 (0%)	1 (2.6%)	1 (1.2%)
DOLOR DEFECATORIO	0 (0%)	1 (2.6%)	1 (1.2%)
DOLOR VESICAL	0 (0%)	1 (2.6%)	1 (1.2%)
DPC	0 (0%)	4 (10.3%)	4 (5.0%)
ENDOMETRIOSIS	0 (0%)	1 (2.6%)	1 (1.2%)
ENDOMETRIOSIS-DISMENORREA PRIEMARIA	0 (0%)	1 (2.6%)	1 (1.2%)
ESTENOSIS VAGINAL-INCONTINENCIA RECTAL	0 (0%)	1 (2.6%)	1 (1.2%)
IU-HIPERTONIA VAGINAL-RECTAL	0 (0%)	1 (2.6%)	1 (1.2%)
NEURALGIA PUDENDO	0 (0%)	1 (2.6%)	1 (1.2%)
VULVODINIA	0 (0%)	2 (5.1%)	2 (2.5%)
Missing	41 (100%)	0 (0%)	41 (51.2%)
MEDICAMENTOS
0	0 (0%)	16 (41.0%)	16 (20.0%)
1	0 (0%)	13 (33.3%)	13 (16.2%)
AINE OCASIONAL	0 (0%)	1 (2.6%)	1 (1.2%)
AMITRIPTILINA-TRAMADOL	0 (0%)	1 (2.6%)	1 (1.2%)
ANALGESICOS-AINES	0 (0%)	1 (2.6%)	1 (1.2%)
BLISSEL-GEL-ESTRÓGENOS	0 (0%)	1 (2.6%)	1 (1.2%)
CIDINE, NEXIUM	0 (0%)	1 (2.6%)	1 (1.2%)
ESTROGENOS	0 (0%)	1 (2.6%)	1 (1.2%)
HIDROFEROL,LORAZEPAM,VITC,GINGO BILOBA,ALPRAZOLAM,EXCITALOPRAM	0 (0%)	1 (2.6%)	1 (1.2%)
Missing	41 (100%)	3 (7.7%)	44 (55.0%)
IU
0	0 (0%)	32 (82.1%)	32 (40.0%)
1	0 (0%)	7 (17.9%)	7 (8.8%)
Missing	41 (100%)	0 (0%)	41 (51.2%)

Esos IMC on ridículos. Solo habéis mirado gordas? No me creo que solo hayáis cogido a mujeres que no perdonan una merienda.

En fin, vamos a ponerles algunos “p” a la comparación:

      rndr <- function(x, name, ...) {
        if (length(x) == 0) {
          y <- df[[name]]
          s <- rep("", length(render.default(x=y, name=name, ...)))
          if (is.numeric(y)) {
            p <- t.test(y ~ df$grupo)$p.value
          } else {
            p <- chisq.test(table(y, droplevels(df$grupo)))$p.value
          }
          s[2] <- sub("<", "&lt;", format.pval(p, digits=3, eps=0.001))
          s
        } else {
          render.default(x=x, name=name, ...)
        }
      }
      
      rndr.strat <- function(label, n, ...) {
        ifelse(n==0, label, render.strat.default(label, n, ...))
      }
      
      
      df$grupo=factor(df$grupo,levels=c(0,1,2),labels=c("0","1","p-valor"))  
      table1( ~ edad+ESTATURA+PESO+IMC+PIE+MANO+ANOD+ANAI+VAGD+VAGI+CLIT|grupo,
             data=df, droplevels=F, render=rndr, render.strat=rndr.strat, overall=F)

	0 (n=41)	1 (n=39)	p-valor
edad
Mean (SD)	27.9 (4.24)	41.8 (10.5)	<0.001
Median [Min, Max]	27.0 [23.0, 43.3]	40.5 [20.4, 63.8]
ESTATURA
Mean (SD)	1.65 (0.0491)	1.65 (0.0610)	0.999
Median [Min, Max]	1.64 [1.57, 1.77]	1.64 [1.54, 1.78]
PESO
Mean (SD)	60.1 (7.01)	59.6 (6.96)	0.74
Median [Min, Max]	59.0 [49.0, 75.0]	60.0 [41.0, 73.5]
IMC
Mean (SD)	36.4 (4.06)	36.1 (4.17)	0.746
Median [Min, Max]	35.4 [29.3, 46.5]	36.0 [26.1, 46.5]
PIE
Mean (SD)	0.135 (0.187)	0.0369 (0.0280)	0.00186
Median [Min, Max]	0.0700 [0.0200, 1.00]	0.0400 [0.00800, 0.160]
MANO
Mean (SD)	0.0279 (0.0259)	0.0448 (0.158)	0.512
Median [Min, Max]	0.0200 [0.00800, 0.160]	0.0200 [0.00800, 1.00]
ANOD
Mean (SD)	0.0248 (0.0211)	0.0343 (0.0353)	0.155
Median [Min, Max]	0.0200 [0.00600, 0.0700]	0.0200 [0.00800, 0.160]
ANAI
Mean (SD)	0.0247 (0.0219)	0.0234 (0.0278)	0.818
Median [Min, Max]	0.0200 [0.00800, 0.0700]	0.0200 [0.00800, 0.160]
VAGD
Mean (SD)	0.0212 (0.0268)	0.0171 (0.0139)	0.391
Median [Min, Max]	0.00800 [0.00800, 0.160]	0.00800 [0.00800, 0.0700]
VAGI
Mean (SD)	0.0170 (0.0147)	0.0169 (0.0133)	0.974
Median [Min, Max]	0.00800 [0.00800, 0.0700]	0.00800 [0.00800, 0.0700]
CLIT
Mean (SD)	0.0221 (0.0185)	0.0174 (0.0119)	0.179
Median [Min, Max]	0.0200 [0.00800, 0.0700]	0.00800 [0.00800, 0.0400]

Sensibilidad

31/12/2019