Package 'zipsae'

Sample Data for Small Area Estimation with Zero-Inflated Poisson model

Description

A Dataset which is generate with Zero-Inflated Poisson method for Small Area Estimation purpose

This data is generated based on Zero-Inflated Poisson with EBLUP based model

Usage

dataSAEZIP

Format

A data frame with 300 rows and 3 variables:

y

Direct Estimation of y

x1

Auxiliary variable of x1

vardir

Sampling Variance of y

---

EBLUPs under Zero-Inflated Poisson Model

Description

This function produce EBLUPs for Zero-Inflated data and its Relative Standard Error. Small Area Estimation with Zero-Inflated Model (SAE-ZIP) is a model developed for Zero-Inflated data that can lead us to overdispersion situation. To handle this kind of situation, this model is created. The model in this package is based on Small Area Estimation with Zero-Inflated Poisson model proposed by Dian Christien Arisona (2018). For the data sample itself, we use combination method between Roberto Benavent and Domingo Morales (2015) and Sabine Krieg, Harm Jan Boonstra and Marc Smeets (2016).

Usage

zipsae(data, vardir, formula, PRECISION = 1e-04, MAXITER = 100)

Arguments

data	The data frame with vardir, response, and explanatory variables included with Zero-Inflated situation also.
vardir	Sampling variances of direct estimations, if it is included in data frame so it is the vector with the name of sampling variances.if it is not, it is a data frame of sampling variance in order : var1, cov12,.,cov1r,var2,cov23,.,cov2r,.,cov(r-1)(r),var(r)
formula	List of formula that describe the fitted model
PRECISION	Limit of Fisher-scoring convergence tolerance. We set the default in 1e-4
MAXITER	Maximum number of iterations in Fisher-scoring algorithm. We set the default in 100

Value

This function returns a list of the following objects:

estimate	A Vector with a list of EBLUP with Zero-Inflated Poisson model
dispersion	A list containing the following objects:

• rse : A dataframe with the values of relative square errors of estimation

coefficient

A list containing the following objects:

• a : EBLUP estimator for Non-Zero data

• b : EBLUP estimator for Complete Data

Examples

##load the dataset in package
data(dataSAEZIP)

##Extract the vardir (sampling error)
dataSAEZIP$vardir -> sError

##Compute the data with SAE ZIP model
formula = (y~x1)
zipsae(data = dataSAEZIP, vardir = sError, formula) -> saezip

saezip$estimate        #to see the result of Small Area Estimation with Zero-Inflated Model
saezip$dispersion$rse  #to see the relative standard error from the estimation
saezip$coefficient$a   #to see the a coefficient which is gained from the non-zero compilation data
saezip$coefficient$b   #to see the b coefficient which is gained from the complete compilation data.

head(saezip)

---