This function treats measure-specific variance as reliable.
The Mosier composite formula is computed as:
$$\rho_{XX}=\frac{\mathbf{w}^{T}\left(\mathbf{r}\circ\mathbf{s}\right)+\mathbf{w}^{T}\mathbf{S}\mathbf{w}-\mathbf{w}^{T}\mathbf{s}}{\mathbf{w}^{T}\mathbf{S}\mathbf{w}}$$
where \(\rho_{XX}\) is a composite reliability estimate, \(\mathbf{r}\) is a vector of reliability estimates, \(\mathbf{w}\) is a vector of weights, \(\mathbf{S}\) is a covariance matrix, and \(\mathbf{s}\) is a vector of variances (i.e., the diagonal elements of \(\mathbf{S}\)).