R implementation of generalized estimating equations