Extract Degrees of Freedom from Linear Hypotheses of Fixed Effects or LS Means.

Share:

Description

Determine degrees of freedom for custom linear hypotheses of fixed effects or LS Means using one of three possible approximation methods.

Usage

1
getDF(obj, L, method = c("contain", "residual", "satterthwaite"), ...)

Arguments

obj

(VCA) object

L

(matrix) specifying one or multiple linear hypothese, as returned by function getL

method

(character) the method to be used to determine the degrees of freedom for a linear hypothesis

...

additional parameters

Details

This is a convenience function to determine the DFs for linear hypotheses in the same way as function test.fixef. Only the "DF" part is returned here which can be passed to other functions expecting DFs as input.

Value

(numeric) vector with the DFs for each row of 'L'

Author(s)

Andre Schuetzenmeister andre.schuetzenmeister@roche.com

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
## Not run: 
data(VCAdata1)
tmpDat <- VCAdata1[VCAdata1$sample==1,]
tmpDat <- tmpDat[-c(11,51,73:76),]
fitMM <- anovaMM(y~(lot+device)/(day)/(run), tmpDat)
fitMM
L <- getL(fitMM, c("lot1-lot2", "device1-device2"))
getDF(fitMM, L)						# method="contain" is Default
getDF(fitMM, L, method="res")

getDF(fitMM, L, method="satt")		# takes quite long for this model

## End(Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.