inst/extexamples/qdrg-biglm.R

In emmeans: Estimated Marginal Means, aka Least-Squares Means

# Post hoc analysis of a "biglm" object -- not supported by emmeans
bigmod <- biglm(log(conc) ~ source + factor(percent), data = pigs)

rg1 <- qdrg(log(conc) ~ source + factor(percent), data = pigs, 
    coef = coef(bigmod), vcov = vcov(bigmod), df = bigmod$df.residual)

emmeans(rg1, "source", type = "response")

## But in this particular case, we could have done it the easy way:
##     rg1 <- qdrg(object = bigmod)