This dataset contains annual peak-flow data, PeakFQ specifications, and results of flood-frequency analyses of annual peak flows for 368 selected streamflow gaging stations (streamgages) operated by the U.S. Geological Survey (USGS) in the Great Lakes and Ohio River basins. "PeakFQinput_all.txt" contains annual peak-flow data, ending in water year 2013, for all 368 streamgages in the study area. Annual peak-flow data were obtained from the USGS National Water Information System (NWIS) database (https://nwis.waterdata.usgs.gov/usa/nwis/peak). "PeakFQspec_all.psf" contains PeakFQ specifications for all 368 streamgages in the study area. The specifications were developed by hydrologists in the various USGS Water Science Centers that participated in the study. "PeakFQoutput_all.PRT" contains the results of flood-frequency analyses of annual peak-flow data, for each of the 368 streamgages in the study area, that were conducted using the Expected Moments Algorithm (England and others, 2018). Using the annual peak-flow data in "PeakFQinput_all.txt" and the specifications in "PeakFQspec_all.psf", "PeakFQoutput_all.PRT" was generated in version 7.2 of USGS flood-frequency analysis software PeakFQ (https://water.usgs.gov/software/PeakFQ/; Veilleux and others, 2014). Results of the flood-frequency analyses were used to estimate regional skew for the study area using Bayesian Weighted Least Squares / Bayesian Generalized Least Squares (B-WLS / B-GLS) regression.