Code:
. forvalues i=1(27)2052 {
2. matrix num = J(1,1,0)
3. matrix C=A[`i'..(`i'+27-1),1...]
4. svmat C
5. matrix accum cov = C1-C24, deviations noconstant
6. matrix cov=cov/(r(N)-1)
7. matrix ST=ST+cov
8. drop C1-C24
9. forvalues k = 1/24 {
10. mat num = num + (e`k" * cov * e`k')
11. }
12. mat I_`t' = num/trace(cov)
13. mat list I_`t'
14. local t = `t'+1
15. }
(obs=27)
too few quotes
r(132);
I am on version 13

Comment