powered by
Traditional EBLUP Estimator of Theta
RM_theta_eblup(X, Y, D, theta.var.est)
Returns a vector of estimates of theta.
observed independent data to be analyzed
observed dependent data to be analyzed
known precisions of response Y
Y
estimate of variance term for latent model; if NA, will automatically use method-of-moments
NA
Traditional EBLUP estimator of latent values theta.
X <- matrix(1:10, ncol=1) Y <- 2*X + rnorm(10, sd=1.1) D <- rep(1, 10) th.var.est <- 0.1 RM_theta_eblup(X, Y, D, th.var.est) RM_theta_eblup(X, Y, D)
Run the code above in your browser using DataLab