A novel radio occultation (RO) technique, polarimetric RO (PRO), has recently been developed to measure differential polarimetric phase shift together with traditional RO products such as temperature and moisture. PRO observations have been shown to be associated with the vertical structure of cloud hydrometeors. With this unique measurement capability, the PRO soundings could potentially be used to evaluate model microphysics. This study compared PRO observations with WRF simulations of three typhoon cases in 2019 and 2021, initialized with ERA5 and NCEP FNL global analysis, respectively, with five microphysics parameterizations (Purdue Lin, WSM6, Goddard, Thompson, and Morrison). There is notable variability in the distribution of the model's hydrometeors, which could be affected by the initial conditions, microphysics parameterization schemes, typhoon locations, and circulation rainbands. The results show that WRF simulation initialized with ERA5 and using the Goddard microphysics scheme performs better in synoptic-scale verification and comparisons with PRO observations. The ensemble mean from 36 ensemble forecasts also exhibits consistent results with the deterministic run. The comparative results demonstrate that PRO data have the potential to evaluate the performance of different microphysics schemes in numerical models.