We designed this package to provide several functions for area level of small area estimation using hierarchical Bayesian (HB) method. This package provides model using Rao-Yu Model for variable interest.This package also provides a dataset produced by a data generation. The “rjags” package is employed to obtain parameter estimates. Model-based estimators involves the HB estimators which include the mean and the variation of mean. For the reference, see Rao and Molina (2015) and Torabi and Shokoohi (2012).
Velia Tri Marliana, Azka Ubaidillah
Velia Tri Marliana 221810642@stis.ac.id
RaoYuAr1()
This function gives estimation of y using
Hierarchical Bayesian under Rao Yu ModelPanel()
This function gives estimation of y using
Hierarchical Bayesian under Rao Yu Model when rho = 0You can install the development version of saeHB.panel from GitHub with:
# install.packages("devtools")
::install_github("Veliatrimarliana/saeHB.panel") devtools
This is a basic example which shows you how to solve a common problem:
library(saeHB.panel)
data(dataAr1)
= ydi ~ xdi1 + xdi2
formula = max(dataAr1[, "area"])
area = max(dataAr1[,"period"])
period = dataAr1[,4]
vardir <- Raoyu.Ar1(formula, area, period, vardir = vardir, data = dataAr1) result
Extract area mean estimation
$Est result
Extract coefficient estimation
$coefficient result
Extract area random effect variance
$refVar result
##References * Rao, J.N.K & Molina. (2015). Small Area Estimation 2nd Edition. New York: John Wiley and Sons, Inc. * Torabi, M., & Shokoohi, F. (2012). Likelihood inference in small area estimation by combining time-series and cross-sectional data. Journal of Multivariate Analysis, 111, 213–221. https://doi.org/10.1016/j.jmva.2012.05.016