robust: Computes a coefficient summary table for glms with robust...

In bmckuw/UWbe536: UW BIOST/EPI 536 Utility Functions

Description

Computes a coefficient summary table for glms with robust standard errors and robust CIs

Usage

robust(model, type = c("HC1"), df = Inf, expo = TRUE, level = 0.95,
  digits = 3)

Arguments

model	The fit of a generalized linear model (lm or glm)
type	Character string giving type of robust standard error from sandwich package
df	Degrees of freedom for reference t distribution (Inf means Normal distribution)
expo	Logical expression indicating whether exponentiated linear combination of coefficients is requiested
level	Numeric confidence coefficient for CIs (as proportion)
digits	Number of digits to round to in output

Value

A table with coefficient estimates (exponentiated or not), CIs and P-values based on robust standard errors

Examples

x <- rbinom(100, 2, .5)
expit <- function(x) exp(x)/(1 + exp(x))
y <- rbinom(100, 1, p = expit(x + .1))
model <- glm(y~x, family = "binomial")
robust(model)

bmckuw/UWbe536 documentation built on May 31, 2019, 1:54 p.m.