График процентили VGAM в R
Я запускаю следующий код из файлов справки VGAM:
library(VGAM)
fit4 <- vgam(BMI ~ s(age, df = c(4, 2)), lms.bcn(zero = 1), data = bmi.nz, trace = TRUE)
qtplot(fit4, percentiles = c(5,50,90,99), main = "Quantiles", las = 1, xlim = c(15, 90), ylab = "BMI", lwd = 2, lcol = 4)
Как я могу просто запретить размещение точек на графике, чтобы на графике отображались только процентильные кривые? Есть ли какая-то опция в qtplot, где я могу подавить точки на этом графике, чтобы мне не пришлось прибегать к длинному маршруту ggplot, как это сделано на этой странице: Percentiles from VGAM? В моем предыдущем вопросе были и другие проблемы, поэтому этот момент был проигнорирован. Спасибо за вашу помощь.
1 ответ
Решение
Здесь нет qtplot
страница справки, поэтому я пошел к индексу помощи пакета и увидел qtplot.lmscreg
в списке. У него был pcol.arg для управления цветом точек, поэтому я установил его "прозрачным":
qtplot(fit4, percentiles = c(5,50,90,99), main = "Quantiles", las = 1,
xlim = c(15, 90), ylab = "BMI", lwd = 2, lcol = 4,
pcol.arg="transparent")