ここで与えられたFama-MacBeth 共同検定統計への入力として、価格設定エラーの共分散行列を計算しようとしています。
これを正しく行っているかどうかはわかりませんが、助けを求めていました。
これらの数式の写真を投稿するつもりでしたが、画像を投稿するには少なくとも 10 の評判ポイントが必要です...だから imageur を使用します
ここに私が計算しようとしている共同テストがあります:
共分散式は次のとおりです。
私のコードは次のとおりです
alpha.diff1 <- full.resid-full.mean.resid
alpha.diff2 <- t(alpha.diff1)
alpha.sum <- alpha.diff1%*%alpha.diff2
cov.term <- alpha.sum/(full.num.t^2)
full.joint.test <- t(full.mean.resid)%*%ginv(cov.term)%*%full.mean.resid
変数名はかなり自明だと思います。
- full.num.t : 数式の T (月ごとの総観測数)
- full.resid : alpha(it) [25 行の資産、192 列の月次誤差] を格納します
full.joint.test変数の結果は20.86672 です
私がこれを正しく行ったかどうかを誰かが確認できますか?