Hi,

I am using the gam function in the mgcv package, I have random effects in
my model (bs="re") this has worked fine, but after I updated the mgcv
package to version 1.7-19 I recive an error message when I run the model.

>
fit1<-gam(IV~s(RUTE,bs="re")+s(T13)+s(H40)+factor(AAR)+s(V3)+s(G1)+s(H1)+s(V1)+factor(LEDD),data=data5,method="ML")
> summary.gam(fit1)
Error in mroot(Ve[ind, ind]) : Supplied matrix not symmetric

This same model worked fine before I updated the mgcv package, does anyone
know why I now get the error message?

Another question I have is if anyone knows how I can compare two gam models
that includes "re" terms? As far as I understand an anova is not the best
way when I have "re" terms.

Best regards
Silje Skår

	[[alternative HTML version deleted]]

