0.1 load package

library(Brq)
library(ggplot2)

載入分析需要的套包

0.2 open data

data(USgirl, package="Brq")

0.3 plot age and weight

plot(USgirl, pch='.', bty="n",
     xlim=c(0, 25),
     ylim=c(0, 150),
     xlab="Age (in years)",
     ylab="Body weight (in kg)")

繪製年齡對體重的散佈圖。結果看起來,隨著年齡增加體重有上升的趨勢,但有些發散。 ## Quantile regression

ggplot(USgirl, aes(Age, Weight)) + 
 geom_point(alpha=.5, size=rel(.3)) + 
 geom_quantile(quantiles=1:9/10, 
               formula=y ~ x)+
 labs(x="Age (in years)", 
      y="Body weight (in kg)") +
 theme_minimal()

使用分量迴歸分析,分為九個分量進行繪圖