computes the sample covariance matrix of the items,
uses the classical formula
$$\alpha = \frac{k}{k-1} \left(1 - \frac{\sum \sigma_i^2}{\sigma_X^2}\right),$$
where \(k\) is the number of items, \(\sigma_i^2\) are item variances,
and \(\sigma_X^2\) is the variance of the total score,
computes SEM as \(\text{SD}(X) \sqrt{1 - \alpha}\).