# mse_FHme: Mean Squared Error Estimator of the EBLUP under a Fay-Herriot... In saeME: Small Area Estimation with Measurement Error

## 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.