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)
下面展示一副我的图: