proc genmod data=skinny ;
      class personid ;
      model  sample1_totalSpermCount_1 = samplePerson_1_byr
            deepgen1935c totalchildren_1935 /  dist=normal link = identity;
      repeated  subject=personid / type=exch;
      where deepgen_1935>=3;
   run;

skinny.gee3 <- gee(sample1_totalSpermCount_1 ~ samplePerson_1_byr + deepgen1935.c + totalchildren_1935,
                   id = PersonID,
                   data = clean[clean$DeepGen_1935 >= 3, ],
                   family = gaussian,
                   corstr = "exchangeable")

Coefficients:
                       Estimate   Naive S.E.    Naive z  Robust S.E.   Robust z
(Intercept)        1376.0792841 1659.6117221  0.8291574 1906.1089688  0.7219311
samplePerson_1_byr   -0.5671856    0.8434850 -0.6724312    0.9694925 -0.5850335
deepgen1935.c       -12.4639418    7.0401018 -1.7704207    4.3967397 -2.8348146
totalchildren_1935    0.6572002    0.6458183  1.0176240    0.4139229  1.5877359