R code for GEE with bias corrected covariance for small sample is here.
The bias correction for variance of GEE estimates is performed based on the following reference: