求助R分位数回归作图

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我在做分位数回归的时候，由于自变量较多，使用plot(summary(fit))的时候，一张图分布太密，想把这些图放到多张纸上，并且让他显示的更加精确一些，如何做？
fit<-rq(y ~x1+x2+x3+x4+x5+x6,tau=2:98/100,data=data)
summary(fit)
plot(summary(fit))

一张上的图太多，如何设置每张显示2*2的图片，谢谢！

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：分位数回归 分位数 Summary Summa Mary 自变量 如何

plot(summary(fit),mar = c(5.1, 4.1, 2.1, 2.1), parm = c(1,2,3,4), main = "", xlab = "tau",
  ylab =c("intercept","Air.Flow","Water.Temp","Acid.Conc"), cex = 1, pch = 19)

谢谢！！解决了！！

mar = c(5.1, 4.1, 2.1, 2.1),代表什么意思呀

我要显示剩余的的那几个变量的图怎么办呢？？？

parm = c(5,6,7)
修改这个就可以，其他不用动

ermutuxia 发表于 2015-1-27 18:24
parm = c(5,6,7)
修改这个就可以，其他不用动

高手 可不可以帮我解答个问题  就是我进行分位数回归
> fit1<-rq(y~x1+x2+x3+x4+x5+x6+x7,tau=c(0.1,0.25,0.5,0.75,0.95),data=rt)
> summary(fit1,se='boot')

这些都能出结果    我也想画和楼主一样的图  就用了命令plot(summary(fit1))   但是却显示错误
Error in plot.window(...) : infinite axis extents [GEPretty(-inf,inf,5)]   我也不知道哪错了？？？ 同样的命令  为什么我的就错了  求高手指点  谢谢

plot(summary(fit),mar = c(5.1, 4.1, 2.1, 2.1), se=`boot')试试

ermutuxia 发表于 2015-2-10 09:17
plot(summary(fit),mar = c(5.1, 4.1, 2.1, 2.1), se=`boot')试试

请问  是不是和数据数量有关呢？我用别人的例子做（有300组数值）没有问题  而我的只有20组数据  是不是因为数据太少不能做呢

ermutuxia 发表于 2015-2-10 09:17
plot(summary(fit),mar = c(5.1, 4.1, 2.1, 2.1), se=`boot')试试

我按照你的命令做 结果
> plot(summary(fit),mar = c(5.1, 4.1, 2.1, 2.1), se=`boot')
+
自动弹出“+”  就好像命令没完成一样


非常感谢