c04 region-cluster variant
Source:R/otis_all_analyze.R
morie_otis_analyze_c04_ruhela_aggregate_region_cluster.Rdc04 region-cluster variant
Usage
morie_otis_analyze_c04_ruhela_aggregate_region_cluster(data, out_dir = NULL)
morie_otis_analyze_c04_mrm_aggregate_region_cluster(data, out_dir = NULL)Examples
# \donttest{
otis_c04 <- morie_synth_otis("c04", n = 120L, seed = 1L)
morie_otis_analyze_c04_ruhela_aggregate_region_cluster(otis_c04)
#> OTIS c04 -- Aggregate RF, region-clustered: Indigenous -> Individuals in RC
#> ===========================================================================
#> Source file c04 (region-clustered variant)
#> Dataset id c04-RC
#> Treatment column T_indigenous
#> Outcome (count) NumberIndividuals_RestrictiveConfinement
#> Year FE EndFiscalYear
#> Covariate FEs none
#> Cluster group (GEE) Region_MostRecentPlacement
#> Cells 120
#> Total outcome count 4395
#> Overdispersion (Poisson AIC - NB AIC) 1764.2
#>
#> Aggregate Ruhela formulation -- Poisson + NB IRR + GEE clustered for treatment effect on count outcome:
#> Family Coefficient IRR 95% CI p-value AIC Notes
#> --------------------------------------------------------------------------------------------------------------------------------------------
#> Poisson factor(T_indigenous)1 1.0171 [0.926, 1.117] 7.23e-01 2882.9 n=120, total_y=4395
#> NB factor(T_indigenous)1 1.0169 [0.620, 1.669] 9.47e-01 1118.6 n=120, total_y=4395
#> GEE-Poisson (cluster:Region_MostRecentPlacement, Exch) factor(T_indigenous)1 1.0183 [0.850, 1.220] 8.44e-01 -- groups=5, n_obs=120
#>
#> Region-clustered GEE variant of c04.
#>
#> 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.
# }