HDtweedie_kfold_aug function

Fit a penalized gamma augmentation model via cross fitting