| 4 | 1/1 | 返回列表 |
| 查看: 1065 | 回復(fù): 3 | ||||
| 【懸賞金幣】回答本帖問(wèn)題,作者kk1424將贈(zèng)送您 50 個(gè)金幣 | ||||
kk1424主管區(qū)長(zhǎng) (著名寫(xiě)手)
中藥研究者
|
[求助]
求高手幫忙為小弟做一個(gè)四元回歸曲線,高分懸賞! 已有1人參與
|
|||
|
遇到棘手的問(wèn)題,我有MATLAB軟件,但不是很會(huì)。求高手達(dá)人,幫忙給個(gè)編程。。。。感激不盡. 這回歸方程: y1=0.227220-0.002859x1+0.018825x2-0.003369x3-0.021391x4 拜托了! |

榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
榮譽(yù)版主 (文壇精英)
![]() |
專家經(jīng)驗(yàn): +518 |
|
Matlab多元線性回歸 [ b , bint , r , rint , stats ]=regress ( y , x ) , 其中b 是回歸方程中的參數(shù)估計(jì)值,bint 是b 的置信區(qū)間,r 和rint 分別表示殘差及殘差對(duì)應(yīng)的置信區(qū)間。StatS 數(shù)組包含三個(gè)數(shù)字,分別是相關(guān)系數(shù),F(xiàn) 統(tǒng)計(jì)量及對(duì)應(yīng)的概率p 值。 擬合結(jié)果: Y=b(1)x(1)+b(2)x(2)+b(3)x(3)+…+b(n)x(n) b(1)是系數(shù),x(1)為全1的一個(gè)列向量。 注意:不是插值。 x=[1097 1284 1502 1394 1303 1555 1917 2051 2111 2286 2311 2003 2435 2625 2948 3155 3372];%因變量時(shí)間序列數(shù)據(jù) y=[698 872 988 807 738 1025 1316 1539 1561 1765 1762 1960 1902 2013 2446 2736 2825];%自變量時(shí)間序列數(shù)據(jù) X=[ones(size(x')),x']; [b,bint,r,rint,stats]=regress(y',X,0.05);%調(diào)用一元回歸分析函數(shù) rcoplot(r,rint)%畫(huà)出在置信度區(qū)間下誤差分布. 舉例: x = 1 2 4 9 1 4 3 7 1 5 9 0 1 9 1 8 >> y=[10 3 90 48]'; >> [ b , bint , r , rint , stats ]=regress ( y , x ) 得到的結(jié)果 b = -186.8333 16.0238 21.8571 8.5952 bint = NaN NaN NaN NaN NaN NaN NaN NaN r = 1.0e-013 * -0.5684 -0.9948 -0.1421 0.1421 rint = NaN NaN NaN NaN NaN NaN NaN NaN stats = 1 NaN NaN NaN 另外,用b=inv(x)*y 得到的結(jié)果和上面用regression得到的一樣?梢(jiàn),求逆的問(wèn)題也是用了最小二乘的原理。 但是,regress更優(yōu)于inv,體現(xiàn)在當(dāng)x的元素存在缺陷時(shí)。舉例: x = 1 0 0 4 1 0 0 3 1 0 0 5 1 0 0 8 y = 3 4 3 9 [ b , bint , r , rint , stats ]=regress ( y , x ) Warning: X is rank deficient to within machine precision. b = -0.9643 0 0 1.1429 bint = -12.0010 10.0724 0 0 0 0 -0.9245 3.2102 r = -0.6071 1.5357 -1.7500 0.8214 rint = -9.2316 8.0173 -1.9141 4.9855 -7.4997 3.9997 0.2465 1.3964 stats = 0.7388 5.6575 0.1405 3.2321 |
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 289求調(diào)劑 +6 | 步川酷紫123 2026-03-11 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 268求調(diào)劑 +6 | 一定有學(xué)上- 2026-03-14 | 7/350 |
|
|
[考研] 0703化學(xué)336分求調(diào)劑 +3 | zbzihdhd 2026-03-15 | 4/200 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[考研] 311求調(diào)劑 +5 | 26研0 2026-03-15 | 5/250 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 一志愿華中師范071000,325求調(diào)劑 +6 | RuitingC 2026-03-12 | 6/300 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢(mèng)微 2026-03-15 | 3/150 |
|
|
[考研] 344求調(diào)劑 +3 | knight344 2026-03-16 | 3/150 |
|
|
[基金申請(qǐng)] 現(xiàn)在如何回避去年的某一個(gè)專家,不知道名字 +3 | zk200107 2026-03-12 | 6/300 |
|
|
[考研] 328求調(diào)劑 +3 | 5201314Lsy! 2026-03-13 | 6/300 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 308 085701 四六級(jí)已過(guò)求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 330求調(diào)劑 +3 | ?醬給調(diào)劑跪了 2026-03-13 | 3/150 |
|
|
[考研] 材料工程,326分,求調(diào)劑 +6 | KRSLSR 2026-03-10 | 6/300 |
|
|
[考研] 341求調(diào)劑 +4 | 番茄頭--- 2026-03-10 | 4/200 |
|
|
[考研] 315求調(diào)劑 +9 | 小羊小羊_ 2026-03-11 | 10/500 |
|
|
[考研] 281求調(diào)劑 +9 | Koxui 2026-03-12 | 11/550 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
[考研] 420求調(diào)劑 +4 | 莫向外求11 2026-03-10 | 6/300 |
|