A method for determining luminosity functions incorporating both flux measurements and flux upper limits, with applications to the average X-ray to optical luminosity ratio for quasars