The Veytsman statistics for hydrogen bonding (HB) contribution is combined with the nonrandom lattice fluid (NLF) model developed recently by the present authors. The physical contribution is characterized by two temperature dependent molecular parameters representing molecular size and interaction energy for a pure fluid, and a binary interaction parameter for a binary mixture. The chemical part is represented by internal energy and entropy of HE. The resulting NLF-HB equation of state was applied to describe thermodynamic properties of pure fluids and phase equilibrium of mixtures. The results for alkane-alcohol mixtures showed significant improvements over those of the NLF EOS (equation of state) and HM (hole quasichemical group-contribution) model.