wenzimu

新蟲 (初入文壇)

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

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

現(xiàn)在有三組143*4的矩陣數(shù)據(jù)，分別命名為GDP，LGDP，GRO�，F(xiàn)在想通過regress以GDP為因變量，LGDP,GRO為自變量進(jìn)行回歸。如果是整體回歸，那直接調(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)這樣分別進(jìn)行批量回歸，并畫出實(shí)際值GDP與擬合值GDP的圖；然后第二列、第三列、第四列數(shù)據(jù)也同樣進(jìn)行這樣的擬合（如下圖所示），不知這樣該如何實(shí)現(xiàn)，還請高人幫忙~!

示意圖

回復(fù)此樓

» 猜你喜歡

311求調(diào)劑已經(jīng)有9人回復(fù)
302求調(diào)劑已經(jīng)有4人回復(fù)
被我言中：新模板不強(qiáng)調(diào)格式了，假專家開始管格式了已經(jīng)有4人回復(fù)
303求調(diào)劑已經(jīng)有5人回復(fù)
考研求調(diào)劑已經(jīng)有4人回復(fù)
085601專碩，總分342求調(diào)劑，地區(qū)不限已經(jīng)有4人回復(fù)
085601材料工程專碩求調(diào)劑已經(jīng)有5人回復(fù)
一志愿天津大學(xué)化學(xué)工藝專業(yè)（081702）315分求調(diào)劑已經(jīng)有6人回復(fù)
能源材料化學(xué)課題組招收碩士研究生8-10名已經(jīng)有6人回復(fù)
296求調(diào)劑已經(jīng)有5人回復(fù)

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

關(guān)于matlab中循環(huán)讀取多個(gè)excel文件的問題？已經(jīng)有3人回復(fù)
如何用SAS 進(jìn)行回歸系數(shù)的顯著性檢驗(yàn)。。。。。已經(jīng)有3人回復(fù)
DFT Scan之后使用MP2計(jì)算能量，如何建立批量輸入文件？已經(jīng)有6人回復(fù)
matlab中如何將雙坐標(biāo)系圖中的另外一個(gè)坐標(biāo)系設(shè)置為當(dāng)前坐標(biāo)系呢？已經(jīng)有10人回復(fù)
求助：非線性方程系數(shù)回歸問題？已經(jīng)有10人回復(fù)
非線性回歸參數(shù)..急求幫助��！謝謝已經(jīng)有8人回復(fù)
求助：哪位大俠懂多重logistic回歸的指導(dǎo)我分析一個(gè)數(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小時(shí)
蟲號(hào): 1823664
注冊: 2012-05-18
專業(yè): 人文地理學(xué)

我自己用隨機(jī)數(shù)做了一下，但最后發(fā)現(xiàn)回歸參數(shù)不能一起出現(xiàn)在同一個(gè)矩陣中，最后矩陣中只有某一組的回歸參數(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

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

【答案】應(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小時(shí)
蟲號(hào): 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版主，我按照你的提示做了一遍，可是出錯(cuò)了，錯(cuò)誤：“??? Subscripted assignment dimension mismatch.”  能幫我再看一下嗎？謝謝！

贊一下

回復(fù)此樓

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

【答案】應(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

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

【答案】應(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小時(shí)
蟲號(hào): 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ù)。按照您的方法，可以實(shí)現(xiàn)了。有個(gè)小疑問想再麻煩您一下，就是如果代碼中的Xj，Yj等這些變量也要像b這樣輸出13組值應(yīng)該怎么改代碼呢，因?yàn)榘凑者@個(gè)運(yùn)算的話，好像只能得到Xj和Yj最后一組的值，而其他的值得不到。再次感謝，真的太麻煩您了！

贊一下

回復(fù)此樓

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

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

csgt0

榮譽(yù)版主 (著名寫手)

彩色掛圖

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

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

引用回帖:

7樓: Originally posted by wenzimu at 2012-08-17 14:55:39
謝謝csgt0版主熱心的答復(fù)。按照您的方法，可以實(shí)現(xiàn)了。有個(gè)小疑問想再麻煩您一下，就是如果代碼中的Xj，Yj等這些變量也要像b這樣輸出13組值應(yīng)該怎么改代碼呢，因?yàn)榘凑者@個(gè)運(yùn)算的話，好像只能得到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小時(shí)
蟲號(hào): 2197926
注冊: 2012-12-20
性別: MM
專業(yè): 金融工程

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

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

贊一下

回復(fù)此樓

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

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

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

[考研] 268求調(diào)劑 +7	好運(yùn)連綿不絕 2026-03-12	8/400	2026-03-17 20:28 by xilongliang
[考研] 277調(diào)劑 +5	自由煎餅果子 2026-03-16	6/300	2026-03-17 19:26 by 李leezz
[考研] 301求調(diào)劑 +4	A_JiXing 2026-03-16	4/200	2026-03-17 17:32 by ruiyingmiao
[考研] 312求調(diào)劑 +4	陌宸希 2026-03-16	5/250	2026-03-17 17:09 by ruiyingmiao
[考研] 26考研求調(diào)劑 +6	丶宏Sir 2026-03-13	6/300	2026-03-17 16:13 by 醉在風(fēng)里
[考研] 梁成偉老師課題組歡迎你的加入 +8	一鴨鴨喲 2026-03-14	10/500	2026-03-17 15:07 by 一鴨鴨喲
[考研] 290求調(diào)劑 +6	孔志浩 2026-03-12	11/550	2026-03-17 14:41 by 周舟舟77
[考研] 278求調(diào)劑 +3	Yy7400 2026-03-13	3/150	2026-03-17 08:24 by laoshidan
[考研] [導(dǎo)師推薦]西南科技大學(xué)國防/材料導(dǎo)師推薦 +3	尖角小荷 2026-03-16	6/300	2026-03-16 23:21 by 尖角小荷
[考研] 環(huán)境工程調(diào)劑 +6	大可digkids 2026-03-16	6/300	2026-03-16 17:16 by barlinike
[考研] 一志愿華中師范071000，325求調(diào)劑 +6	RuitingC 2026-03-12	6/300	2026-03-16 14:50 by 可淡不可忘
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7	膩膩gk 2026-03-14	7/350	2026-03-16 10:12 by houyaoxu
[考研] 機(jī)械專碩調(diào)劑 +3	笨笨兔子 2026-03-12	3/150	2026-03-15 20:02 by 栗子粥?
[考研] 297一志愿上交085600求調(diào)劑 +5	指尖八千里 2026-03-14	5/250	2026-03-14 17:26 by a不易
[考研] 265求調(diào)劑 +4	威化餅07 2026-03-12	4/200	2026-03-14 17:23 by userper
[考研] 復(fù)試調(diào)劑 +4	z1z2z3879 2026-03-14	5/250	2026-03-14 16:30 by JourneyLucky
[考研] 材料專碩288分求調(diào)劑一志愿211 +4	在家想你 2026-03-11	4/200	2026-03-13 22:49 by JourneyLucky
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 307求調(diào)劑 +5	超級(jí)伊昂大王 2026-03-12	5/250	2026-03-13 15:56 by 棒棒球手
[考研] 求調(diào)劑資源與環(huán)境 285 +3	未名考生 2026-03-10	3/150	2026-03-13 10:31 by houyaoxu