The WLW model can be implemented using the SAS PHREG procedure, selecting the STRATA option to allow different baseline hazards function for each HPV type. A robust variance is requested by choosing the option COVS(AGGREGATE).