data("grapes", package = "sae")
data("grapesprox", package = "sae")
fitRFH <- rfh(
grapehect ~ area + workdays - 1,
data = grapes,
samplingVar = "var"
)
# The variance component of a mixed linear model:
matV <- variance(fitRFH)
# The full variance matrix:
matV$V()
# The sampling error component
matV$Ve()
# the random effects component
matV$Vu()
Run the code above in your browser using DataLab