scale.value = scale.value))
W <- geemod$naive.variance
fname <-
- if is.function(family) deparse(substitute(family)) else family
+ if (is.function(family)) deparse(substitute(family)) else family
if (fname == "binomial")
W <- summary(glm(formula, family = quasibinomial, data = data))$cov.scaled
N <- geemod$nobs