# NOT RUN {
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')
# }
Run the code above in your browser using DataLab