1

I am running following code from help files of 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) 

enter image description here

How can I just prevent putting points on the plot so that the graph shows only percentile curves? Is there some option in qtplot where I can suppress points on this graph, so that I do not need to resort to long ggplot route as done on this page: Percentiles from VGAM ? In my earlier question there were other issues also so this point got ignored. Thanks for your help.

Community
  • 1
  • 1
rnso
  • 23,686
  • 25
  • 112
  • 234

1 Answers1

0

There is no qtplot help page so I went to the package help Index and saw qtplot.lmscreg listed. It had a 'pcol.arg' to control points color so I set it to "transparent":

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")

enter image description here

IRTFM
  • 258,963
  • 21
  • 364
  • 487
  • I am surprised they did not keep the option of not putting points on plot. Thanks. – rnso Dec 19 '14 at 01:09
  • I was puzzled about that as well. Also puzzled by no help page for a generic function. The class of the fit4 object did not match any of the qtplot methods in the Index so I just took a guess that since you were making it with lms.bcn that the nearest match to a method name might have a listing for parameters. – IRTFM Dec 19 '14 at 01:12
  • Please help me in another question on vgam: http://stackoverflow.com/questions/27558687/how-to-suppress-vgam-loglikelihood-output-in-r – rnso Dec 19 '14 at 01:54
  • You should be embarrassed. – IRTFM Dec 19 '14 at 02:04
  • I had definitely tried this option, but at that time it did not work. I may have misspelled something. I have deleted that question. – rnso Dec 19 '14 at 02:19