作者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)