R help - Nov 2012 - Plot with residuals in mgcv

Hi,

I am using the mgcv package (version 1.7-22.) running the model works fine,
but when I want to have a plot with residuals I get an error.

fit29<-gam(IV~s(G3)+s(V3)+factor(AAR)+s(D3)+s(RUTE,bs="re"),data=subsf,gamma=1.4,method="ML")


plot(fit29,residuals=T)
Error in X[, first:last] %*% object$coefficients[first:last] :
  non-conformable arguments

does some one know what this error means? the subsf matrix is 35x27.

Silje

	[[alternative HTML version deleted]]

On 2012-11-28 05:02, silje sk?r wrote:
> Hi,
>
> I am using the mgcv package (version 1.7-22.) running the model works fine,
> but when I want to have a plot with residuals I get an error.
>
>
fit29<-gam(IV~s(G3)+s(V3)+factor(AAR)+s(D3)+s(RUTE,bs="re"),data=subsf,gamma=1.4,method="ML")
>
>
> plot(fit29,residuals=T)
> Error in X[, first:last] %*% object$coefficients[first:last] :
>    non-conformable arguments
>
> does some one know what this error means? the subsf matrix is 35x27.
>
> Silje
Well, the error means that the matrix X[, first:last] and
the vector object$coefficients[first:last] don't have
'matching' dimensions. I don't see why, but since the
'residuals' argument to plot.gam() can be an array
_of the correct length_ (see ?plot.gam), it just might be
that you have an object called 'T' hanging around, in which
case using 'T' in place of 'TRUE' is a bad idea.
Actually, it's _always_ a bad idea.

Peter Ehlers

Hi Silje,

Thanks for this. I guess RUTE is a numeric variable in this model and 
hence only has one associated random coefficient? This then causes a 
problem when calling predict.gam as part of processing 'residuals=T'. 
I've fixed the problem for the next release, but did you really want a 
single random slope for RUTE, or should it really have been declared as 
a factor variable? btw gamma doesn't do anything with ML and REML 
smoothness selection (it's only used in prediction error criteria 
smoothness selection to increase the penalization per effective degrees 
of freedom).

best,
Simon

On 28/11/12 13:02, silje sk?r wrote:
> Hi,
>
> I am using the mgcv package (version 1.7-22.) running the model works fine,
> but when I want to have a plot with residuals I get an error.
>
>
fit29<-gam(IV~s(G3)+s(V3)+factor(AAR)+s(D3)+s(RUTE,bs="re"),data=subsf,gamma=1.4,method="ML")
>
>
> plot(fit29,residuals=T)
> Error in X[, first:last] %*% object$coefficients[first:last] :
>    non-conformable arguments
>
> does some one know what this error means? the subsf matrix is 35x27.
>
> Silje
>
> 	[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.