RM_theta_eblup {agfh}

Traditional EBLUP Estimator of Theta

Description

Traditional EBLUP Estimator of Theta

Usage

  RM_theta_eblup(X, Y, D, theta.var.est)

Arguments

X

observed independent data to be analyzed

Y

observed dependent data to be analyzed

D

known precisions of response Y

theta.var.est

estimate of variance term for latent model; if NA, will automatically use method-of-moments

Details

Traditional EBLUP estimator of latent values theta.

Value

Returns a vector of estimates of theta.

Source

Marten Thompson [e-mail]

Examples

  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)