代码之家 › 专栏 › 技术社区 › Leosar

使用晶格和quantreg绘制多条分位数回归线

quantreg lattice plot r

Leosar · 技术社区 · 12 年前

我想画几条分位数回归线( quantreg 包装),如下所示:

 library(quantreg)
 data(engel)
 attach(engel)
 plot(income, foodexp, cex = 0.25, type = "n",
      xlab = "Household Income", ylab = "Food Expenditure")
 points(income, foodexp, cex = 0.5, col = "blue")
 abline(rq(foodexp ~ income, tau = 0.5), col = "blue")
 abline(lm(foodexp ~ income), lty = 2, col = "red")
 taus <- c(0.05, 0.1, 0.25, 0.75, 0.9, 0.95)
 for (i in 1:length(taus)) {
   abline(rq(foodexp ~ income, tau = taus[i]),
          col = "gray")
 }
 detach(engel)

但是使用 lattice ,我只尝试了一次回归,但没有成功:

require(lattice)
xyplot(foodexp~income , data =engel, 
   type = c("g"),
   auto.key=list(x=.8,y=.35,cex=.8,cex.title=.8, title="", points=TRUE), 
   scales=list(tck=-1),ylab=list("Food Expenditure",font=3),
   xlab=list("Household Income",font=3),
   panel=function(x,y){
     panel.xyplot(x,y)
     panel.grid()
     panel.abline(rq(foodexp ~ income, tau = 0.5))}
)

1 回复 | 直到 9 年前

agstudy 12 年前

使用 lattice 你应该使用 y 和 x 在面板功能内。

enter image description here

xyplot(foodexp~income , data =engel, 
       type = c("g"),
       auto.key=list(x=.8,y=.35,cex=.8,cex.title=.8, title="", points=TRUE), 
       scales=list(tck=-1),ylab=list("Food Expenditure",font=3),
       xlab=list("Household Income",font=3),
       panel=function(x,y,...){
         panel.xyplot(x,y)
         panel.grid()
         panel.abline(rq(y ~ x, tau = 0.5))
         panel.points(x, y, cex = 0.5, col = "blue")
         panel.abline(rq(y ~ x, tau = 0.5), col = "blue")
         panel.abline(lm(y ~ x), lty = 2, col = "red")
         taus <- c(0.05, 0.1, 0.25, 0.75, 0.9, 0.95)
         for (i in 1:length(taus)) {
           panel.abline(rq(y ~ x, tau = taus[i]),
                  col = "gray")
         }

       }
)

推荐文章

Vint · 时间序列级图上的覆盖点

8 年前

LuckyProgrammer · 使用as绘制splom()图形。数据框架

8 年前

Geo-sp · 在水平图面板区域外添加文本

10 年前

degopwn · 自定义图例R绘图(晶格)

11 年前

Carlos Cinelli · 模型的三维绘图(三维散点图+模型表面+连接点到表面)

11 年前

orizon · 组合点阵图在x轴上产生空间隙

11 年前

Leosar · 使用晶格和quantreg绘制多条分位数回归线

12 年前

user2123601 · 无法标记点图中的点

12 年前