我正在使用分位数进行协方差分析,我想ggplot2
按组使用分位数线进行绘图,但如果它们在某个水平上显着,例如 p<0.1,则更改样式
library(quantreg)
library(ggplot2)
data(barro)
barro$group=factor(round(runif(nrow(barro),1,3)))
fit0 <- rq(y.net ~ lgdp2 * group , data = barro,tau=c(.25, .5, .75))
sf <-summary(fit0,se="boot")
#extract coefficients
as.data.frame(lapply(sf, coef))
ggplot(barro,aes(lgdp2,y.net)) +facet_wrap( ~ group)+ geom_point(aes(colour=group)) + theme_bw() + stat_quantile()
在这个图中,所有线条都具有相同的样式,我不知道如何更改分位数线条的样式,以了解它们中的哪些在指定的 p 级别上是有意义的。