所以我试图计算分位数回归并在ggplot
中绘制结果由于某些原因,在ggplot
]中绘制结果时,由于某些原因,无法显示虚拟变量
示例mtcars
数据集上的代码,如下所示:
library(dplyr) library(ggplot2) library(qre) library(quantreg) library(fastDummies) dataset <- mtcars dataset <- dummy_cols(dataset, select_columns = "gear") dataset rq(data=dataset, tau= 1:9/10, formula = hp ~ disp + mpg + qsec + gear_4 + gear_5) %>% broom::tidy() %>% #filter(term!="(Intercept)") %>% ggplot(aes(x=tau,y=estimate))+ geom_point(color="#27408b", size = 3)+ geom_ribbon(aes(ymin=conf.low,ymax=conf.high),alpha=0.25, fill="#27408b")+ geom_line(color="#27408b", size = 1)+ geom_smooth(method= "lm", colour = "red", se = T)+ my_theme + facet_wrap(~term,scales="free",ncol=2) QR.2 <- rq(hp ~ disp + mpg + qsec + gear_4 + gear_5, data = dataset, tau = 1:9/10) plot(summary(QR.2, se = "boot"))
[使用
plot(summary(QR.2, se = "boot"))
绘制结果时,一切都很好,但是由于某些原因,使用ggplot会显示错误。
因此,出于某种原因,我试图在ggplot中计算分位数回归并绘制结果,由于在某些示例mtcars上的ggplot代码中绘制结果而由于某些原因而无法显示伪变量...
ggplot不会为两个gear
列绘制任何内容,因为至少一个分位数具有基本为无限的标准误差范围。请注意,在下面的输出中,对于每个.Machine$double.xmax
级别,至少一个分位数的最大绝对值约为1.8E308(gear
)。