In this work, we outline a new method for estimating dayside effective temperatures of exoplanets and associated uncertainties using Gaussian process (GP) regression. By applying our method to simulated observations, we show that the GP method estimates uncertainty more robustly than other model-independent approaches. We find that unbiased estimates of effective temperatures can be made using as few as three broad-band measurements (white-light HST WFC3 and the two warm Spitzer IRAC channels), although we caution that estimates made using only IRAC can be significantly biased. We then apply our GP method to the twelve hot Jupiters in the literature whose secondary eclipse depths have been measured by WFC3 and IRAC channels 1 and 2: CoRoT-2 b; HAT-P-7 b; HD 189733 b; HD 209458 b; Kepler-13A b; TrES-3 b; WASP-4 b; WASP-12 b; WASP-18 b; WASP-33 b; WASP-43 b; and WASP-103 b. We present model-independent dayside effective temperatures for these planets, with uncertainty estimates that range from $\pm$66 K to $\pm$136 K.