Hello,

I don't understand what went wrong or how to fix this. How do I set qr=TRUE
for gam?

When I produce a fit using gam like this:

fit = gam(y~s(x),data=as.data.frame(l_yx),family=family,control =
list(keepData=T))

...then try to use predict:
(see #1 below in the traceback() )

> traceback()
6: stop("lm object does not have a proper 'qr' component.\n Rank zero or
should not have used lm(.., qr=FALSE).") at #81
5: qr.lm(object) at #81
4: summary.glm(object, dispersion = dispersion) at #81
3: summary(object, dispersion = dispersion) at #81
2: predict.glm(fit, data.frame(x = xx), type = "response", se.fit = T,
       col = prediction_col, lty = prediction_ln) at #81
1: predict(fit, data.frame(x = xx), type = "response", se.fit = T,
       col = prediction_col, lty = prediction_ln) at #81

...I get this error:

Error in qr.lm(object) : lm object does not have a proper 'qr' component.
 Rank zero or should not have used lm(.., qr=FALSE).

I read this post: http://tolstoy.newcastle.edu.au/R/devel/06/04/5133.html

So I tried adding qr=T to the gam call but it didn't make any difference.
This is how I did it:

fit = gam(y~s(x),data=as.data.frame(l_yx),family=family,control =
list(keepData=T),qr=T)

Its all very strange because I've produced fits with this data may times
before with no issues (and never having to do anything with the qr
parameter. I don't understand why this is coming up or how to fix it.

PS - I don't think this matters, but I am calling a script called
FunctionGamFit.r like this:
err = system(paste('"C:\\Program Files\\R\\R-2.14.1\\bin\\R.exe"', 'CMD
BATCH FunctionGamFit.r'), wait = T)
...to produce the fit.

Thanks for any help!

ben

	[[alternative HTML version deleted]]

