# NOT RUN {
## load dataset
data(datamsaeRBns)
# Compute MSE EBLUP and Ratio Benchmark
# This is the long running example
## Using parameter 'data'
Fo = list(f1 = Y1 ~ X1 + X2,
f2 = Y2 ~ X1 + X2,
f3 = Y3 ~ X1 + X2)
vardir = c("v1", "v12", "v13", "v2", "v23", "v3")
weight = c("w1", "w2", "w3")
cluster = c("c1", "c2", "c3")
mse_msae = mse_msaeRBns(Fo, vardir, weight, cluster, data = datamsaeRBns)
## Without parameter 'data'
Fo = list(f1 = datamsaeRBns$Y1 ~ datamsaeRBns$X1 + datamsaeRBns$X2,
f2 = datamsaeRBns$Y2 ~ datamsaeRBns$X1 + datamsaeRBns$X2,
f3 = datamsaeRBns$Y3 ~ datamsaeRBns$X1 + datamsaeRBns$X2)
vardir = datamsaeRBns[, c("v1", "v12", "v13", "v2", "v23", "v3")]
weight = datamsaeRBns[, c("w1", "w2", "w3")]
cluster = datamsaeRBns[, c("c1", "c2", "c3")]
mse_msae = mse_msaeRBns(Fo, vardir, weight, cluster)
## Return
mse_msae$pbmse.eblupRB # to see the MSE of Ratio Benchmark
# }
Run the code above in your browser using DataLab