knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

STEP 1: Load package and load the data

library(saeHB.zinb)
data("dataZINB")

STEP 2: Fitting HB Model

result <- ZinbHB(formula = y ~ x1 + x2, data = dataZINB)

STEP 3 Extract mean estimation

Small Area mean Estimates

result$Est 

Coefficient Estimation

result$coefficient  

Random effect variance estimation

result$refVar 


hayunbuto/saeHB.zinb documentation built on June 22, 2022, 12:02 p.m.