ggplot中带有虚拟变量的分位数回归错误

问题描述 投票:0回答:1

所以我试图计算分位数回归并在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代码中绘制结果而由于某些原因而无法显示伪变量...

r ggplot2 quantile-regression
1个回答
0
投票

ggplot不会为两个gear列绘制任何内容,因为至少一个分位数具有基本为无限的标准误差范围。请注意,在下面的输出中,对于每个.Machine$double.xmax级别,至少一个分位数的最大绝对值约为1.8E308(gear)。

© www.soinside.com 2019 - 2024. All rights reserved.