多元線性回歸

• 一元線性回歸
• • 繪圖
  • 求參
• 多元線性回歸
• • regress函數(shù)
• matlab 文檔非常有用

一元線性回歸

一個(gè)舊東西隨著時(shí)間價(jià)格也在變動(dòng)，數(shù)據(jù)如下：

12345678910

265	1942	1493	1086	766	539	485	291	224	202

繪圖

clear all clc x=1:10; y=[2650,1942,1493,1086,766,539,485,291,224,202]; for i=1:10plot(x(i),y(i),'or');hold on end xlabel('x'); ylabel('y');

由圖可發(fā)現(xiàn)，x和y呈現(xiàn)指數(shù)關(guān)系，于是我們假設(shè)$z=\ln y$,記作$z_i=\ln y_i$,重新繪圖：

clear all clc x=1:10; y=[2650,1942,1493,1086,766,539,485,291,224,202]; z=zeros(size(y)); N=length(y); hold on for i =1:Nz(i)=log(y(i));plot(x(i),z(i),'ok'); end xlabel('x'); ylabel('y');

分析發(fā)現(xiàn)現(xiàn)在變得更擬合了。各點(diǎn)基本處于一條直線上，可以認(rèn)為$z=a+bx$

求參

clear all clc x=1:10; y=[2650,1942,1493,1086,766,539,485,291,224,202]; z=zeros(size(y)); N=length(y); for i =1:Nz(i)=log(y(i)); end [p,s]=polyfit(x,z,1)

返回結(jié)果：

可以得出：a=8.1671,b=-0.2984,即$z=8.1671?0.2984x$

多元線性回歸

regress函數(shù)

%調(diào)用格式 [b,bint,r,rint,stats]=regress(y,x，alpha)

alpha為顯著性水平，缺省設(shè)定為0.05, b表示為輸出量, bint為回歸系數(shù)估計(jì)值和他們的置信區(qū)間，r為殘差, rint為置信區(qū)間, stats適用于檢驗(yàn)回歸模型的統(tǒng)計(jì)量。
案例如下:

$x_1$1.3761.3751.3871.4011.4121.4281.4451.477

$x_2$	0.450	0.475	0.485	0.5	0.535	0.545	0.55	0.575
$x_3$	2.17	2.554	2.676	2.713	2.823	3.088	3.122	3.262
$x_4$	0.8922	1.161	0.5346	0.9589	1.0239	1.0499	1.1065	1.1387
$y$	5.19	5.3	5.6	5.82	6	6.06	6.45	6.95

%開(kāi)始多元回歸 clear all clc x1=[1.376,1.375,1.387,1.401,1.412,1.428,1.445,1.477]; x2=[0.450,0.475,0.485,0.5,0.535,0.545,0.55,0.575]; x3=[2.170,2.554,2.676,2.713,2.823,3.088,3.122,3.262]; x4=[5.19,1.161,0.5346,0.9589,2.0239,1.0499,1.1065,1.1387]; y=[5.19,5.3,5.6,5.82,6,6.06,6.45,6.95]; save data x1 x2 x3 x4 y%保存數(shù)據(jù) load data %取出數(shù)據(jù) y=[y'];%轉(zhuǎn)置矩陣 x=[ones(size(x1')),x1',x2',x3',x4']; [b,bint,r,rint,stats]=regress(y,x)

根據(jù)b可以解得該函數(shù)關(guān)系式為：
$$y=?20.5297+19.1269x_1+8.0045x_2?1.5867x_3?0.1465x_4$$

matlab 文檔非常有用

文檔連接：點(diǎn)擊這里

總結(jié)

以上是生活随笔為你收集整理的【Matlab】多元线性回归的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。