Produces estimates of domain mean or total along with the variance estimate based on estimated regression tree
domain_est(t1, domdat, weights = rep(1, nrow(domdat)), theta = "mean")
rpms object
data frame rows inside desired domain estimate
vector containing the sample weights for each observation
string of either "mean" or "total"
list containing F-statistic and corresponding p-value