全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 MATLAB等数学软件专版
3556 2
2013-02-06
我想利用Matlab工具箱中的genetic algorithm 來求解hansen and seo(2002)' Testing for two-regime threshold cointegration in vector error correction models'  中的門檻值和cointegration vector

但genetic algorithm要設定 fitness function(文獻建議設上面這篇的likelihood function)
所以我的 fitness function 設定如下 但跑不出來  有人可以跟我說正確的程式寫法嗎??

function y = store(betab,gam)
global x;
global y;
global xlag;
load aniexc.txt;
n=length(dat(:,1));
y=dat(2+k:n,:)-dat(1+k:n-1,:);
t=length(y(:,1));
xlag=dat(1+k:n-1,:);
x=ones(t,1);
y = (t/2)*log(det((u-zzj*bz)'*(u-zzj*bz)/t))
xx=inv(x'*x);
u=y-x*xx*(x'*y);
w=xlag*[1;-betab];
z=[w,x];
d1=(w<=gam);
zj=[(z.*(d1*ones(1,length(z(1,:))))),w];
zzj=zj-x*xx*(x'*zj);
warning off;
lastwarn(' ');
bz=(u'/zzj')';
[mw,idw] = lastwarn;
lastwarn(' ');
warning on;
if (1-mw==' ');
bz=pinv(zzj'*zzj)*(zzj'*u);
end;

目前很急需高手來幫我解答這問題
我是論壇新手  沒有甚麼金幣送解答者  
我相信版上一定有熱心高手能幫我這個大忙  感謝
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2013-2-6 01:25:33
楼主一个币都没有啊? 呵呵。 还不赶紧上传点好资料赚点币,网上已经有的资料就不用浪费精力了。
有好资料就可赚币,有币可使鬼推磨。 知道不?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2013-2-6 19:46:45
tigerwolf 发表于 2013-2-6 01:25
楼主一个币都没有啊? 呵呵。 还不赶紧上传点好资料赚点币,网上已经有的资料就不用浪费精力了。
有好资料 ...
這就是我的問題了,  我不是不想給幣 是我剛進論壇不久 根本沒幣 且時間蠻急迫了 我自己試了很久 可是還是做不出來  請高手幫幫忙 幫我解答一下
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群