Concatenating data from the millimetre regime to the infrared, we have performed spectral energy distribution (SED) modelling for 227 of the 405 millimetre continuum sources of Hill et al. (2005, Cat. J/MNRAS/363/405) which are thought to contain young massive stars in the earliest stages of their formation. Three main parameters are extracted from the fits: temperature, mass and luminosity. The method employed was the Bayesian inference, which allows a statistically probable range of suitable values for each parameter to be drawn for each individual protostellar candidate. This is the first application of this method to massive star formation.