simdata3: Dataset for Model 3

In mme: Multinomial Mixed Effects Models

Description

Dataset used by the multonomial mixed effects model with two independent random effects in each category of the response variable: one domain random effect and another correlated time and domain random effect (Model 3). This dataset contains ten small areas and four periods. The response variable has three categories. The last is the reference category. The variables are as follows:

Usage

simdata3

Format

A data frame with 40 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

## Not run: 

k=3 #number of categories of the response variable
pp=c(1,1) #vector with the number of auxiliary variables in each category
mod=3 #type of model
data(simdata3) #data
datar=data.mme(simdata3,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
msef=msef.ct(pp,datar$X,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)
## End(Not run)

Example output

Loading required package: MASS

 Package mme: Multinomial Mixed Effects Models 
 Version 0.1-6 (built on 2019-01-27) is now loaded.
 Copyright E. Lopez-Vizcaino, M.J. Lombardia and D. Morales 

mme documentation built on May 2, 2019, 10:46 a.m.