Monte Carlo studies of the critical behavior of site-dilute two-dimensional Ising models