2013-02-20 44 views
0

我在R中使用quantreg包對一組數據運行分位數迴歸(95%)。分位數迴歸中的測試斜率

我想將分位數迴歸的斜率設置爲值1.4,就像之前的分析中所做的那樣,我想比較我的結果。如果在lm()這對於函數offset()是可能的,對於固定分位數(例如0.025)使用rq(),則這不起作用。

該代碼不會給出錯誤,但1.4的值對我的結果沒有影響。

fit.0.025<-rq(y~offset(1.4*x),tau=0.025, data=mydataframe) 

回答

0

這不應該是這樣的測試代碼沒有得到充分測試嗎? (除非攔截列被抑制,否則它仍然是一個需要估計的參數,從而破壞了修復它的努力)。(編輯,需要重複偏移並使用-1代替+0)

dfrm <- data.frame(x=runif(1000, 1,100), y=runif(1000, 1,2)) 
    fit.0.025 < -rq(y ~ x*(1 + offset(rep(1.4, 1000))) -1, tau=0.025, data=dfrm) 

# Same as: 
fit.0.025<-rq(y ~ x + offset(rep(1.4, 1000)) -1 , tau=0.025, data=dfrm) 

很坦白地說,我在質疑這是否有任何統計意義。運行數學運算並不總是產生可解釋的輸出。

我前面鍵入:

fit.0.025<-rq(y ~ x+ offset(1.4), tau=0.025, data=mydataframe) 

....但這隻會是適當的乘法鏈接。