分位数を使用して共分散の分析を行っています。グループごとに分位線を使用してプロットしたいと思います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 レベルでどの線が有意であるかを知る方法がわかりません。