| datos | R Documentation |
A simulated dataset created by the authors in order to check the correct operation of the saery package
data(datos)
A simulated data frame with 2000 observations on the following 6 variables.
Areaa numeric vector with the area (domain) of the data
Perioda numeric vector with the period (subdomain) of the data
ydia numeric vector with the direct estimator of the indicator of interest for area (domain)
onesa numeric vector. This is only needed to include the intercept parameter in the model
xdia numeric vector containing the aggregated (population) values of an auxiliary variable
sigma2edia numeric vector with the known variance of the error term
sigma2edi <- datos[,6]
X <- as.matrix(datos[,5])
ydi <- datos[,3]
D <- length(unique(datos[,1]))
md <- rep(length(unique(datos[,2])), D)
output.fit.ar1 <- fit.saery(X, ydi, D, md, sigma2edi, "AR", 0.9)
output.fit.ar1
#For computational reasons B is too low. We recomend to increase up to 100
eblup.output.ar1 <- eblup.saery(X, ydi, D, md, sigma2edi, "a", plot = TRUE, B = 2)
eblup.output.ar1
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.