Regional regression equations were calculated in Puerto Rico with generalized least squares techniques to estimate flood frequency statistics at ungaged locations using drainage area as the only explanatory variable. The island was divided into 2 regions to minimize residuals. The region division that resulted in lower and more balanced residuals runs primarily north-south near the center of the island, mostly along an 8-digit hydrologic unit code (HUC8) boundary. The division line runs through a HUC8 polygon on the southern end of the island, but care was taken to include entire watersheds and consideration was given where hydrologic and physiographic properties differed. This data release includes geographic information system files that define the polygons for both regions.