作者henry7448 (P=MC)
看板MATLAB
标题[问题]急~关於genetic algorithm中fitness 写法
时间Wed Feb 6 09:02:35 2013
我想利用Matlab optimization工具箱中的(GA)genetic algorithm 来求解hansen
and seo(2002)' Testing for two-regime threshold cointegration in vector
error correction models' 中的门槛值和共积值
但GA要设定 fitness function(文献建议设成上面这篇的likelihood function)
所以我的 fitness function 设定如下 但跑不出来
我想扩充至4变数2门槛,所以出来结果应该要有2门槛值3共积值
拜托高手帮忙解答,若有高手解答出来,我愿意赠送至少100枚以上P币 谢谢
function r = 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);
r = (t/2)*log(det((u-zzj*bz)'*(u-zzj*bz)/t)) @这就是论文的likelihood
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;
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 114.38.114.64
1F：→ henry7448:我刚看了一下 若有高手能帮忙解答出来 我愿意赠送1000P 02/06 09:05