先贴一下GRS论文中的stata源程序, 这个GRS程序是网上公开的,网址是ideas的working paper 附带程序
https://ideas.repec.org/c/boc/bocode/s457786.html
先说说程序中的几个错误
第一
matrix R[2,5] = `meanSE'
matrix R[1,6] = `meanabsalpha'
matrix R[2,6] = `meanabsalpha'
qui correlate `d', covariance
matrix F=r(C)
*--- Wald test (J0) ----*
matrix T=(1+ M' * invsym(F) * M)
matrix H=1/T[1,1]
matrix R[1,2] = e(N)* H * A' * invsym(S
这个
qui correlate `d', covariance
matrix F=r(C)
是错误的,不能用因子平均收益的相关系数矩阵,应该改为
(F-F的平均数)‘ *(F-F的平均数)
这是第一个错误
下面这段也不对
qui mat accum RR = Res*, dev noconstant
mat S =RR/(e(N)-1)
drop Res*
mat drop RR
应该改为
mat S =RR/(e(N)-'b'-1)