〖前人指路〗非线性、逐步回归,MATLAB实现-经管之家官网！

1. stepwise(x,y)

——MATLAB技术论坛

Fantuanxiaot 坛友 补充的非线性回归中的稳健回归！！

函数

1. function S=Nonlinearmodel(beta,x)
   
2. S=beta(1)*exp(beta(2)*x);

主程序

1. x=[10:0.5:20]';
   
2. y=2*exp(0.5*x)+normrnd(0,1,21,1);
   
3. beta0=[1 1]';
   
4. options=statset('MaxIter',200,'Robust','on','WgtFun','huber');
   
5. [beta,resid,jacobimatrix,covbeta,MSE]=nlinfit(x,y,'Nonlinearmodel',beta0,options);%非线性稳健回归
   
6. [yy,delta]=nlpredci('Nonlinearmodel',x,beta,resid,'jacobian',jacobimatrix);
   
7. ci=nlparci(beta,resid,jacobimatrix);
   
8. nlintool(x,y,'Nonlinearmodel',beta0,0.05);

基于内联函数的例子

1. [free]x=poissrnd(2,30,1);
   
2. y=poissrnd(3,30,1);
   
3. X=[x y];
   
4. Z=2*exp(x+y)+rand(30,1);
   
5. beta0=[1.6 1.2 1.4];
   
6. EXPFUN=inline('b(1)*exp(b(2)*G(:,1)+b(3)*G(:,2))','b','G');
   
7. options=statset('MaxIter',200,'Robust','on','WgtFun','huber');
   
8. [beta,resid,jacobimatrix,covbeta,MSE]=nlinfit(X,Z,EXPFUN,beta0,options);
   
9. ci=nlparci(beta,resid,jacobimatrix);
   
10. [yy,delta]=nlpredci(EXPFUN,X,beta,resid,'jacobian',jacobimatrix);[/free]

「经管之家」APP：经管人学习、答疑、交友，就上经管之家！
免流量费下载资料----在经管之家app可以下载论坛上的所有资源，并且不额外收取下载高峰期的论坛币。
涵盖所有经管领域的优秀内容----覆盖经济、管理、金融投资、计量统计、数据分析、国贸、财会等专业的学习宝库，各类资料应有尽有。
来自五湖四海的经管达人----已经有上千万的经管人来到这里，你可以找到任何学科方向、有共同话题的朋友。
经管之家（原人大经济论坛），跨越高校的围墙，带你走进经管知识的新世界。
扫描下方二维码下载并注册APP