Objective. Numerical simulations are largely adopted to estimate dosimetric quantities, e.g. specific absorption rate (SAR) and temperature increase, in tissues to assess the patient exposure to the radiofrequency (RF) field generated during magnetic resonance imaging (MRI). Simulations rely on reference anatomical human models and tabulated data of electromagnetic and thermal properties of biological tissues. However, concerns may arise about the applicability of the computed results to any phenotype, introducing a significant degree of freedom in the simulation input data. In addition, simulation input data can be affected by uncertainty in relative positioning of the anatomical model with respect to the RF coil. The objective of this work is the to estimate the variability of SAR and temperature increase at 3 T head MRI due to different sources of variability in input data, with the final aim to associate a global uncertainty to the dosimetric outcomes. Approach. A stochastic approach based on arbitrary Polynomial Chaos Expansion is used to evaluate the effects of several input variability’s (anatomy, tissue properties, body position) on dosimetric outputs, referring to head imaging with a 3 T MRI scanner. Main results. It is found that head anatomy is the prevailing source of variability for the considered dosimetric quantities, rather than the variability due to tissue properties and head positioning. From knowledge of the variability of the dosimetric quantities, an uncertainty can be attributed to the results obtained using a generic anatomical head model when SAR and temperature increase values are compared with safety exposure limits. Significance. This work associates a global uncertainty to SAR and temperature increase predictions, to be considered when comparing the numerically evaluated dosimetric quantities with reference exposure limits. The adopted methodology can be extended to other exposure scenarios for MRI safety purposes.