search for: niterpql

...he 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 added some noise to lon and lat, modelled the density instead of using a count with model offset, and we normali...

...ed Additive Mixed Model, as covariates I have the geographical longitude and latitude locations of indexed data. I include a random effect for each district (dist) so the code is fit <- gamm(y~s(lon,lat,bs="tp", m=2)+offset(log(exp.)), random=list(dist=~1), family="poisson", niterPQL=30) when I run the gamm function I obtain the next error message: %%% Maximum number of PQL iterations: 30 iteration 1 iteration 2 ... iteration 8 iteration 9 iteration 10 Error in chol(XVX + S) : the leading minor of order 29 is not positive definite %%% Could be any problem in gamm() ??? --...

...ernally to glmmPQL (see results of traceback()) [large data vector scrolls off screen] 6: eval(expr, envir, enclos) 5: eval(mcall) 4: glmmPQL(y ~ X - 1, random = rand, data = strip.offset(mf), family = family, correlation = correlation, control = control, weights = weights, niter = niterPQL, verbose = verbosePQL) 3: eval(expr, envir, enclos) 2: eval(parse(text = paste("ret$lme<-glmmPQL(", deparse(fixed.formula), ",random=rand,data=strip.offset(mf),family=family,correlation=correlation,co ntrol=control,", "weights=weights,niter=niterPQL,verbose=ve...

...station (StnF), and a spatial correlation structure, for Depth_m. When trying to fit an equivalent GAMM to start with, I used the formula: SB.gamm1<-gamm(count~offset(vol_offset)+ s(Depth_m,by=StnF,bs="cs")+StageF*RegionF, family=negative.binomial(theta=1.41), niterPQL=50, data=Zoop_2011[Zoop_2011$SpeciesF=='SB',]) The theta=1.41 came from SB.gam4, using the theta.ml() command. However, this model does not converge, even with 5000 iterations. When I increase theta to 6, the model converges, but I'm not sure why. I have also tried theta...

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