График процентиля 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, как это сделано на этой странице: Процентили от VGAM ? В моем предыдущем вопросе были и другие проблемы, поэтому этот момент был проигнорирован. Спасибо за вашу помощь.


r quantile curve-fitting vgam
person rnso    schedule 18.12.2014    source источник

Ответы (1)


arrow_upward
0
arrow_downward

Страницы справки 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")

введите здесь описание изображения

person IRTFM    schedule 18.12.2014
comment
Я удивлен, что они не сохранили возможность не ставить точки на сюжете. Спасибо. - person rnso; 19.12.2014
comment
Я тоже был озадачен этим. Также озадачен отсутствием справочной страницы для универсальной функции. Класс объекта fit4 не соответствует ни одному из методов qtplot в индексе, поэтому я просто предположил, что, поскольку вы делаете это с помощью lms.bcn, ​​ближайшее совпадение с именем метода может иметь список параметров. - person IRTFM; 19.12.2014
comment
Пожалуйста, помогите мне в другом вопросе о vgam: stackoverflow.com/questions/27558687/ - person rnso; 19.12.2014
comment
Вам должно быть стыдно. - person IRTFM; 19.12.2014
comment
Я определенно пробовал этот вариант, но в то время он не работал. Я может что-то неправильно написал. Я удалил этот вопрос. - person rnso; 19.12.2014