123456myresiduals <- 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:
(DOI: 10.1126/science.aan3842 , Fig. 2)