Convert General Estimating Equations (marginal means) from sas code to R code?

Question

I would like to convert the following SAS GEE model code to R. How can I do that for instance with the R packages (geepack OR gee)?

proc genmod data=inds descending; 
class SMA male disdur2 skill_isced id; 
model dayslost2 = SMA|disdur age male skill_isced / dist = gamma link=log; 
repeated subject=id / type=exch within=disdur2 ; 
title 'gee with exchangeable corr structure';
run;