minque: Minimum norm quadratic unbiased estimation
In varComp: Variance Component Models

Description Usage Arguments Value Author(s) References See Also Examples

This function computes the minimum norm quadratic unbiased estimate of variance components. Typically this is used as starting values of REML.

1	minque(y, varcov, start = rep(0, length(k)), lower.bound = -Inf, restricted = TRUE)

`y`	A numeric vector of zero-mean response variable.
`varcov`	A list of variance-covariance matrices, the same as in `varComp` with no `random` argument.
`start`	A numeric vector of prior values of the ratio of variance components to the error variance.
`lower.bound`	A numeric value of the ratio of variance components to the error variance. If this is -Inf (default), minque solves the linear equation as in Rao (1972). If this is zero, quadratic programming is used to minimize the squared error between the two sides of the linear equation.
`restricted`	Not used. Currently only restricted likelihood is supported.

A numeric vector of estimates of ratio of variance components to error variance.

Long Qu

C. Radhakrishna Rao (1972) Estimation of Variance and Covariance Components in Linear Models. Journal of the American Statistical Association 67, 112–115.

varComp.fit

library(nlme)
data(Oxide)
vcf0 = varComp(Thickness~Source, Oxide, ~Lot/Wafer, 
      control=varComp.control(nlminb=nlminb.control(iter.max=0L)))
coef(vcf0,'var.ratio')
(st=minque(vcf0$residual.contrast, vcf0$working.cor, lower.bound = 0))
(vcf = varComp(Thickness~Source, Oxide, ~Lot/Wafer, control=varComp.control(start=st)))
coef(vcf,'var.ratio')