cmseRHNERM function

Conditional mean squared error estimation of the empirical Bayes estimators under random heteroscedastic nested error regression models

Calculates the conditional mean squared error estimates of the empirical Bayes estimators under random heteroscedastic nested error regression models based on the parametric bootstrap.

cmseRHNERM(y, X, ni, C, k=1, maxr=100, B=100)

Arguments

• y: N*1 vector of response values.
• X: Np matrix containing N1 vector of 1 in the first column and vectors of covariates in the rest of columns.
• ni: m*1 vector of sample sizes in each area.
• C: m*p matrix of area-level covariates included in the area-level parameters.
• k: area number in which the conditional mean squared error estimator is calculated.
• maxr: maximum number of iteration for computing the maximum likelihood estimates.
• B: number of bootstrap replicates.

Returns

conditional mean squared error estimate in the kth area.

References

Kubokawa, K., Sugasawa, S., Ghosh, M. and Chaudhuri, S. (2016). Prediction in Heteroscedastic nested error regression models with random dispersions. Statistica Sinica, 26, 465-492.

Author(s)

Shonosuke Sugasawa

Examples

#generate data
set.seed(1234)
beta=c(1,1); la=1; tau=c(8,4)
m=20; ni=rep(3,m); N=sum(ni)
X=cbind(rep(1,N),rnorm(N))

mu=beta[1]+beta[2]*X[,2]
sig=1/rgamma(m,tau[1]/2,tau[2]/2); v=rnorm(m,0,sqrt(la*sig))
y=c()
cum=c(0,cumsum(ni))
for(i in 1:m){
  term=(cum[i]+1):cum[i+1]
  y[term]=mu[term]+v[i]+rnorm(ni[i],0,sqrt(sig[i]))
}

#fit the random heteroscedastic nested error regression
C=cbind(rep(1,m),rnorm(m))
cmse=cmseRHNERM(y,X,ni,C,B=10)
cmse