| 4 | 1/1 | 返回列表 |
| 查看: 1244 | 回復(fù): 3 | ||||
dota_sao木蟲 (正式寫手)
|
[求助]
用1stopt非線性回歸,擬合度較差 已有2人參與
|
|
代碼如下,擬合度較差,且不唯一,不知該如何解決。 Variable m1,m2,la,m3,lr; Function lr= 0.125*( - c - 8*b1*m1 - 8*b2*m2 + 2*b2*m1 + 2*b3*m2 - 4*Ln(m1) + Ln(m2) - 8*Ln(m3) - 8*la ) ; Data; 0.019166667 0.068333333 -0.006515836 0.1913 -0.263498937 0.03 0.065625 -0.007161192 0.21 -0.267435897 0.006944444 0.071388889 -0.016739675 0.4777 -0.276948724 0.05590863 0.150554092 -9.48905E-05 0.002681523 -0.055947827 0.071875 0.1465625 -0.000821925 0.0238 -0.142264245 0.045625 0.153125 -0.010000846 0.2913 -0.277623889 0.0321875 0.156484375 -0.013227869 0.3817 -0.280268169 0.04578125 0.153085938 -0.020932431 0.5892 -0.267907337 0.440977131 0.171809289 -0.000114726 0.003247193 -0.061066124 0.46 0.167053571 -0.000855986 0.0248 -0.144514772 0.089464286 0.2596875 -0.003454672 0.1016 -0.22953092 0.102857143 0.256339286 -0.007095494 0.2081 -0.267069379 0.047321429 0.270223214 -0.01301493 0.3758 -0.280294355 0.953333333 0.200416667 -0.000493757 0.0142 -0.11595175 0.409791667 0.336302083 -0.001290656 0.0376 -0.168382974 0.256041667 0.374739583 -0.002142865 0.0628 -0.199826385 0.08125 0.4184375 -0.006048192 0.1777 -0.26013016 0.038333333 0.429166667 -0.045993692 1.1953 -0.170924146 0.199922496 0.608144376 -0.000282509 0.008068999 -0.091455914 1.677 0.238875 -0.000582886 0.0168 -0.124114221 1.044 0.397125 -0.001639299 0.0479 -0.183036058 0.815 0.454375 -0.001889458 0.0553 -0.191908338 0.2715 0.59025 -0.040700232 1.0749 -0.194270321 |
鐵桿木蟲 (職業(yè)作家)
|
已經(jīng)都解釋過了,過擬合參數(shù)不唯一,自己公式的問題。 ******************************* 公式中“8*b1*m1+8*b2*m2-2*b2*m1-2*b3*m2” 可以寫成:(8*b1-2*b2)*m1+(8*b2-2*b3)*m2 等價于:b1*m1+b2*m2 兩個參數(shù)均可,要用三個參數(shù),只能是完全多組解。 |
主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +1059 |
|
dingd大神正解。 結(jié)果不穩(wěn)定問題,把原公式中 - 8*b1*m1 - 8*b2*m2 + 2*b2*m1 + 2*b3*m2 ,替換成 b1*m1+b2*m2,減少參數(shù)個數(shù)即可得到唯一、穩(wěn)定的結(jié)果。 擬合度不好問題,說明原公式和數(shù)據(jù)不匹配,選用其他其他公式試試。 |

主管區(qū)長 (職業(yè)作家)
![]() |
專家經(jīng)驗: +1059 |
|
如果公式換成:lr= 0.125*( - c+ b1*m1+b2*m2+ d*Ln(m1) + e*Ln(m2) - f*Ln(m3) - g*la ) 可得到?jīng)Q定系數(shù)0.98+的結(jié)果: 相關(guān)系數(shù)(R): 0.992093956081203 相關(guān)系數(shù)之平方(R^2): 0.984250417692851 決定系數(shù)(DC): 0.984250417692851 卡方系數(shù)(Chi-Square): -0.00807575711571369 F統(tǒng)計(F-Statistic): 1312.36869451249 參數(shù) 最佳估算 ---------- ------------- c 3.09591250724687 b1 0.0156358088159255 b2 0.306871885892378 d 0.0527773560656148 e -0.0585877018241693 f 0.458262804355688 g 36.569403287351 |

| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0703化學(xué)336分求調(diào)劑 +5 | zbzihdhd 2026-03-15 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 材料專碩306英一數(shù)二 +8 | z1z2z3879 2026-03-16 | 10/500 |
|
|
[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3 | 脫穎而出 2026-03-16 | 6/300 |
|
|
[考研] 070300化學(xué)319求調(diào)劑 +4 | 錦鯉0909 2026-03-17 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑 +5 | 204818@lcx 2026-03-17 | 5/250 |
|
|
[考研] 工科材料085601 279求調(diào)劑 +3 | 困于星晨 2026-03-17 | 3/150 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 化學(xué)調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +7 | Ncdx123456 2026-03-13 | 8/400 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[教師之家] 焦慮 +7 | 水冰月月野兔 2026-03-13 | 9/450 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
[考研] 工科調(diào)劑 +4 | Jiang191123! 2026-03-11 | 4/200 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|
|
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|