We perform a global sensitivity analysis of the binding energy and the charge radius of the nucleus ^{16}O to identify the most influential low-energy constants in the next-to-next-to-leading order chiral Hamiltonian with two- and three-nucleon forces. For this purpose, we develop a subspace-projected coupled-cluster method using eigenvector continuation Frame D. et al., Phys. Rev. Lett. 121, 032501 (2018)PRLTAO0031-900710.1103/PhysRevLett.121.032501. With this method, we compute the binding energy and charge radius of ^{16}O at more than 10^{6} different values of the 16 low-energy constants in one hour on a standard laptop computer. For relatively small subspace projections, the root-mean-square error is about 1% compared to full-space coupled-cluster results. We find that 58(1)% of the variance in energy can be apportioned to a single contact term in the ^{3}S_{1} wave, whereas the radius depends sensitively on several low-energy constants and their higher-order correlations. The results identify the most important parameters for describing nuclear saturation and help prioritize efforts for uncertainty reduction of theoretical predictions. The achieved acceleration opens up an array of computational statistics analyses of the underlying description of the strong nuclear interaction in nuclei across the Segrè chart.