powered by
Fiellers(model, Which = 1:length(which(b)), sample = "Sample", factor = paste0("factor(", sample, ")"), independent = "Z", df = summary(model)$df[2], alpha = 0.05)