r - 编程QQ图

我有一个男女生数学考试成绩的样本。我想为每个性别绘制QQ图，以查看每个性别是否均呈正态分布。我知道如何绘制整个样本的QQ图，但是如何分别绘制它们呢？

最佳答案

这是使用base图形的简单解决方案:

scores <- rnorm(200, mean=12, sd=2)
gender <- gl(2, 50, labels=c("M","F"))
opar <- par(mfrow=c(1,2))
for (g in levels(gender))
  qqnorm(scores[gender==g], main=paste("Gender =", g))
par(opar)

然后，一个更优雅的lattice解决方案:

qqmath(~ scores | gender, data=data.frame(scores, gender), type=c("p", "g"))

有关更多讨论和可能的自定义示例，请参见qqmath的在线帮助。