Displaying 4 results from an estimated 4 matches for "gamm2".
Did you mean:
gamma
2008 Nov 19
2
GAMM and anove.lme question
...) in 1,000
meter long intervals (elementary sampling units) from multiple transects
(each 20-30 km long, tranf) in two different lakes and three different
years.
bloat.gamm1 <- gamm(lbloat ~ s(depth), correlation=corSpher(c(30000,
0.01),form = ~ x+y|tranf, nugget=TRUE), data=fish3)
bloat.gamm2 <- gamm(lbloat ~ lakef + s(depth),
correlation=corSpher(c(30000, 0.01),form = ~ x+y|tranf, nugget=TRUE),
data=fish3)
bloat.gamm3 <- gamm(lbloat ~ lakef + s(depth, by=lakef),
correlation=corSpher(c(30000, 0.01),form = ~ x+y|tranf, nugget=TRUE),
data=fish3)
> anova.lme(bloat.gamm1$lm...
2011 Oct 05
2
gamm: problems with corCAR1()
...when analyzing the data of one replicate, I received following
error:
/Error in MEestimate(lmeSt, grps) :
Singularity in backsolve at level 0, block 1
/
I have no idea how to solve this. It seems to be related with the complexity
of the model because no error occurred when running a simpler gamm (gamm2):
gamm2 <- gamm(richness~
s(Week,k=6,fx=FALSE) + factor(Treatment),
correlation=corCAR1(form=~Week|conc.f), data=test,family=gaussian)
Any help would be well appreciated!
With kind regards,
Karel
--
View this message in context: http://r.789695.n4.nabble.com/gamm-problems-with-corCAR1-t...
2009 Jan 22
1
convergence problem gamm / lme
...(arg 1)
We tried several things. We added some
noise to lon and lat, modelled the density instead of using a count with model
offset, and we normalized the explanatory variables. We also changed several
settings (see models below).
Interestingly, we do manage to fit a more
complex model:
gamm2<-gamm(count~offset(offsetter)+
s(lat,lon,year,dayofyear), random=list(code_tripnr=~1),family="poisson",
correlation = corGaus(0.1, form=~lat + lon))
The models are fitted using mgcv 1.4-1 and
R 2.7.1 on a 64Bits Debian OS.
So there seems to be a convergence problem, correct? A...
2009 Jan 28
2
t.test in a loop
Hi All,
I've been having a little trouble with creating a loop that will run a a
series of t.tests for inspection,
Below is the code i've tried, and some checks i've looked at.
I've used the get(paste()) idea as i was told previously that the use of the
eval should try and be avoided.
I've run a single syntax to check that my systax is correct and works
without any problems