simdata {mme}

Dataset for Model 1

Description

Dataset used by the multinomial mixed effects model with one independent random effect in each category of the response variable (Model 1). This dataset contains 15 small areas. The response variable has three categories. The last is the reference category. The variables are as follows:

Usage

simdata

Format

A data frame with 15 rows and 9 variables in columns

Details

• area: area indicator.

• Time: time indicator.

• sample: the sample size of each domain.

• Population: the population size of each domain.

• Y1: the first category of the response variable.

• Y2: the second category of the response variable.

• Y3: the third category of the response variable.

• X1: the covariate for the first category of the response variable.

• X2: the covariate for the second category of the response variable.

Examples

k=3 #number of categories of the response variable
pp=c(1,1) #vector with the number of auxiliary variables in each category
data(simdata)  #data
mod=1 # type of model
datar=data.mme(simdata,k,pp,mod)
# Model fit
result=model(datar$d,datar$t,pp,datar$Xk,datar$X,datar$Z,datar$initial,datar$y[,1:(k-1)],
       datar$n,datar$N,mod)

#Analytic MSE
mse=msef(pp,datar$X,datar$Z,result,datar$N,datar$n)

B=1    #Bootstrap iterations
ss=12345 #SEED
set.seed(ss)

##Bootstrap parametric BIAS and MSE
mse.pboot=mseb(pp,datar$Xk,datar$X,datar$Z,datar$n,datar$N,result,B,mod)