This function wraps around several sandwich and lmtest functions to calculate robust standard errors and returns them in a useful format.

get_robust_se(model, type = "HC3", cluster = NULL,
  data = model.frame(model))

Arguments

model

A regression model, preferably of class lm or glm

type

One of "HC3", "const", "HC", "HC0", "HC1", "HC2", "HC4", "HC4m", "HC5". See sandwich::vcovHC() for some more details on these choices. Note that some of these do not work for clustered standard errors (see sandwich::vcovCL()]).

cluster

If you want clustered standard errors, either a string naming the column in data that represents the clusters or a vector of clusters that is the same length as the number of rows in data.

data

The data used to fit the model. Default is to just get the model.frame from model.

Value

A list with the following:

  • coefs: a coefficient table with the estimates, standard errors, t-statistics, and p-values from lmtest.

  • ses: The standard errors from coefs.

  • ts: The t-statistics from coefs.

  • ps: The p-values from coefs.

  • type: The argument to robust.

  • use_cluster: TRUE or FALSE indicator of whether clusters were used.

  • cluster: The clusters or name of cluster variable used, if any.

  • vcov: The robust variance-covariance matrix.