c12 aggregate Ruhela: Female -> median RC days (by region)

Usage

morie_otis_analyze_c12_ruhela_aggregate(data, out_dir = NULL)

morie_otis_analyze_c12_mrm_aggregate(data, out_dir = NULL)

Arguments

data: c12 data.frame.
out_dir: Optional.

Value

morie_otis_analysis_result.

Examples

# \donttest{
otis_c12 <- morie_synth_otis("c12", n = 120L, seed = 1L)
morie_otis_analyze_c12_ruhela_aggregate(otis_c12)
#> OTIS c12 -- Aggregate RF: Female -> median RC days, by region
#> =============================================================
#>   Source file                            c12_segregation_and_restrictive_confinement_aggregate_durations_by_region.csv
#>   Dataset id                             c12
#>   Treatment column                       T_female
#>   Outcome (count)                        TotalAggregatedDays_RestrictiveConfinement
#>   Year FE                                EndFiscalYear
#>   Covariate FEs                          Region_MostRecentPlacement
#>   Cluster group (GEE)                    --
#>   Cells                                  37
#>   Total outcome count                    1467
#>   Overdispersion (Poisson AIC - NB AIC)  275.9
#> 
#> Aggregate Ruhela formulation -- Poisson + NB IRR for treatment effect on count outcome:
#>   Family   Coefficient        IRR     95% CI          p-value   AIC    Notes
#>   ---------------------------------------------------------------------------------------
#>   Poisson  factor(T_female)1  1.0199  [0.886, 1.174]  7.84e-01  634.4  n=37, total_y=1467
#>   NB       factor(T_female)1  1.0324  [0.579, 1.839]  9.14e-01  358.5  n=37, total_y=1467
#> 
#> Aggregate RF on c12: region-level companion to c10.
#> 
#> IRR > 1 ==> treatment increases the count rate; IRR < 1 ==> treatment decreases it. Concordance between Poisson and NB indicates equidispersion; large gap (Poisson AIC > NB AIC) indicates overdispersion -- trust NB. 
# }