我想绘制分位数回归的结果,但我无法:
这是我的代码:
require(quantreg)
data(engel)
attach(engel)
xx <- income - mean(income)
zz <- c(120, diff(income))
fit1 <- summary(rq(foodexp~xx+zz, tau=2:98/100))
然后:
png('res.png')
plot(fit1, mfrow=c(1,2))
只有zz图被保存到res.png文件中。有什么办法可以将图保存在单独的文件中(两个和一个)? 以及如何控制图的宽度/高度?当我将它们保存到.png文件时,我喜欢所有单个图的宽度=高度(正方形)?
答案 0 :(得分:24)
您可以按png
参数控制图像尺寸。
png("image.png", width = 800, height = 600)
plot(...)
dev.off()
要“完成”图片,请使用dev.off
。
答案 1 :(得分:5)
用于细分图:
plot(fit1,parm=1:2)
plot(fit1,parm=3)
请注意,您可以通过仔细阅读?plot.summary.rqs
找到答案,但这可能不是很明显:为了知道在哪里看,您需要class(fit1)
来确定正在使用plot
方法。
罗曼的回答照顾了图像维度。