This is a bug. The `prior.weights' element of the the faked `gam' object
(i.e. `test$gam$prior.wieghts' below) has been set to the varIdent()
variance function, rather than the weights that this eventually
represents. I'll fix this for the next patch release, (as soon as I get
any time to do
Hello,
Why would I get an error message with the following code for gamm? I
want to fit the a gam with different variances per stratum.
library(mgcv)
library(nlme)
Y-rnorm(100)
X-rnorm(100,sd=2)
Z-rep(c(T,F),each=50)
test-gamm(Y~s(X),weights=varIdent(form=~1|Z))
summary(test$lme) #ok