The molecular motion in heterogeneous media displays anomalous diffusion by the mean-squared displacement $\langle X^2(t) \rangle = 2 D t^α$. Motivated by experiments reporting populations of the anomalous diffusion parameters $α$ and $D$, we aim to disentangle their respective contributions to the observed variability when this last is due to a true population of these parameters and when it arises due to finite-duration recordings. We introduce estimators of the anomalous diffusion parameters on the basis of the time-averaged mean squared displacement and study their statistical properties. By using a copula approach, we derive a formula for the joint density function of their estimations conditioned on their actual values. The methodology introduced is indeed universal, it is valid for any Gaussian process and can be applied to any quadratic time-averaged statistics. We also explain the experimentally reported relation $D\propto\exp(αc_1+c_2)$ for which we provide the exact expression. We finally compare our findings to numerical simulations of the fractional Brownian motion and quantify their accuracy by using the Hellinger distance.