# NOT RUN {
## load dataset
data(datamsaeRBns)
# Compute EBLUP and Ratio Benchmark using auxiliary variables X1 and X2 for each dependent variable
## Using parameter 'data'
est_sae = est_saeRBns(Y1 ~ X1 + X2, v1, w1, c1, data = datamsaeRBns)
## Without parameter 'data'
est_sae = est_saeRBns(datamsaeRBns$Y1 ~ datamsaeRBns$X1 + datamsaeRBns$X2,
datamsaeRBns$v1, datamsaeRBns$w1, datamsaeRBns$c1)
## Return
est_sae$eblup$est.eblupRB # to see the Ratio Benchmark estimators
# }
Run the code above in your browser using DataLab