Quantification of fat has been investigated using images acquired from multiple gradient echoes. The evolution of the signal with echo time and flip angle was measured in phantoms of known fat and water composition and in 21 research subjects with fatty liver. Data were compared to different models of the signal equation, in which each model makes different assumptions about the T1 and/or T2* relaxation effects. A range of T1, T2*, fat fraction and number of echoes was investigated to cover situations of relevance to clinical imaging. Results indicate that quantification is most accurate at low flip angles (to minimize T1 effects) with a small number of echoes (to minimize spectral broadening effects). At short echo times, the spectral broadening effects manifest as a short apparent T2 for the fat component.