avgPseudo {saePseudo}

Small Area Estimation using Averaging Pseudo Area Level Model

Description

Provides function for small area estimation at area level using averaging pseudo area level model for variables of interest. A dataset produced by data generation are also provided. This package estimates small areas at the village level and then aggregates them to the sub-district, region, and provincial levels.

Usage

avgPseudo(prov, reg, sub, vill, y, x, var, N, method = "REML")

Arguments

prov

Vector containing information of province

reg

Vector containing information of region

sub

Vector containing information of subdistrict

vill

Vector containing information of village

y

Direct estimation for each area

x

Auxiliary variable for each area

var

Sampling variances of direct estimators for each domain

N

Number of population in each area

method

Method used to fit the Fay-Herriot model, which can be either "ML", "REML" or "FH" methods. Default is "REML" method

Value

This function returns a list of the following objects:

Est_Area3

A dataframe with the values of Small Area Estimation with averaging pseudo area level model for sub-district level

Est_Area2

A dataframe with the values of Small Area Estimation with averaging pseudo area level model for region level

Est_Area1

A dataframe with the values of Small Area Estimation with averaging pseudo area level model for provincial level

Examples

# Load Dataset
data(dataVill)
saeAVG.Pseudo <- avgPseudo(prov = dataVill$Area1, reg = dataVill$Area2, sub = dataVill$Area3,
                 vill = dataVill$Area4, y = dataVill$ydir_area4, x = dataVill$X1,
                 var = dataVill$vardir_area4, N = dataVill$N, method="REML")

# Result
saeAVG.Pseudo$Est_Area3
saeAVG.Pseudo$Est_Area2
saeAVG.Pseudo$Est_Area1