作者ching0629 (Syameroke)
看板Statistics
標題Re: [程式] R分量迴歸結果顯示
時間Thu Jan 3 20:56:47 2013
※ 引述《Ennail5566 (小妹)》之銘言：
: [軟體程式類別]:
: R
: [程式問題]:分量迴歸結果輸出
: [軟體熟悉度]:
: 低(1~3個月)
: [問題敘述]:
: 目前我可以將我想輸出分量迴歸的表顯示出，
: 如:
: 0.01分量
: 變數 value std. error .....
: x1 . . .....
: x2 . . .....
: x3 . . .....
: 但我希望單獨顯示所有分量個別的係數
: 如:
: x1變數
: 分量 value
: 0.01 3
: 0.05 -2
: 0.10 -1
: 0.25 -10
: . .
: . .
: . .
: 請問我該輸入什麼指令
: 感謝幫忙
A=c(0.01,0.05,0.1,0.25,0.5,0.75,0.9,0.95,0.99)
X1變數=matrix(NA,nrow=length(A),ncol=2)
for (分量 in 1:length(A)) {
X1變數[分量,1]=A[分量]
X1變數[分量,2]=rq(y~x,A[分量])$coefficients[2]
}
colnames(X1變數)=c("分量","value")
其中A是分量的總數(我想你應該看的懂)
第五行的y與x擺上自己想擺的,但x1(或是該自變數)要放第一個
由於你給你資訊很少,我甚至不知道你的自變項共有多少個,所以你就按照自己需要的呈現
吧
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 27.105.2.50
※ 編輯: ching0629 來自: 27.105.2.50 (01/03 20:59)
1F：→ ching0629:資料儲存在"X1變數"這個物件內 01/03 21:00
2F：推 Ennail5566:謝謝 再請教一個問題 如果想輸出標準誤 01/03 21:31
3F：→ Ennail5566:該將coefficients改成什麼 01/03 21:31
4F：→ ching0629:手上沒R,你用ls(rq(y~x,0.5))去試試看裡面有哪些 01/04 11:01
5F：→ ching0629:應該會看到標準誤是哪一個名字 01/04 11:01
6F：推 q6261901:std 但輸出出來是null 我再研究看看好了 謝謝 01/04 12:45
7F：→ ching0629:我查了一下,rq(y~x,0.5)裡面似乎沒有標準誤,而 01/04 15:36
8F：→ ching0629:summary(rq(y~x,0.5))裡面也只有上下界,沒有標準誤 01/04 15:36
9F：→ ching0629:下界是這樣叫的: 01/04 15:37
10F：→ ching0629:summary(rq(y~x,0.5))$coefficients[2,2] 01/04 15:38
11F：→ ching0629:這是X1的下界,你自己在for迴圈中計算標準誤吧 01/04 15:38
12F：推 q6261901:謝謝!! 01/04 18:52