R語言：畫q-q圖

1. 首先安裝畫圖的package叫做qqman：

install.packages('qqman')
library(qqman)

2. 導入需要繪圖的數據：

data1 <-read.table(file.choose(),header=T,sep="")

data2 <-read.table(file.choose(),header=T,sep="")

data3<- read.table(file.choose(),header = T,sep = "")

3. 將需要作圖的未校正P值變量（如UNADJ）拿來作圖即可：

qq(data1$UNADJ)

如果想把三幅圖放到一起橫向比較，則：

par(mfrow=c(1,3))
qq(dataage1$UNADJ)
qq(dataage2$UNADJ)
qq(dataage3$UNADJ)

下面展示一副我的圖：