We have discovered a correlation between the observed peak spectral energy E_(pk,obs) and the Euclidean value of of gamma-ray bursts (GRBs). We present the evidence for the correlation in the GUSBAD catalog and use it to derive the luminosity function of GRBs without using any redshifts. The procedure involves dividing GUSBAD GRBs into five spectral classes based on their E_(pk,obs). The overall luminosity function is derived assuming that each of the spectral classes contributes a Gaussian luminosity function. Their central luminosity is derived from the observed Euclidean . We explore various forms for the GRB rate function GR(z) in predicting redshift distributions of GRBs detected by Swift. We find that GR(z) peaks at a higher redshift than the typical star formation history currently favored in the literature. We consider two examples of GR(z) that successfully predict the observed redshift distribution of Swift GRBs. With the luminosity functions in hand, we convert the E_(pk,obs) – V/V_(max) correlation into an E_(pk,obs) – L_(iso) correlation and a rest-frame E_(pk) – L_(iso) correlation. By comparing the E_(pk) – L_(iso) correlation with a published correlation based on GRBs with known E_(pk,obs) and redshifts, we discuss the effect of Malmquist bias
