| 4 | 1/1 | 返回列表 |
| 查看: 1066 | 回復(fù): 3 | ||||
| 【懸賞金幣】回答本帖問題,作者kk1424將贈送您 50 個(gè)金幣 | ||||
kk1424主管區(qū)長 (著名寫手)
中藥研究者
|
[求助]
求高手幫忙為小弟做一個(gè)四元回歸曲線,高分懸賞! 已有1人參與
|
|||
|
遇到棘手的問題,我有MATLAB軟件,但不是很會。求高手達(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 分別表示殘差及殘差對應(yīng)的置信區(qū)間。StatS 數(shù)組包含三個(gè)數(shù)字,分別是相關(guān)系數(shù),F(xiàn) 統(tǒng)計(jì)量及對應(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)%畫出在置信度區(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得到的一樣?梢姡竽娴膯栴}也是用了最小二乘的原理。 但是,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ā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +6 | zjl的號 2026-03-16 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林,給自己招師弟中~ +3 | TqlXswl 2026-03-16 | 6/300 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 285化工學(xué)碩求調(diào)劑(081700) +9 | 柴郡貓_ 2026-03-12 | 9/450 |
|
|
[考研] 268求調(diào)劑 +5 | 好運(yùn)連綿不絕 2026-03-12 | 5/250 |
|
|
[考研]
|
zhouzhen654 2026-03-16 | 3/150 |
|
|
[考研] 梁成偉老師課題組歡迎你的加入 +8 | 一鴨鴨喲 2026-03-14 | 9/450 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
[基金申請]
有必要更換申報(bào)口嗎
20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +9 | Copy267 2026-03-10 | 9/450 |
|
|
[考研] 337一志愿華南理工0805材料求調(diào)劑 +7 | mysdl 2026-03-11 | 9/450 |
|
|
[考研] 一志愿中科院,化學(xué)方向,295求調(diào)劑 +4 | 一氧二氮 2026-03-11 | 4/200 |
|
|
[考研] 一志愿211化學(xué)學(xué)碩310分求調(diào)劑 +8 | 努力奮斗112 2026-03-12 | 9/450 |
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|
|
[考研] 08食品或輕工求調(diào)劑,本科發(fā)表3篇sci一區(qū)top論文,一志愿南師大食品科學(xué)與工程 +3 | 我是一個(gè)兵, 2026-03-10 | 3/150 |
|
|
[考研] 270求調(diào)劑 085600材料與化工專碩 +3 | YXCT 2026-03-11 | 3/150 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|