Compute covariance of fixed effect parameters as a function of variance parameters of a linear mixed model