mse_FHme: Mean Squared Error Estimator of the EBLUP under a Fay-Herriot...

Description Usage Arguments Details Value Examples

View source: R/mse_FHme.R

Description

This function gives the mean squared error estimator of the EBLUP based on Fay-Herriot model with measurement error using jackknife method.

Usage

1
mse_FHme(formula, vardir, var.x, type.x = "witherror", data)

Arguments

formula

an object of class formula (or one that can be coerced to that class): a symbolic description of the model to be fitted. The variables included formula must have a length equal to the number of domains m. This formula can provide auxiliary variable either measured with error or without error or combination between them. If the auxiliary variable are combination between noerror and witherror variable, input all witherror variable first then noerror variable.

vardir

vector containing the m sampling variances of direct estimators for each domain. The values must be sorted as the Y.

var.x

vector containing mean squared error of X . The values must be sorted as the X. if you use optional data, input this parameter use c(""), example: var.x = c("c1") or var.x = c("c1","c2").

type.x

type of auxiliary variable used in the model. Either source measured with noerror, witherror and mix. Default value is witherror.

data

optional data frame containing the variables named in formula, vardir, and var.x.

Details

A formula has an implied intercept term. To remove this use either y ~ x - 1 or y ~ 0 + x. See formula for more details of allowed formulae.

Value

The function returns a list with the following objects:

mse

vector with the values of the mean squared errors of the EBLUPs for each domain.

Examples

1
2
3
4
5
6
data(dataME)
data(datamix)

mse.sae.me <- mse_FHme(formula = y ~ x.hat, vardir = vardir, var.x = c("var.x"), data = dataME)
mse.sae.mix <- mse_FHme(formula = y ~ x.hat1 + x.hat2 + x3 + x4,
                vardir = vardir, var.x = c("var.x1", "var.x2"), type.x = "mix", data = datamix)

saeME documentation built on Jan. 13, 2021, 11:03 a.m.