Re: [程式] R分量迴歸結果顯示

看板Statistics作者 (Syameroke)時間11年前 (2013/01/03 20:56), 編輯推噓3(309)
留言12則, 3人參與, 最新討論串2/2 (看更多)
※ 引述《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)

01/03 21:00, , 1F
資料儲存在"X1變數"這個物件內
01/03 21:00, 1F

01/03 21:31, , 2F
謝謝 再請教一個問題 如果想輸出標準誤
01/03 21:31, 2F

01/03 21:31, , 3F
該將coefficients改成什麼
01/03 21:31, 3F

01/04 11:01, , 4F
手上沒R,你用ls(rq(y~x,0.5))去試試看裡面有哪些
01/04 11:01, 4F

01/04 11:01, , 5F
應該會看到標準誤是哪一個名字
01/04 11:01, 5F

01/04 12:45, , 6F
std 但輸出出來是null 我再研究看看好了 謝謝
01/04 12:45, 6F

01/04 15:36, , 7F
我查了一下,rq(y~x,0.5)裡面似乎沒有標準誤,而
01/04 15:36, 7F

01/04 15:36, , 8F
summary(rq(y~x,0.5))裡面也只有上下界,沒有標準誤
01/04 15:36, 8F

01/04 15:37, , 9F
下界是這樣叫的:
01/04 15:37, 9F

01/04 15:38, , 10F
summary(rq(y~x,0.5))$coefficients[2,2]
01/04 15:38, 10F

01/04 15:38, , 11F
這是X1的下界,你自己在for迴圈中計算標準誤吧
01/04 15:38, 11F

01/04 18:52, , 12F
謝謝!!
01/04 18:52, 12F
文章代碼(AID): #1GvO0H_X (Statistics)
文章代碼(AID): #1GvO0H_X (Statistics)