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 TeV, corresponding to an integrated luminosity of 35.9 $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 W or Z boson, or a top quark-antiquark pair) and the following decay modes: H $\rightarrow$ γγ, ZZ, WW, ττ, bb, and μμ. 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 ± 0.10$, assuming a Higgs boson mass of 125.09 GeV. Additional results are given for parametrizations with varying assumptions on the scaling behavior of the different production and decay modes, including generic ones based on ratios of cross sections and branching fractions or coupling modifiers. The results are compatible with the standard model predictions in all parametrizations considered. In addition, constraints are placed on various two Higgs doublet models.