Glob observations w/ residuals above a certain magnitude

glob_obs_mat(X, r, thresh)

Arguments

X

design matrix

r

vector of residuals

thresh

threshold to use when computing globs