The paper presents a numerical approach to the analysis of the statistical effect of functional groups on the sorption of methanol and water on hard coal samples. The material used for the analysis was obtained from numerous samples of hard subbituminous and bituminous coals up to anthracite from different Polish coal mines and includes sorption isotherms of water and methanol vapors, as well as carbon dioxide and methane on these samples. Measurements were made of the sorption isotherms of water and methanol vapors and the data set has been supplemented with the sorption isotherms of water, methane and carbon dioxide taken from the literature for the precise estimation of the model parameters. More precise estimation is reached by using a strong setting of the coal structure parameters for a bigger number of sorption system, and in each case the same parameters of coal geometry are constant with an exact fitting of sorption isotherm. The adsorption-absorption model of the sorption in coal (Multiple Sorption Model-MSM) is used in the numerical experiments and the parameters of hard coal structure and the sorption systems are estimated. It has been stated that water as a polar substance together with methanol, as well as carbon dioxide and methane give good estimates of coal structure and let us quantify the polar effect of surface groups present in hard coal. The polar effect is introduced in the model in the range of adsorption and expansion subprocesses. The presence of oxygen groups in the bulk of coal matter has no significant effect and can be neglected. A weak decreasing tendency of polar factors ratio for water and methanol is discovered.