quantreg

3熱度

1回答

我有一個大約250年的最大降雨量測量數據系列，maxima [，]，並且想要一次對所有系列應用分位數迴歸，並獲得每個迴歸模型的顯着性R. library(quantreg) qmag <- array(NA, c(250,4)) taus <- c(0.05, 0.1, 0.95, 0.975) for(igau in 1:250){ qure <-

4熱度

1回答

如何通過引導默認分位數迴歸來獲得置信區間

在我的統計類中，我們使用Stata，由於我是R用戶，我想在R中做同樣的事情。我得到了正確的結果，但似乎就像獲得像置信區間這樣簡單的事情有點尷尬。這裏是我的粗液： library(quantreg) na = round(runif(100, min=127, max=144)) f <- rq(na~1, tau=.5, data=ds) s <- summary.rq(f, se="bo

0熱度

1回答

分位數迴歸中的測試斜率

我在R中使用quantreg包對一組數據運行分位數迴歸（95％）。我想將分位數迴歸的斜率設置爲值1.4，就像之前的分析中所做的那樣，我想比較我的結果。如果在lm()這對於函數offset()是可能的，對於固定分位數（例如0.025）使用rq()，則這不起作用。該代碼不會給出錯誤，但1.4的值對我的結果沒有影響。 fit.0.025<-rq(y~offset(1.4*x),tau=0.025,

2熱度

2回答

如何在R中使用texreg進行分位數迴歸時設置se = boot？

我正在運行分位數迴歸（包quantreg）並使用texreg來創建我的模型的乳膠輸出。我感興趣的是bootstrapped s.e.並在摘要的選項中設置se =「boot」，但是當我使用texreg時，我得到「n.i.d.」 S.E. 如何更改該選項？下面是我在做什麼： tm3 <- rq(nback ~ cara + mat + dut + e_brown + e_green + e_bl

1熱度

1回答

LAPACK不掛編譯軟件包時，quantreg

包quantreg無法安裝，因爲它沒有鏈接到LAPACK： install.packages("quantreg") [剪斷] gcc -std=gnu99 -shared -o quantreg.so akj.o boot.o brute.o chlfct.o cholesky.o combos.o crq.o crqfnb.o dsel05.o etime.o extract.o idmi

1熱度

1回答

[R quantreg在協

我有我的quantreg包的R如下運行位數迴歸的每個級別獲得的效果： read19 <- rq(comppercentile~q19*mothered+ethnicity+gender,tau=.5,data=medmodel) Q19是定量和地照顧是一個2級因素。一旦我在q19和mothered之間得到了顯着的互動，我如何在每個motled的水平上得到q19的效果？