R/eblup.mse.f.c1.R

eblup.mse.f.c1 <-
function(lme.obj, n.i, gamma.i, ...){
    var.e <- as.numeric(VarCorr(lme.obj)[,1])[2]
    ret <- gamma.i * (var.e / n.i)
    #class(ret) <- "eblup.mse.f"
    return(ret)
}

Try the JoSAE package in your browser

Any scripts or data that you put into this service are public.

JoSAE documentation built on May 2, 2019, 2:16 a.m.