# NOT RUN {
## load dataset
data(datamsaeRBns)
# Compute EBLUP and Ratio Benchmark using auxiliary variables X1 and X2 for each dependent variable
## 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")
est_msae = est_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")]
est_msae = est_msaeRBns(Fo, vardir, weight, cluster)
## Return
est_msae$eblup$est.eblupRB # to see the Ratio Benchmark estimators
# }
Run the code above in your browser using DataLab