Программирование сюжета QQ

Question

У меня есть образец результатов теста по математике для студентов мужского и женского пола.Я хочу нарисовать график QQ для каждого пола, чтобы увидеть, распределен ли каждый из них нормально.Я знаю, как нарисовать график QQ для всего образца, но как я могу нарисовать их отдельно?

Answer 1

Вот простое решение с использованием 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 для более подробного обсуждения и примера возможной настройки.