數(shù)模regress函數(shù)在Matlab統(tǒng)計(jì)工具箱中使用命令實(shí)現(xiàn)多元線性

1、在Matlab統(tǒng)計(jì)工具箱中使用命令regress()實(shí)現(xiàn)多元線性回歸，調(diào)用格式為 或b，bint，r，rint，statsl=其中因變量數(shù)據(jù)向量yx出向量b，bint為回歸系數(shù)估計(jì)值和它們的置信區(qū)間，r，rint為殘差及其置信區(qū)間，stats 是用于檢驗(yàn)回歸模型的統(tǒng)計(jì)量，有三個(gè)數(shù)值，第一個(gè)是R2R時(shí)拒絕 H0，回歸模型成立畫(huà)出殘差及其置信區(qū)間，用命令1COD濃度實(shí)測(cè)值(y)輸入數(shù)x1=1.376,1.375,1.387,1.401,1.412,1.428,1.445,x2=0.450,0.475,0.485,0.500,0.535,0.545,0.550,x3=2.170,2.554,2.67

2、6,2.713,2.823,3.088,3.122,x4=0.8922,1.1610,0.5346,0.9589,1.0239,1.0499,1.1065,y=5.19,5.30,5.60，5.82，6.00,6.06，6.45，保存數(shù)據(jù)(以數(shù)據(jù)文件.mat形式保存，便于以后調(diào)用) save data x1 x2 x3 x4 yloaddata取出數(shù)據(jù)) x =ones(8,1)， ； b，bint，r，rint，statsregressb=(-16.5283，15.7206，2.0327，-0.2106，-0.1991) stats = (0.9908，80.9530，0.0022)即=-1

3、6.5283+15.7206xl+2.0327x2-0.2106x3+0.1991x4 R2 = 0.9908，F(xiàn) = 80.9530，P = 0.0022非線性回歸可由命令nlinfitbeta,r,j=其中，輸人數(shù)據(jù)x，y分別為nmn維列向量，對(duì)一元非線性回歸，x為n維列向量modelm-文件定義的非線性函數(shù)，beta0是回歸系數(shù)的初值， beta 是估計(jì)出的回歸系數(shù)，r 是殘差，jJacobian 矩陣，它們是估計(jì)y,delta= (1)對(duì)所要擬合的非線性模型建立的m-文件mode1.m如下： function yhat=model(beta,t) 輸人數(shù)loaddatay(在data.

4、maty) 即預(yù)測(cè)及作yy，delta=nlprodei(model，t，beta,r，j)； stepwise，它提供了一個(gè)交互式畫(huà)面，通過(guò)此工具可以自其中x是自變量數(shù)據(jù)，y是因變量數(shù)據(jù)，分別為nm和nl矩陣，inmodel是矩陣的列數(shù)指標(biāo)(缺省時(shí)為全部自變量)，alpha,為顯著性水平(0.5) stepwise Table 窗口中(RMSE)，相關(guān)系數(shù)(R-square)，F(xiàn)值和PMATLAB中進(jìn)行回歸分析bregress(y，x)y和回歸矩陣xb，bint，r,rint，statsregress(y,x)則給出系數(shù)的估計(jì)值b；系數(shù)估計(jì)值的置信度為95的置信區(qū)間bintr及各殘差的置信區(qū)

5、間rint;statsR2統(tǒng)計(jì)量和FP值b，bint，r,rint，statsregess(y，alpha)給出置信度為 1-alpha 的結(jié)果，他符號(hào)意義同上1結(jié)果說(shuō)明：b 為回歸模型中的常數(shù)項(xiàng)及回歸系數(shù). Bint 為各系數(shù)的 95%置信區(qū)間r和rint為對(duì)應(yīng)每個(gè)實(shí)際值的殘差和殘差置信區(qū)間。Stats向量的值分別為擬合優(yōu)度、F值和顯著性概率p. . b,bint,r,rint,ststs=regress(y,x,alpha) 可用help查閱此命令的具體用法殘差及置信區(qū)間可以b,bint,r,rint,stats=regress(y,X,0.05); b,bint,stats,rcoplot(r,rint) 為 0=27.0269 1=140.61940的置信區(qū)間是22.3226,31.73131的111.7842,169.4546建模軟件建模軟件matlab R2=0.9219F=118.0670,p10-4.R是衡量y與x的相關(guān)程度的指標(biāo),稱為相關(guān)系數(shù).R越大,x與y關(guān)系越密切.R0.9才認(rèn)為相關(guān)關(guān)系成立. F 是一統(tǒng)計(jì)指標(biāo),p