Displaying 1 result from an estimated 1 matches for "resultreml".
Did you mean:
resultreg
2017 Mar 07
0
Potential clue for Bug 16975 - lme fixed sigma - inconsistent REML estimation
...ct, update, data)
object
}
###############################################################################
# Example bug report 16975
###############################################################################
data(milk)
milk$var<-milk$SD^2
#Fay-Herriot model using sae library
attach(milk)
resultREML <- eblupFH(yi ~ as.factor(MajorArea), SD^2)
resultREML
detach(milk)
# $fit
# $fit$method
# [1] "REML"
#
# $fit$convergence
# [1] TRUE
#
# $fit$iterations
# [1] 4
#
# $fit$estcoef
# beta std.error tvalue pvalue
# (Intercept) 0.9681890 0.06936208 13.958476 2.793443e-...