search for: gamm3

...ish. We have 2135 records, from 75 vessels (code_tripnr) and 7 to 39 observations for each vessel, hence the random effect for code_tripnr. The offset (‘offsetter’) accounts for the haul duration and sub sampling factor. There are no extreme outliers in lat/lon. The model we try to fit is: > gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~1),family="poisson", niterPQL=200) Maximum number of PQL iterations: 200 iteration 1 iteration 2 Error in MEestimate(lmeSt, grps) : NA/NaN/Inf in foreign function call (arg 1) We tried several things. We add...

...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$lme, bloat.gamm2$lme, bloat.gamm3$lme) Model df AIC BIC logLik Test L.Ratio p-value bloat.gamm1$lme...

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