Thank you Simon,
that works well.
I had tried previously
m1 <-gamm(H ~ Tillage + s(Year, by =Tillage), random=list(Block=~ Tillage),
na.action=na.omit, data =mydata, correlation = corARMA(form =~
Year|Tillage/Block, p = 1, q = 0))
and that gave me convergence problems.
Thank you!
eva,
The problem is that the random effects and correlation structure that
you have specified don't meet the restrictions required by lme (which
gamm calls). You can see this by trying
M0 <-lme(H ~ Tillage , random=list(Block=~1), na.action=na.omit, data =
mydata, correlation = corARMA(form
eva,
I can't manage to replicate this by simulation. Is there any chance you
could send me the data off-list, and I can take a look (if so I'll only
use the data for the purpose of finding out what's wrong, of course).
best,
Simon
On 17/02/12 09:35, HERNANDEZ PLAZA, MARIA EVA wrote:
HELLO
HELLO ALL,
I AM GETTING AN ERROR MESSAGE WHEN TRYING TO RUN A GAMM MODEL LIKE THE ONE
BELOW.
I AM USING R VERSION 2.14.1 (2011-12-22) AND MGCV 1.7-12.
M1 <-gamm(DepVar ~ Treatment + s(Year, by =Treatment), random=list(Block=~1),
na.action=na.omit, data = mydata, correlation = corARMA(f
4 matches
Mail list logo
