# NOT RUN {
##load dataset
data(datamsaeDBns)
#Note : Make sure yout dataset does not contain NA Values
# you can set 0 in Direct estinations and vardir for non-sampled areas.
#Compute Fitted model for Y1, Y2, and Y3
#Y1 ~ X1 + X2
#Y2 ~ X1 + X2
#Y3 ~ X1 + X2
##Using parameter 'data'
formula = list(f1 = Y1~X1+X2,
f2 = Y2~X1+X2,
f3 = Y3~X1+X2)
vardir = c("v1","v12","v13","v2","v23","v3")
cluster = c("clY1","clY2","clY3")
nonsample = "nonsample"
saeFHns <- saefhns(formula,vardir,cluster, nonsample, data=datamsaeDBns)
#to calculate only one response variable
saeFHns1 <- saefhns(formula=list(f1=Y1~X1+X2),vardir ="v1", cluster = "clY1",
nonsample = "nonsample",data=datamsaeDBns )
# }
Run the code above in your browser using DataLab