作者q6261901 (幸福的尾巴)
看板Statistics
標題[程式] [R]quantile regression and Pseudo R2
時間Tue Jan 29 17:08:45 2013
[軟體程式類別]: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
1F：→ bmka:把pseudo R2 的公式找出來,寫一個小function自己算 01/29 20:33
2F：→ bmka:這個function應該會用到fit$resid, fit$tau 01/29 20:34
3F：→ bmka:Professor Koenker的回答在此: 01/29 20:36
4F：→ bmka:r.789695.n4.nabble.com/Pseudo-R-for-Quant-Reg-td805101.ht 01/29 20:36
5F：→ bmka:ml (接上一行) 01/29 20:37
6F：→ q6261901:他寫rho函數 裡面的u一直看不懂是甚麼 如果是rho函數 01/29 22:49
7F：→ q6261901:r裡面就有了 fit$rho 01/29 22:49
8F：→ q6261901:偏偏查來查去 也只能找到相同的回應ˊˋ 01/29 22:50
有查到一個公式 1-(fit$resid/原模型殘差)
^^^^^^^^^^
不懂意思，而且出處非常不可靠...
※ 編輯: q6261901 來自: 111.254.229.9 (01/29 22:54)