reml {qle} | R Documentation |
Compute the value of the REML function (without constant term) for a given covariance model and data
reml(models, pars, data, Xs, verbose = FALSE)
models |
object of class |
pars |
covariance parameter vector (including global scalar nugget value) |
data |
data frame of simulated statistics, each column corresponds to a
single covariance model in ' |
Xs |
matrix of sample points |
verbose |
logical, if |
Given a list of covariance models the function calculates the REML function values at
the covariance parameter 'pars
'.
List of REML function values.
M. Baaske
data(normal) # extract the sample points (the design) X <- as.matrix(qsd$qldata[1:2]) # get the observed statistic T <- qsd$qldata[c("mean.T1")] reml(qsd$covT[1],pars=c(1e-4,1.5,0.1),T,X)