下載本文檔
版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請進(jìn)行舉報或認(rèn)領(lǐng)
文檔簡介
1、多元 garch 模型預(yù)測的matlab 程序function parameters,loglikelihood, ht, likelihoods, stdresid,stderrors,a, b, scores= full_bekk_mvgarch(data,p,q,bekkoptions);% purpose:%to estimatea full bekk multivariate garch model.% usage:%parameters,loglikelihood, ht, likelihoods, stdresid,stderrors,a, b, scores= full_bek
2、k_mvgarch(data,p,q,options);% inputs:%data- a t by k matrix of zeromeanresiduals%p- the lag lengthof the innovation process%q- the lag lengthof the ar process%options- (optional) options for the optimization(fminunc)% outputs:%parameters- a (k*(k+1)/2+p*k2+q*k2vector of estimatedparameteters.f%or an
3、yk2 setof innovation or ar parametersx,%reshape(x,k,k)will give thecorrect matrix%to recoverc, useivech(parmaeters(1:(k*(k+1)/2)%loglikelihood - the loglikelihood of thefunction atthe optimum%ht- a k x k x t 3 dimension matrix of conditional covariances%likelihoods- a t by 1 vector of individual lik
4、elihoods%stdresid- a t by k matrix of multivariate standardizedresiduals%stderrors- a numparams2 square matrix of robust standaderrors(a(-1)*b*a(-1)*t(-1)%a- theestimatedinverseof thenon-robuststandarderrors%b- the estimatedcovarianceof tehscores%scores- a t by numparamsmatrix of individual scores%
5、needto try andgetsomesmartstartgin valuesif size(data,2) size(data,1)data=data;endt k=size(data);k2=k*(k+1)/2;scalaropt=optimset(fminunc);scalaropt=optimset(scalaropt,tolfun,1e-1,display,iter,diagnostics,on,diffmaxchange,1e-2);startingparameters=scalar_bekk_mvgarch(data,p,q,scalaropt);cchol=starting
6、parameters(1:(k*(k+1)/2);c=ivech(startingparameters(1:(k*(k+1)/2)*ivech(startingparameters(1:(k*(k+1)/2);newa=;newb=;for i=1:pnewa=newa diag(ones(k,1)*startingparameters(k*(k+1)/2)+i);endfor i=1:qnewb=newb diag(ones(k,1)*startingparameters(k*(k+1)/2)+i+p);endnewa=reshape(newa,k*k*p,1);newb=reshape(n
7、ewb,k*k*q,1);startingparameters=cchol;newa;newb;if nargin=6a=hessian_2sided(full_bekk_mvgarch_likelihood,parameters,data,p,q,k,k2,t);h=max(abs(parameters/2),1e-2)*eps(1/3);hplus=parameters+h;hminus=parameters-h;likelihoodsplus=zeros(t,length(parameters);likelihoodsminus=zeros(t,length(parameters);fo
8、r i=1:length(parameters)hparameters=parameters;hparameters(i)=hplus(i);holder,indivlike=full_bekk_mvgarch_likelihood(hparameters,data,p,q,k,k2,t);likelihoodsplus(:,i)=indivlike;endfor i=1:length(parameters)hparameters=parameters;hparameters(i)=hminus(i);holder,indivlike=full_bekk_mvgarch_likelihood(hparameters,data,p,q,k,k2,t);likelihoodsminus(:,i
溫馨提示
- 1. 本站所有資源如無特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁內(nèi)容里面會有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫網(wǎng)僅提供信息存儲空間,僅對用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時也不承擔(dān)用戶因使用這些下載資源對自己和他人造成任何形式的傷害或損失。
最新文檔
- 2024下半年工作計劃
- 2024年3月上半年四川省^會辦公廳直屬事業(yè)單位招考聘用工作人員筆試歷年典型考點解題思路附帶答案詳解
- 2024屆普通高等學(xué)校招生全國統(tǒng)一考試青桐鳴大聯(lián)考理科綜合化學(xué)試題(含答案解析)
- 文學(xué)作品情節(jié)分析模擬試卷
- 2.4洋務(wù)運動和邊疆危機(jī)課件 2024-2025學(xué)年統(tǒng)編版八年級歷史上冊
- 北師大版歷史七年教案大全
- 2024年《民法典》學(xué)習(xí)題庫及答案(共668題)
- 軟件開發(fā)技術(shù)基礎(chǔ) 第4版 課件 6-2 圖像處理編程
- 湖北省黃石市(2024年-2025年小學(xué)四年級語文)人教版期末考試(下學(xué)期)試卷及答案
- 高一化學(xué)《電解質(zhì)的電離》分層練習(xí)含答案解析
- 廣東學(xué)校衛(wèi)生制度匯編
- 運輸公司車輛掛靠協(xié)議
- 足療店承包合同
- 盲板的厚度計算與使用
- 名著片段賞析《匆匆》內(nèi)容簡介教育培訓(xùn)PPT模板
- 培訓(xùn)效果滿意度調(diào)查表.doc
- 遷墳協(xié)議書范本(精選5篇)
- 消防檢測報價表
- 人教部編版四年級數(shù)學(xué)上冊《角的度量單元復(fù)習(xí)》PPT教學(xué)課件
- 產(chǎn)品合格證模板
- 電纜電線配電工程量自動計算表Excel
評論
0/150
提交評論