Martijn,
It's a bit hard to know without seeing the full model structure, but
it's possible that the issue is related to an undesirable side effect of
the handling of identifiability constraints on smooth terms, prior to
mgcv 1.7-13: the standard side constraint approach used for smooths
could lead to unexpected constraints being applied to s(...,bs="re")
terms in some cases.
So, could you sent me the gam call that generates the problem, and
perhaps try out if it still happens in 1.7-13?
best,
Simon
On 23/04/12 18:26, Martijn Wieling wrote:
Dear useRs,
I have used using the excellent mgcv package (version 1.7-12) to
create a generalized additive model (gam) including random effects -
represented with s(...,bs="re") - on the basis of dialect data.
My model contains two random-effect factors (Word and Key - the latter
representing a speaker) and I have added both random intercepts and
various random slopes for these random-effect factors. There is no
missing data in my dataset. When I try to extract the by-word random
intercepts from my model, using coef(model), I find 357 values, equal
to the number of words in my dataset. Using coef(model) I get
uninformative names: s(Word,1) until s(Word,357), but I'm assuming (I
might be wrong though?) that I can link the labels of the words to
these values by obtaining the 357 labels from the original dataset:
unique(dat[,c("Word")])
Unfortunately, I cannot use this procedure to label the by-word random
slopes, because I find a varying number of values for these (ranging
from 346 to 356) which is always less than 357. (The number of
by-speaker random slopes does equal the number of speakers, though.)
Does anybody i) have an idea why I obtain fewer by-word random slopes
than words, and/or ii) how I can link the random slopes which are
present to the correct labels of the words?
(I did not include the model as it is>300 MB in size, but let me know
if this is necessary.)
Any help would be greatly appreciated!
With kind regards,
Martijn Wieling
University of Groningen
http://www.martijnwieling.nl
______________________________________________
R-help@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.
--
Simon Wood, Mathematical Science, University of Bath BA2 7AY UK
+44 (0)1225 386603 http://people.bath.ac.uk/sw283
______________________________________________
R-help@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.
