c01 region-cluster variant (year-clustered GEE).

Usage

morie_otis_analyze_c01_ruhela_aggregate_region_cluster(data, out_dir = NULL)

morie_otis_analyze_c01_mrm_aggregate_region_cluster(data, out_dir = NULL)

Arguments

data: c01 data.frame.
out_dir: Optional.

Value

morie_otis_analysis_result.

Examples

# \donttest{
otis_c01 <- morie_synth_otis("c01", n = 120L, seed = 1L)
morie_otis_analyze_c01_ruhela_aggregate_region_cluster(otis_c01)
#> OTIS c01 -- Aggregate RF, year-clustered: Female -> Individuals in RC
#> =====================================================================
#>   Source file                            c01 (year-clustered variant)
#>   Dataset id                             c01-RC
#>   Treatment column                       T_female
#>   Outcome (count)                        NumberIndividuals_RestrictiveConfinement
#>   Year FE                                EndFiscalYear
#>   Covariate FEs                          none
#>   Cluster group (GEE)                    EndFiscalYear
#>   Cells                                  120
#>   Total outcome count                    4579
#>   Overdispersion (Poisson AIC - NB AIC)  1481.8
#> 
#> Year-clustered GEE variant of c01. With only 3 fiscal years (= 3 clusters) the GEE inference is unreliable; shipped for transparency, not as the primary inference.
#> 
#> 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. 
# }