Correction to mira code

rvlenth · May 28, 2024 · decccfe · decccfe
1 parent 3af99bb
commit decccfe
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/R/multiple-models.R b/R/multiple-models.R
@@ -177,7 +177,7 @@ emm_basis.mira = function(object, trms, xlev, grid, ...) {
              b = sum(a * (B %*% a))
              t = sum(a * (T %*% a))
              lambda = (1 + 1/m) * b / t
-             dfold = (m - 1)/lambda*2
+             dfold = (m - 1)/lambda^2
              dfobs = (dfcom + 1)/(dfcom + 3) * dfcom * (1 - lambda)
              ifelse(is.infinite(dfcom), dfold, 
                     dfold * dfobs/(dfold + dfobs))