[程式] [R]quantile regression and Pseudo R2

看板Statistics作者 (幸福的尾巴)時間13年前 (2013/01/29 17:08), 編輯推噓0(008)
留言8則, 2人參與, 最新討論串1/1
[軟體程式類別]:R [程式問題]:Pseudo R2 計算 [軟體熟悉度]: 低(1~3個月) [問題敘述]: 想把分量迴歸中的Pseudo R2計算出來, 在網路上找很多資料 但就是找不到... 完全不抱任何希望上來問問看> < [程式範例]: xx <- read.table("E:/rdata.csv", header = TRUE, sep = ",") qr1 <- function(aa,bb){ fit <- rq(LNSALES~SOURCE+ENTROPY+INASS+LEVERAGE+EFFICENC+ MD1+MD2+MD3+MD4+MD5+RAILWAY+ROAD+WAT+ELE, tau=aa, data=bb) } fit1<-qr1(1:99/100,xx) result1<-summary(fit1) -- ※ 發信站: 批踢踢實業坊(ptt.cc) ◆ From: 111.254.229.9

01/29 20:33, , 1F
把pseudo R2 的公式找出來,寫一個小function自己算
01/29 20:33, 1F

01/29 20:34, , 2F
這個function應該會用到fit$resid, fit$tau
01/29 20:34, 2F

01/29 20:36, , 3F
Professor Koenker的回答在此:
01/29 20:36, 3F

01/29 20:36, , 4F
r.789695.n4.nabble.com/Pseudo-R-for-Quant-Reg-td805101.ht
01/29 20:36, 4F

01/29 20:37, , 5F
ml (接上一行)
01/29 20:37, 5F

01/29 22:49, , 6F
他寫rho函數 裡面的u一直看不懂是甚麼 如果是rho函數
01/29 22:49, 6F

01/29 22:49, , 7F
r裡面就有了 fit$rho
01/29 22:49, 7F

01/29 22:50, , 8F
偏偏查來查去 也只能找到相同的回應ˊˋ
01/29 22:50, 8F
有查到一個公式 1-(fit$resid/原模型殘差) ^^^^^^^^^^ 不懂意思,而且出處非常不可靠... ※ 編輯: q6261901 來自: 111.254.229.9 (01/29 22:54)
文章代碼(AID): #1H1v6X7i (Statistics)