library(qgam); library(MASS)
quSeq <- c(0.4, 0.6)
set.seed(737)
fit <- mqgam(accel~s(times, k=20, bs="ad"), data = mcycle, qu = quSeq)
qdo(fit, 0.4, summary)
invisible(qdo(fit, 0.4, plot, pages = 1))
# Return the object for qu = 0.6 and then plot it
tmp <- qdo(fit, 0.6)
plot(tmp)
Run the code above in your browser using DataLab