因子模型代码.pdf

  1. 1、本文档共18页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
Data • Time-series Regressions • CAPM R −R  + (R −R ) + it f i i mt f it • FF3F Model R −R  + (R −R ) +s smb +h hml + it f i i mt f i t i t it Part I CAPM % loading data x = load('25_Portfolios_5x5.txt'); T1 = size(x,1); labels = { 's1v5' 's1v4' 's1v3' 's1v2' 's1v1' ... 's2v5' 's2v4' 's2v3' 's2v2' 's2v1' ... 's3v5' 's3v4' 's3v3' 's3v2' 's3v1' ... 's4v5' 's4v4' 's4v3' 's4v2' 's4v1' ... 's5v5' 's5v4' 's5v3' 's5v2' 's5v1'}; x = load('F-F_Research_Data_Factors.txt'); x = x(st:st-1+size(r,1),:); rmrf = x(1:en2,2); smb = x(1:en2,3); hml = x(1:en2,4); rf = x(1:en2,5); dates = x(1:en2,1); dates = floor(dates/100)+(dates-100*floor(dates/100))*1/12; % turns 192603 in to 1926+3/12; Part I CAPM % loading data T = length(rmrf); N = size(r,2); rx = r-rf*ones(1,N); %CAPM f = [rmrf]; [alpha, beta, R2, siga, sigb] = tsregress(rx,f); Part I CAPM • tsregress.m function: function [alpha, beta, R2, siga, sigb] = tsregress(rx,f); T = size(rx,1); N = size(rx,2); K = size(f,2); if size(f,1) ~= size(rx,1); disp('tsregress error: factors and returns must be of same length'); end; X = [ones(T,1) f]; b = X\rx; u = rx-X*b; s2 = var(u)'*(T-1)/(T-K-1); % produces a vector, variance of the 25 errors Part I CAPM • tsregress.m function( ): mx = inv(X'*X); dmx = diag(mx); % we're interested in standard errors, the diagonals of the covariance matrix of bs siga = (s2*dmx(1)).^0.5; % standard error

文档评论(0)

183****7931 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档

相关课程推荐