Abstract Combined measurements of the production and decay rates of the Higgs boson, as well as its couplings to vector bosons and fermions, are presented. The analysis uses the LHC proton–proton collision data set recorded with the CMS detector in 2016 at $$\sqrt{s}=13\,\text {Te}\text {V} $$ s=13Te , corresponding to an integrated luminosity of 35.9$${\,\text {fb}^{-1}} $$ fb-1 . The combination is based on analyses targeting the five main Higgs boson production mechanisms (gluon fusion, vector boson fusion, and associated production with a $$\mathrm {W}$$ W or $$\mathrm {Z}$$ Z boson, or a top quark-antiquark pair) and the following decay modes: $$\mathrm {H} \rightarrow \gamma \gamma $$ H→γγ , $$\mathrm {Z}\mathrm {Z}$$ ZZ , $$\mathrm {W}\mathrm {W}$$ WW , $$\mathrm {\tau }\mathrm {\tau }$$ ττ , $$\mathrm {b} \mathrm {b} $$ bb , and $$\mathrm {\mu }\mathrm {\mu }$$ μμ . Searches for invisible Higgs boson decays are also considered. The best-fit ratio of the signal yield to the standard model expectation is measured to be $$\mu =1.17\pm 0.10$$ μ=1.17±0.10 , assuming a Higgs boson mass of $$125.09\,\text {Ge}\text {V} $$ 125.09Ge . Additional results are given for various assumptions on the scaling behavior of the production and decay modes, including generic parametrizations based on ratios of cross sections and branching fractions or couplings. The results are compatible with the standard model predictions in all parametrizations considered. In addition, constraints are placed on various two Higgs doublet models.