R/eblup.mse.f.gamma.i.R

Defines functions eblup.mse.f.gamma.i

Documented in eblup.mse.f.gamma.i

eblup.mse.f.gamma.i <-
function(lme.obj, n.i,...){
        #extract variance components: first for random, second for residual
    varcomp <- as.numeric(VarCorr(lme.obj)[,1])
    ret <- varcomp[1]/(varcomp[1]+varcomp[2]/n.i)#checked and ok
    #class(ret) <- "JoSAE.gamma.i.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 6, 2018, 5:04 p.m.