When observations are excluded from a covariance adjustment model fit (as tracked in stats::na.action(SandwichLayer@fitted_covariance_model)), the following code produces a vector of ID's that doesn't align with the row dimension of SandwichLayer@keys:
|
C_cls <- Reduce( |
|
function(l, r) paste(l, r, sep = "_"), |
|
as.list(stats::expand.model.frame(cmod, cluster)[, cluster, drop=FALSE]) |
|
) |
|
} |
|
jk_units <- unique(C_cls[in_Q]) |
This leads to some units being jackknifed that need not be.
When observations are excluded from a covariance adjustment model fit (as tracked in
stats::na.action(SandwichLayer@fitted_covariance_model)), the following code produces a vector of ID's that doesn't align with the row dimension ofSandwichLayer@keys:propertee/R/teeMod.R
Lines 489 to 494 in 2b726c0
This leads to some units being jackknifed that need not be.