where I is an nixni identity matrix, VN is the “naïve” or model-based variance estimator and Hi=D^iVND^iTV^i−1. Mancl and DeRouen justify this correction on the grounds that the true expected value is expressed as E[Cov(Yi)]≈(I−Hi)(Yi−μi)(Yi−μi)T(I−HiT) rather than E[Cov(Yi)]=(Yi−μi)(Yi−μi)T.