mfasiolo / qgam

Additive quantile regression R package
http://mfasiolo.github.io/qgam/
30 stars 7 forks source link

Memory error in PIRLS with zero weights? #25

Open mfasiolo opened 6 years ago

mfasiolo commented 6 years ago

This code:

library(qgam)
data("UKload")
form <- NetDemand~s(wM,k=20,bs='cr') + s(wM_s95,k=20,bs='cr') +
  s(Posan, k=30, bs = "ad") + Dow + s(Trend,k=4,bs='cr') +
  NetDemand.48 + Holy

qu <- 0.7
fit <- qgam(form = form, data = UKload, qu = qu, err = 0.05)

occasionally generates a warning. But sometimes it doesn't. Given that nothing is random here, this makes me think that there might be a memory error in the C code. This does not happen if we set err = 0.1.