Compara, por sexo y banda, la varianza posterior de log(prevalencia) con la del ajuste solo con el ancla:
contracción c = 1 - Var_post / Var_prior. Cada celda recibe la etiqueta prior-driven, prior-informed o
data-driven; con varios valores de rho se queda la más conservadora. Las celdas subnacionales heredan la
etiqueta de su celda nacional.
Usage
dl_etiquetas(
ajuste,
ajuste_prior,
insumos,
grilla_rho = as.numeric(unlist(insumos$cfg$sensibilidad$rho)),
semilla,
opciones = NULL,
cascada = NULL
)Arguments
- ajuste
Ajuste de
dl_ajustar().- ajuste_prior
Ajuste de
dl_ajustar_solo_prior()de los mismos insumos.- insumos
Insumos de
dl_insumos().- grilla_rho
Valores de la correlación por edad del ancla (rho) con que se recalculan las etiquetas; por defecto, la grilla
sensibilidad.correlacion_edadde la configuración.- semilla
Semilla (entero), obligatoria: el resultado es reproducible con la misma semilla.
- opciones
Opciones de
dl_opciones_mcmc()de los reajustes por rho (por defecto, cadenas cortas: 2 de 6000 iteraciones con 3000 de calentamiento y adelgazamiento 10, sea cual sea el del ajuste, con el motor del ajuste; guardan 600 simulaciones, y de ellas se toman a lo sumo las del ajuste).- cascada
Cascada de
dl_cascada()(opcional): agrega las celdas subnacionales.
Value
Tabla (data.table) con una fila por ubicación, sexo y banda (location_id, year, age_group_id,
sex_id, cause_id): etiqueta, contraccion (c), lambda_usado (el peso del ancla), rho_usado (el rho
de la etiqueta publicada) y kappa_usado (el de la cascada en las celdas subnacionales; NA en las
nacionales).
Details
En cada celda (sexo y banda de edad) se comparan la contracción c y el desplazamiento d = |mediana_post(log q) - mediana_prior(log q)|, con q la prevalencia de la banda:
prior-driven(la estimación viene del ancla): d < 0.05 o c < 0.05;data-driven(los datos locales redujeron la varianza al menos un 30 %): c >= 0.3;prior-informed: en medio.
Sin datos locales en el ajuste todas las celdas son prior-driven. La contracción depende de la correlación por
edad del ancla (rho): con cada valor de grilla_rho distinto del de la configuración se hacen un ajuste y un
ajuste solo con el ancla nuevos (con opciones) y se publica la etiqueta más conservadora, con el rho que la dio.
Las celdas subnacionales no tienen parámetros propios: heredan la etiqueta de su celda nacional y declaran el
kappa de la cascada, para que ningún resultado subnacional se lea como evidencia local. Los tres valores de la
etiqueta son parte del contrato de las tablas de salida y no se traducen.
See also
dl_ajustar_solo_prior(), dl_sensibilidad() y dl_exportar_corrida() (que las guarda en la corrida).
Diagnóstico:
dl_ess(),
dl_rhat(),
dl_sensibilidad(),
dl_validar_ancla()
Examples
# \donttest{
b <- dl_insumos(dl_proyecto(dl_ejemplo(), causa = 9100))
#> dl_insumos(): 25 fila(s) nacionales de la tabla datos quedan fuera del ajuste porque datos_en_ajuste está vacío en la configuración; las filas subnacionales se usan solo para validar la cascada
op <- dl_opciones_mcmc(simulaciones = 100, cadenas = 2, iteraciones = 2000, calentamiento = 1000)
f <- dl_ajustar(b, op, semilla = 1)
f0 <- dl_ajustar_solo_prior(b, op, semilla = 1, ajuste = f)
# solo con el rho de la configuración (sin reajustes)
et <- dl_etiquetas(f, f0, b, grilla_rho = 0.5, semilla = 1)
head(et)
#> location_id year age_group_id sex_id cause_id etiqueta contraccion
#> <char> <int> <int> <int> <int> <char> <num>
#> 1: 123 2023 13 1 9100 prior-driven 0
#> 2: 123 2023 14 1 9100 prior-driven 0
#> 3: 123 2023 15 1 9100 prior-driven 0
#> 4: 123 2023 16 1 9100 prior-driven 0
#> 5: 123 2023 17 1 9100 prior-driven 0
#> 6: 123 2023 18 1 9100 prior-driven 0
#> lambda_usado rho_usado kappa_usado
#> <num> <num> <num>
#> 1: 1 0.5 NA
#> 2: 1 0.5 NA
#> 3: 1 0.5 NA
#> 4: 1 0.5 NA
#> 5: 1 0.5 NA
#> 6: 1 0.5 NA
# sin datos locales en el ajuste, todas las celdas vienen del ancla
table(et$etiqueta)
#>
#> prior-driven
#> 24
# }