Plot Quantile-Quantile plot with qqline using ggplot2

image

1
2
3
4
5
6
myresiduals <- rnorm(100)
pdf("01.f4.pdf",width=3,height=3)
ggplot(data=as.data.frame(qqnorm( myresiduals , plot=F)), mapping=aes(x=x, y=y)) + geom_point() + geom_smooth(method="lm", se=FALSE)+labs(title="D(Agarli,Urial;Mouflon1,Mouflon2)",x="Expected Z-scores",y="Observed Z-scores")+ theme_classic()
dev.off()

Quantile-Quantile plot can use to analysis D-statistics like this:
image
(DOI: 10.1126/science.aan3842 , Fig. 2)