# NOT RUN {
##load dataset
data(datamsaeDB)
#Compute Fitted model for Y1, Y2, and Y3
#Y1 ~ X1 + X2
#Y2 ~ X2
#Y3 ~ X1
##Using parameter 'data'
formula = list(f1 = Y1~X1+X2,
f2 = Y2~X2,
f3 = Y3~X1)
vardir = c("v1","v12","v13","v2","v23","v3")
msaeFH <- msaefh(formula, vardir, data=datamsaeDB)
#Do not use parameter 'data'
formula = list(f1 = datamsaeDB$Y1~datamsaeDB$X1+datamsaeDB$X2,
f2 = datamsaeDB$Y2~datamsaeDB$X2,
f3 = datamsaeDB$Y3~datamsaeDB$X1)
vardir = datamsaeDB[,c("v1","v12","v13","v2","v23","v3")]
msaeFH_d <- msaefh(formula, vardir)
msaeFH$MSAE_Eblup #to see EBLUP Estimators
msaeFH$MSE_Eblup #to see estimated MSE of EBLUP estimators
# }
Run the code above in your browser using DataLab