%k=0:1e-3:10;
%bb=ridge(y,x,k);
%figure;
%plot(k,bb);
%xlabel('k');
%ylabel('β');
%title('岭迹');
%legend('x1','x2','x3','x4','x5','x6');
通过岭迹图来确定k值。(1)plot(lmridge(GDP~Consume+Investment+IO+Population+Jobless+Goods,data=dat,lambda=seq(0,03,0001)))#和线性回归类似,这个plot可以画出岭迹图,lambda=seq(0,03,0001)设置范围和间隔,可以观察岭迹图,人工选择,但是这样主观性较强。(2)select(lmridge(GDP~Consume+Investment+IO+Population+Jobless+Goods,data=dat,lambda=seq(0,03,0001)))#利用select函数找出最优岭参数lambda,会有三个值,任选一个即可。lmridge(GDP~Consume+Investment+IO+Population+Jobless+Goods,data=dat,lambda=009)#通过(1)或(2)把选取的lmbda参数写到岭回归函数中去,在这里lambda=009。
欢迎分享,转载请注明来源:内存溢出
评论列表(0条)