c11 aggregate Ruhela: long-duration bin (>=16 days) -> RC

Usage

morie_otis_analyze_c11_ruhela_aggregate(data, out_dir = NULL)

morie_otis_analyze_c11_mrm_aggregate(data, out_dir = NULL)

Arguments

data: c11 data.frame.
out_dir: Optional.

Value

morie_otis_analysis_result.

Examples

# \donttest{
otis_c11 <- morie_synth_otis("c11", n = 120L, seed = 1L)
morie_otis_analyze_c11_ruhela_aggregate(otis_c11)
#> OTIS c11 -- Aggregate RF: Long-duration bin (>=16 days) -> Individuals in RC
#> ============================================================================
#>   Source file                            c11_individuals_in_segregation_and_restrictive_confinement_aggregate_lengths.csv
#>   Dataset id                             c11
#>   Treatment column                       T_long_duration
#>   Outcome (count)                        NumberIndividuals_RestrictiveConfinement
#>   Year FE                                EndFiscalYear
#>   Covariate FEs                          none
#>   Cluster group (GEE)                    --
#>   Cells                                  120
#>   Total outcome count                    5101
#>   Overdispersion (Poisson AIC - NB AIC)  1121.4
#> 
#> Aggregate Ruhela formulation -- Poisson + NB IRR for treatment effect on count outcome:
#>   Family   Coefficient               IRR     95% CI          p-value   AIC     Notes
#>   ------------------------------------------------------------------------------------------------
#>   Poisson  factor(T_long_duration)1  1.0113  [0.954, 1.072]  7.04e-01  2241.4  n=120, total_y=5101
#>   NB       factor(T_long_duration)1  1.0307  [0.805, 1.319]  8.10e-01  1120    n=120, total_y=5101
#> 
#> Aggregate RF on c11 testing whether long-duration bins (>=16 days) concentrate more individuals in RC.
#> 
#> 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. 
# }