Climate data sets with both spatial and temporal components are often studied after removing from each time series a temporal mean calculated over a common reference interval, which is generally shorter than the overall length of the data set. The use of a short reference interval affects the temporal properties of the variability across the records, by reducing the standard deviation within the reference interval and inflating it elsewhere. For an annually averaged version of the Climate Research Unit’s (CRU) temperature anomaly product, the mean standard deviation is 0.67 C within the 1961-1990 reference interval, and 0.81 C elsewhere. The calculation of anomalies can be interpreted in terms of a two factor Analysis of Variance model. Within a Bayesian inference framework, any missing values are viewed as additional parameters, and the reference interval is specified as the full length of the data set. This Bayesian scheme is used to re-express the CRU data set as anomalies with respect to means calculated over the entire 1850-2009 interval spanned by the data set. The mean standard deviation is increased to 0.69 C within the original 1961-1990 reference interval, and reduced to 0.76 C elsewhere. The choice of reference interval thus has a predictable and demonstrable effect on the second spatial moment time series of the CRU data set. The spatial mean time series is in this case largely unaffected: the amplitude of spatial mean temperature change is reduced by 0.1 C when using the 1850-2009 reference interval, while the 90% uncertainty interval of (-0.03, 0.23) indicates that the reduction is not statistically significant.