We studied Lyman-α (Lyα) escape in a statistical sample of 43 Green Peas with HST/COS Lyα spectra. Green Peas are nearby star-forming galaxies with strong [OIII]λ5007 emission lines. Our sample is four times larger than the previous sample and covers a much more complete range of Green Pea properties. We found that about 2/3 of Green Peas are strong Lyα line emitters with rest-frame Lyα equivalent width > 20Å. The Lyα profiles of Green Peas are diverse. The Lyα escape fraction, defined as the ratio of observed Lyα flux to intrinsic Lyα flux, shows anti-correlations with a few Lyα kinematic features -both the blue peak and red peak velocities, the peak separations, and FWHM of the red portion of the Lyα profile. Using properties measured from SDSS optical spectra, we found many correlations -Lyα escape fraction generally increases at lower dust reddening, lower metallicity, lower stellar mass, and higher [OIII]/[OII] ratio. We fit their Lyα profiles with the HI shell radiative transfer model and found Lyα escape fraction anti-correlates with the best-fit N HI . Finally, we fit an empirical linear relation to predict f Lyα esc from the dust extinction and Lyα red peak velocity. The standard deviation of this relation is about 0.3 dex. This relation can be used to isolate the effect of IGM scatterings from Lyα escape and to probe the IGM optical depth along the line of sight of each z > 7 Lyα emission line galaxy in the JWST era.