Я пытаюсь воспроизвести отдельные особенности c и отдельные отклонения из пакета plm при использовании случайных эффектов. Мне известно, что в пакете есть разные способы вычисления этих отклонений, а также того, какие степени свободы использовать, но независимо от используемых методов и DF, мои результаты никогда не совпадают с результатами из пакета plm.
Предполагая, что у нас есть 4 человека с 20 наблюдениями, для индивидуальной дисперсии c я вычисляю
(( sum( resid(lm(y1 ~ x1_1 + x2_1, data=data [id==1,])) ^2 ) ) +
( sum( resid(lm(y1 ~ x1_1 + x2_1, data=data [id==2,])) ^2 ) ) +
( sum( resid(lm(y1 ~ x1_1 + x2_1, data=data [id==3,])) ^2 ) ) +
( sum( resid(lm(y1 ~ x1_1 + x2_1, data=data [id==4,])) ^2 ) ) ) / ((20-2-1)*4).