Efficient B3LYP and BHandH density functionals were used to estimate methanol's nuclear magnetic isotropic shieldings and spin-spin coupling constants in the basis set limit. Polarization-consistent pcS-n and pcJ-n (n = 0, 1, 2, 3 and 4), and segmented contracted XZP, where X = D, T, Q and 5, basis sets were used and the results fitted with simple mathematical formulas. The performance of the methods was assessed from comparison with experiment and higher level calculations. 1J(CH) and 3J(HH) values were determined from very diluted solutions in deuterochloroform and compared with theoretical predictions. The agreement between complete basis set (CBS) density functional theory (DFT) predicted isotropic shieldings and spin-spin values and experiment was good. The BHandH/pcS-n methanol shieldings obtained using structures optimized at the same level of theory are approaching the accuracy of the advanced coupled-cluster-singles-doubles-approximate triples (CCSD(T)) calculations.