We present a catalog of all CO (J = 4-3 through J = 13-12), [ C I], and [ N II] lines available from extragalactic spectra from the Herschel SPIRE Fourier Transform Spectrometer (FTS) archive combined with observations of the low-J CO lines from the literature and from the Arizona Radio Observatory. This work examines the relationships between L-FIR, L'(CO), and L-CO/L-CO,L-1-0. We also present a new method for estimating probability distribution functions from marginal signal-to-noise ratio Herschel FTS spectra, which takes into account the instrumental "ringing" and the resulting highly correlated nature of the spectra. The slopes of log(L-FIR) versus log (L'(CO)) are linear for all mid- to high-J CO lines and slightly sublinear if restricted to (ultra) luminous infrared galaxies ((U) LIRGs). The mid-to high-J CO luminosity relative to CO J - 1-0 increases with increasing L-FIR, indicating higher excitement of the molecular gas, although these ratios do not exceed similar to 180. For a given bin in L-FIR, the luminosities relative to CO J = 1-0 remain relatively flat from J = 6-5 through J = 13-12, across three orders of magnitude of L-FIR. A single component theoretical photodissociation region (PDR) model cannot match these flat SLED shapes, although combinations of PDR models with mechanical heating added qualitatively match the shapes, indicating the need for further comprehensive modeling of the excitation processes of warm molecular gas in nearby galaxies.