JoSAE (version 0.1)

JoSAE.gamma.i.f: Returns the gamma_i value for the GREG-like design of the EBLUP estimator

Description

The gamma_i value is also necessary for the variance of the EBLUP estimate

Usage

JoSAE.gamma.i.f(...)

## S3 method for class 'default':
JoSAE.gamma.i.f(lme.obj, n.i, ...)

Arguments

lme.obj

a linear mixed-effects model generated with lme

n.i

the number of samples within domain i

...

further attributes forwarded to other functions; not used so far

Value

• The gamma_i value

See Also

JoSAE-package for more examples

Examples

library(nlme)
data(JoSAE.sample.data)
#fit a lme 
summary(fit.lme <- lme(biomass.ha ~ mean.canopy.ht, data=JoSAE.sample.data
                       , random=~1|domain.ID))
#calculate gamma_i for a domain with 5 samples
JoSAE.gamma.i.f(fit.lme, 5)