Fits penalized generalized estimating equations to longitudinal data with high-dimensional covariates.