compare.fits: Compare Fits
In rlme: Rank-Based Estimation and Prediction in Random Effects Nested Models

Description Usage Arguments Value See Also Examples

Compares two model fits. It returns tdbeta value and cfits values of two fits. The function uses the fixed effects estimates from fit 1 and fit 2 along with the covariance of the rank-based fit.

1	compare.fits(x, fit1, fit2)

`x`	Matrix of covariates
`fit1`	A class of type rlme.
`fit2`	A class of type rlme.

Returns tdbeta and cfits values.

fitdvcov

data(schools)
model = y ~ 1 + sex + age + (1 | region) + (1 | region:school)

# Extract covariants into matrix
cov = as.matrix(data.frame(schools[,"sex"], schools[,"age"]))

# Fit the models using each method
reml.fit = rlme(model, schools, method="reml")
gr.fit = rlme(model, schools, method="gr")

compare.fits(cov, reml.fit, gr.fit)