版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

wenzimu

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學(xué)

[求助] matlab中如何用regress實現(xiàn)批量回歸并出圖

現(xiàn)在有三組143*4的矩陣數(shù)據(jù)，分別命名為GDP，LGDP，GRO。現(xiàn)在想通過regress以GDP為因變量，LGDP,GRO為自變量進行回歸。如果是整體回歸，那直接調(diào)用函數(shù)regress就可以了，但現(xiàn)在遇到的問題是想每列數(shù)據(jù)分區(qū)間求各自的回歸方程并畫圖。比如第一列數(shù)據(jù)中，GDP(1:11,1)與LGDP(1:11,1)，GRO(1:11,1)；GDP(12:22,1)與LGDP(12:22,1)，GRO(12:22,1)；GDP(23:33,1)與LGDP(23:33,1)，GRO(23:33,1)......GDP(133:143,1)與LGDP(133:143,1)，GRO(133:143,1)這樣分別進行批量回歸，并畫出實際值GDP與擬合值GDP的圖；然后第二列、第三列、第四列數(shù)據(jù)也同樣進行這樣的擬合（如下圖所示），不知這樣該如何實現(xiàn)，還請高人幫忙~!

示意圖

回復(fù)此樓

» 猜你喜歡

281求調(diào)劑（0805）已經(jīng)有8人回復(fù)
環(huán)境領(lǐng)域全國重點實驗室招收博士1-2名已經(jīng)有3人回復(fù)
材料專碩306英一數(shù)二已經(jīng)有10人回復(fù)
301求調(diào)劑已經(jīng)有6人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有7人回復(fù)
302求調(diào)劑已經(jīng)有6人回復(fù)
26博士申請已經(jīng)有3人回復(fù)
268求調(diào)劑已經(jīng)有3人回復(fù)
311求調(diào)劑已經(jīng)有10人回復(fù)
被我言中：新模板不強調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

關(guān)于matlab中循環(huán)讀取多個excel文件的問題？已經(jīng)有3人回復(fù)
如何用SAS 進行回歸系數(shù)的顯著性檢驗。。。。。已經(jīng)有3人回復(fù)
DFT Scan之后使用MP2計算能量，如何建立批量輸入文件？已經(jīng)有6人回復(fù)
matlab中如何將雙坐標(biāo)系圖中的另外一個坐標(biāo)系設(shè)置為當(dāng)前坐標(biāo)系呢？已經(jīng)有10人回復(fù)
求助：非線性方程系數(shù)回歸問題？已經(jīng)有10人回復(fù)
非線性回歸參數(shù)..急求幫助！！謝謝已經(jīng)有8人回復(fù)
求助：哪位大俠懂多重logistic回歸的指導(dǎo)我分析一個數(shù)據(jù)（金幣要多少都行）已經(jīng)有6人回復(fù)
請問如何從matlab中導(dǎo)出這張圖數(shù)據(jù) 已經(jīng)有4人回復(fù)

1樓 2012-08-15 20:45:22

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wenzimu

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學(xué)

我自己用隨機數(shù)做了一下，但最后發(fā)現(xiàn)回歸參數(shù)不能一起出現(xiàn)在同一個矩陣中，最后矩陣中只有某一組的回歸參數(shù)，請問這種情況怎么解決��？
>> X=rand(143,4);
>> Y=rand(143,4);
>> C=ones(11,1);
>> for n=0:12
i=11*n
for j=1:4
Xj=X(1+i:11+i,j);
Xj=[Xj,C]
Yj=Y(1+i:11+i,j)
[b,bint,r,rint,stats]=regress(Yj,Xj)
end
end

贊一下

回復(fù)此樓

2樓2012-08-16 16:17:46

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

感謝參與，應(yīng)助指數(shù) +1

改成這樣

CODE:

[b(:,:,n+1),bint(:,:,n+1),r(:,:,n+1),rint(:,:,n+1),stats(:,:,n+1)]=regress(Yj,Xj);

贊一下

回復(fù)此樓

showmethemoney

3樓2012-08-16 18:05:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wenzimu

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學(xué)

引用回帖:

3樓: Originally posted by csgt0 at 2012-08-16 18:05:20
改成這樣

=regress(Yj,Xj);

>> X=rand(143,4);
>> Y=rand(143,4);
>> C=ones(11,1);
>> for n=0:12
   i=11*n
for j=1:4
   Xj=X(1+i:11+i,j)
   Yj=Y(1+i:11+i,j)
   Xj=[Xj,C]
   [b(:,:,n+1),bint(:,:,n+1),r(:,:,n+1),stats(:,:,n+1)]=regress(Yj,Xj)
end
end

csgt0版主，我按照你的提示做了一遍，可是出錯了，錯誤：“??? Subscripted assignment dimension mismatch.”  能幫我再看一下嗎？謝謝！

贊一下

回復(fù)此樓

4樓2012-08-16 18:40:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

★ ★ ★
wenzimu: 金幣+1, ★★★很有幫助 2012-08-17 14:06:19
xzhdty: 金幣+2, 謝謝 2012-08-17 16:50:48

引用回帖:

4樓: Originally posted by wenzimu at 2012-08-16 18:40:35
>> X=rand(143,4);
>> Y=rand(143,4);
>> C=ones(11,1);
>> for n=0:12
   i=11*n
for j=1:4
   Xj=X(1+i:11+i,j)
   Yj=Y(1+i:11+i,j)
   Xj=
   =regre ...

沒問題啊，原版的

CODE:

X=rand(143,4);

Y=rand(143,4);

C=ones(11,1);

for n=0:12

      i=11*n;

    for j=1:4

        Xj=X(1+i:11+i,j);

        Yj=Y(1+i:11+i,j);

        Xj=[Xj,C];

       [b(:,:,n+1),bint(:,:,n+1),r(:,:,n+1),stats(:,:,n+1)]=regress(Yj,Xj); 

 end 

end

比如b的

b(:,:,1) =

-0.1501
0.5187

b(:,:,2) =

-0.3011
0.6585

b(:,:,3) =

0.0059
0.4830

b(:,:,4) =

-0.3984
0.6927

b(:,:,5) =

0.0696
0.4994

b(:,:,6) =

-0.4565
0.6888

b(:,:,7) =

0.1679
0.5196

b(:,:,8) =

-0.4132
0.8050

b(:,:,9) =

-0.4390
0.7496

b(:,:,10) =

-0.8774
0.9798

b(:,:,11) =

0.2749
0.3699

b(:,:,12) =

0.2595
0.2433

贊一下(1人)

回復(fù)此樓

showmethemoney

5樓2012-08-17 09:33:42

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

b(:,:,13) =

0.3210
0.3065

贊一下

回復(fù)此樓

showmethemoney

6樓2012-08-17 09:34:50

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

wenzimu

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 3
帖子: 21
在線: 7小時
蟲號: 1823664
注冊: 2012-05-18
專業(yè): 人文地理學(xué)

引用回帖:

6樓: Originally posted by csgt0 at 2012-08-17 09:34:50
b(:,:,13) =

0.3210
0.3065

謝謝csgt0版主熱心的答復(fù)。按照您的方法，可以實現(xiàn)了。有個小疑問想再麻煩您一下，就是如果代碼中的Xj，Yj等這些變量也要像b這樣輸出13組值應(yīng)該怎么改代碼呢，因為按照這個運算的話，好像只能得到Xj和Yj最后一組的值，而其他的值得不到。再次感謝，真的太麻煩您了！

贊一下

回復(fù)此樓

7樓2012-08-17 14:55:39

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

引用回帖:

7樓: Originally posted by wenzimu at 2012-08-17 14:55:39
謝謝csgt0版主熱心的答復(fù)。按照您的方法，可以實現(xiàn)了。有個小疑問想再麻煩您一下，就是如果代碼中的Xj，Yj等這些變量也要像b這樣輸出13組值應(yīng)該怎么改代碼呢，因為按照這個運算的話，好像只能得到Xj和Yj最后一組的 ...

把維數(shù)往上加就行了

CODE:

clc

clear

X=rand(143,4);

 Y=rand(143,4);

 C=ones(11,1);

 for n=0:12

       i=11*n;

     for j=1:4

         Xj1=X(1+i:11+i,j);

         Yj(:,j,n+1)=Y(1+i:11+i,j);

         Xj(:,:,j,n+1)=[Xj1,C];

        [b(:,:,n+1),bint(:,:,n+1),r(:,:,n+1),stats(:,:,n+1)]=regress(Yj(:,j,n+1),Xj(:,:,j,n+1)); 

end 

end

贊一下

回復(fù)此樓

showmethemoney

8樓2012-08-17 15:28:37

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

xsx小蝦米

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 15.5
帖子: 4
在線: 7.9小時
蟲號: 2197926
注冊: 2012-12-20
性別: MM
專業(yè): 金融工程

【答案】應(yīng)助回帖

如果不是每11個回歸一次，二是滾動回歸，也就是說1到11行回歸一次，2到12,3-13這樣的代碼該怎么寫？？

贊一下

回復(fù)此樓

9樓2012-12-23 18:06:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 wenzimu 的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 能源材料化學(xué)課題組招收碩士研究生8-10名 +3	脫穎而出 2026-03-16	6/300	2026-03-17 21:19 by z1z2z3879
[考研] 070300化學(xué)319求調(diào)劑 +4	錦鯉0909 2026-03-17	4/200	2026-03-17 18:21 by 重科小霸王
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +4	張vvvv 2026-03-15	6/300	2026-03-17 17:15 by ruiyingmiao
[考研] 材料專碩326求調(diào)劑 +6	墨煜姒莘 2026-03-15	7/350	2026-03-17 17:10 by ruiyingmiao
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 085600材料與化工 +4	安全上岸！ 2026-03-16	4/200	2026-03-17 14:02 by 勇敢太監(jiān)王公公
[考研] 274求調(diào)劑 +5	時間點 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +5	Liwangman 2026-03-15	5/250	2026-03-16 17:10 by 我的船我的海
[考研] 一志愿211 0703方向310分求調(diào)劑 +3	努力奮斗112 2026-03-15	3/150	2026-03-16 16:44 by houyaoxu
[考研] 311求調(diào)劑 +5	26研0 2026-03-15	5/250	2026-03-16 16:21 by a不易
[考研] 327求調(diào)劑 +6	拾光任染 2026-03-15	11/550	2026-03-15 22:47 by 拾光任染
[考研] 0703，333分求調(diào)劑一志愿鄭州大學(xué)-物理化學(xué) +3	李魔女斗篷 2026-03-11	3/150	2026-03-13 22:24 by JourneyLucky
[考研] 求調(diào)劑（材料與化工327） +4	愛吃香菜啦 2026-03-11	4/200	2026-03-13 22:11 by JourneyLucky
[考研] 311求調(diào)劑 +3	冬十三 2026-03-13	3/150	2026-03-13 20:41 by JourneyLucky
[考研] 求調(diào)劑 +5	一定有學(xué)上- 2026-03-12	5/250	2026-03-13 18:31 by ms629
[考研] 314求調(diào)劑 +7	無懈可擊的巨人 2026-03-12	7/350	2026-03-13 15:40 by JourneyLucky
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考研] 290求調(diào)劑 +3	ADT 2026-03-13	3/150	2026-03-13 10:19 by peike
[考研] 283求調(diào)劑，材料、化工皆可 +8	蘇打水7777 2026-03-11	10/500	2026-03-13 09:06 by Linda Hu

24小時熱門版塊排行榜

[求助] matlab中如何用regress實現(xiàn)批量回歸并出圖

» 猜你喜歡

» 本主題相關(guān)價值貼推薦，對您同樣有幫助:

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

【答案】應(yīng)助回帖

» 本主題相關(guān)價值貼推薦，對您同樣有幫助: