[R] convergence problem gamm / lme
Simon Wood
s.wood at bath.ac.uk
Thu Jan 29 16:34:41 CET 2009
Geert,
Sorry for slow reply... I don't see any obvious problems with what you've
done, so I guess it's the usual problem that PQL just doesn't *have* to
converge, and the bit of extra flexibility of using a smooth is too much for
it in this case. If you send me the data offline I can dig a little bit more
if you like (I'll only use the data for this purpose etc. etc.)
You are right that PQL does the same thing for Poisson and quasi-poisson. I
don't think there is an easy way to use the values for the reduced dataset
fit in the full dataset fitting, unfortunately.
Another option is to use `gam' to fit the random effects. It'll be a bit slow
with 70+ random effects, as you have, and it's a bit more work to set up, but
it should converge. See ?gam.models which has some examples showing how to do
this.
best,
Simon
On Thursday 29 January 2009 08:20, geert aarts wrote:
> Simon, thanks for your reply and your suggestions.
>
> I fitted the following glmm's
>
> gamm3<-try(glmmPQL(count~offset(offsetter)+poly(lon,3)*poly(lat,3),random=l
>ist(code_tripnr=~1),family="poisson"))
>
> Which worked OK (see summary below)
>
> I also fitted a model using quasipoisson, but that didn't help. I actually
> also thought that glmmPQL and gamm estimate the dispersion parameter and
> hence assumes a quasipoisson distribution, even if you specify poisson. Is
> that correct?
>
> Finally I tried fitting a model to less data, and sometimes gamm managed to
> converge (see summary below). So would it be possible to use the parameter
> estimates from the model fitted to less data as starting values for the
> gamm fitted to the full data set? Or do you have any other suggestions?
>
> Thanks.
> Cheers Geert
>
>
>
>
>
>
> gamm3<-try(glmmPQL(count~offset(offsetter)+poly(lon,3)*poly(lat,3),random=l
>ist(code_tripnr=~1),f
>
> amily="poisson"))
>
>
>
> iteration
> 1
>
> iteration
> 2
>
> iteration
> 3
>
> > detach(Disc_age)
>
> summary(gamm3)
>
> Linear
> mixed-effects model fit by maximum likelihood
>
> Data: NULL
>
> AIC BIC logLik
>
> NA NA
> NA
>
>
>
> Random
> effects:
>
> Formula: ~1 | code_tripnr
>
> (Intercept) Residual
>
> StdDev:
> 0.001391914 231.9744
>
>
>
> Variance
> function:
>
> Structure: fixed weights
>
> Formula: ~invwt
>
> Fixed
> effects: count ~ offset(offsetter) + poly(lon, 3) * poly(lat, 3)
>
> Value
> Std.Error DF t-value p-value
>
> (Intercept) -1.582 11.96 2024 -0.13232174 0.8947
>
> poly(lon,
> 3)1 -4.048 1397.33 2024 -0.00289673 0.9977
>
> poly(lon,
> 3)2 -22.013 699.71 2024 -0.03145996 0.9749
>
> poly(lon,
> 3)3 -8.538 593.87 2024 -0.01437683 0.9885
>
> poly(lat,
> 3)1 -109.624 666.05 2024 -0.16458856 0.8693
>
> poly(lat,
> 3)2 -104.179 381.37 2024 -0.27316977 0.7848
>
> poly(lat,
> 3)3 -10.661 221.93 2024 -0.04803585 0.9617
>
> poly(lon,
> 3)1:poly(lat, 3)1 4290.737 61369.98 2024
> 0.06991589 0.9443
>
> poly(lon,
> 3)2:poly(lat, 3)1 1853.559 36835.63 2024
> 0.05031972 0.9599
>
> poly(lon,
> 3)3:poly(lat, 3)1 -240.521 25771.80 2024 -0.00933272 0.9926
>
> poly(lon,
> 3)1:poly(lat, 3)2 2540.147 41378.38 2024
> 0.06138826 0.9511
>
> poly(lon,
> 3)1:poly(lat, 3)2 2540.147 41378.38 2024
> 0.06138826 0.9511
>
> poly(lon,
> 3)2:poly(lat, 3)2 -1803.911 21522.17
> 2024 -0.08381643 0.9332
>
> poly(lon,
> 3)3:poly(lat, 3)2 1040.858 16352.56 2024
> 0.06365109 0.9493
>
> poly(lon,
> 3)1:poly(lat, 3)3 632.587 12180.28 2024
> 0.05193535 0.9586
>
> poly(lon,
> 3)2:poly(lat, 3)3 -394.339 13088.72 2024 -0.03012818 0.9760
>
> poly(lon,
> 3)3:poly(lat, 3)3 -543.502 6221.71 2024 -0.08735569 0.9304
>
> Correlation:
>
> (Intr) ply(ln,3)1
> ply(ln,3)2 ply(ln,3)3 ply(lt,3)1
>
> poly(lon,
> 3)1 0.889
>
> poly(lon,
> 3)2 0.938 0.878
>
> poly(lon,
> 3)3 0.843 0.981
> 0.792
>
> poly(lat,
> 3)1 -0.829 -0.949 -0.906
> -0.882
>
> poly(lat,
> 3)2 0.859 0.578 0.742
> 0.538 -0.474
>
> poly(lat,
> 3)3 -0.552 -0.783 -0.579
> -0.756 0.837
>
> poly(lon,
> 3)1:poly(lat, 3)1 -0.947 -0.974
> -0.940 -0.940 0.925
>
> poly(lon,
> 3)2:poly(lat, 3)1 -0.934 -0.950
> -0.857 -0.929 0.881
>
> poly(lon,
> 3)3:poly(lat, 3)1 -0.818 -0.963
> -0.866 -0.945 0.931
>
> poly(lon,
> 3)1:poly(lat, 3)2 0.808 0.975
> 0.784 0.968 -0.928
>
> poly(lon,
> 3)2:poly(lat, 3)2 0.737 0.575
> 0.853 0.465 -0.659
>
> poly(lon,
> 3)3:poly(lat, 3)2 0.735 0.896
> 0.647 0.938 -0.765
>
> poly(lon,
> 3)1:poly(lat, 3)3 -0.794 -0.592
> -0.823 -0.518 0.591
>
> poly(lon,
> 3)2:poly(lat, 3)3 -0.542 -0.737
> -0.419 -0.781 0.635
>
> poly(lon,
> 3)3:poly(lat, 3)3 -0.398 -0.383
> -0.534 -0.334 0.425
>
> ply(lt,3)2
> ply(lt,3)3 p(,3)1:(,3)1 p(,3)2:(,3)1
>
> poly(lon,
> 3)1
>
> poly(lon,
> 3)2
>
> poly(lon,
> 3)3
>
> poly(lat,
> 3)1
>
> poly(lat,
> 3)2
>
> poly(lat,
> 3)3 -0.136
>
> poly(lon,
> 3)1:poly(lat, 3)1 -0.708 0.690
>
> poly(lon,
> 3)2:poly(lat, 3)1 -0.701 0.710 0.933
>
> poly(lon,
> 3)3:poly(lat, 3)1 -0.499 0.738 0.956 0.849
>
> poly(lon,
> 3)1:poly(lat, 3)2 0.458 -0.845
> -0.915 -0.934
>
> poly(lon,
> 3)2:poly(lat, 3)2 0.683 -0.344
> -0.719 -0.522
>
> poly(lon,
> 3)2:poly(lat, 3)2 0.683 -0.344
> -0.719 -0.522
>
> poly(lon,
> 3)3:poly(lat, 3)2 0.464 -0.655
> -0.834 -0.884
>
> poly(lon,
> 3)1:poly(lat, 3)3 -0.823 0.241 0.752 0.594
>
> poly(lon,
> 3)2:poly(lat, 3)3 -0.300 0.707 0.612 0.788
>
> poly(lon,
> 3)3:poly(lat, 3)3 -0.266 0.148 0.493 0.250
>
> p(,3)3:(,3)1
> p(,3)1:(,3)2 p(,3)2:(,3)2 p(,3)3:(,3)2
>
> poly(lon,
> 3)1
>
> poly(lon,
> 3)2
>
> poly(lon,
> 3)3
>
> poly(lat,
> 3)1
>
> poly(lat,
> 3)2
>
> poly(lat,
> 3)3
>
> poly(lon,
> 3)1:poly(lat, 3)1
>
> poly(lon,
> 3)2:poly(lat, 3)1
>
> poly(lon,
> 3)3:poly(lat, 3)1
>
> poly(lon,
> 3)1:poly(lat, 3)2 -0.928
>
> poly(lon,
> 3)2:poly(lat, 3)2 -0.637 0.432
>
> poly(lon,
> 3)3:poly(lat, 3)2 -0.851
> 0.935 0.245
>
> poly(lon,
> 3)1:poly(lat, 3)3 0.642 -0.482 -0.894 -0.410
>
> poly(lon,
> 3)2:poly(lat, 3)3 0.609 -0.822 0.007 -0.847
>
> poly(lon,
> 3)3:poly(lat, 3)3 0.551 -0.327 -0.637 -0.291
>
> p(,3)1:(,3)3
> p(,3)2:(,3)3
>
> poly(lon,
> 3)1
>
> poly(lon,
> 3)2
>
> poly(lon,
> 3)3
>
> poly(lat,
> 3)1
>
> poly(lat,
> 3)2
>
> poly(lat,
> 3)3
>
> poly(lon,
> 3)1:poly(lat, 3)1
>
> poly(lon,
> 3)2:poly(lat, 3)1
>
> poly(lon,
> 3)3:poly(lat, 3)1
>
> poly(lon,
> 3)1:poly(lat, 3)2
>
> poly(lon,
> 3)2:poly(lat, 3)2
>
> poly(lon,
> 3)3:poly(lat, 3)2
>
> poly(lon,
> 3)1:poly(lat, 3)3
>
> poly(lon,
> 3)3:poly(lat, 3)1
>
> poly(lon,
> 3)1:poly(lat, 3)2
>
> poly(lon,
> 3)2:poly(lat, 3)2
>
> poly(lon,
> 3)3:poly(lat, 3)2
>
> poly(lon,
> 3)1:poly(lat, 3)3
>
> poly(lon,
> 3)2:poly(lat, 3)3 0.080
>
> poly(lon,
> 3)3:poly(lat, 3)3 0.684 -0.180
>
>
>
> Standardized
> Within-Group Residuals:
>
> Min Q1 Med Q3 Max
>
> -0.504980771 -0.000866948
> 0.028470924 0.078583094
> 33.247831244
>
>
>
> Number
> of Observations: 2113
>
> Number
> of Groups: 74
>
>
>
>
>
>
>
> gamm3<-try(gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=
>~1),family="quasipoisson", niterPQL=200))
>
>
>
>
> summary(gamm3$gam)
>
>
>
> Family:
> quasipoisson
>
> Link
> function: log
>
>
>
> Formula:
>
> count
> ~ offset(offsetter) + s(lon, lat)
>
>
>
> Parametric
> coefficients:
>
> Estimate Std. Error t value Pr(>|t|)
>
> X 1.31370
> 0.09854 13.33
>
>
>
>
> summary(gamm3$lme)
>
> Linear
> mixed-effects model fit by maximum likelihood
>
> Data: data
>
> AIC
> BIC logLik
>
> 2808.398 2837.845 -1398.199
>
>
>
> Random
> effects:
>
> Formula: ~Xr.1 - 1 | g.1
>
> Structure: pdIdnot
>
> Xr.11 Xr.12
> Xr.13 Xr.14 Xr.15
> Xr.16 Xr.17 Xr.18
>
> StdDev:
> 12.49623 12.49623 12.49623 12.49623 12.49623 12.49623 12.49623 12.49623
>
> Xr.19 Xr.110
> Xr.111 Xr.112 Xr.113
> Xr.114 Xr.115 Xr.116
>
> StdDev:
> 12.49623 12.49623 12.49623 12.49623 12.49623 12.49623 12.49623 12.49623
>
> Xr.117 Xr.118
> Xr.119 Xr.120 Xr.121
> Xr.122 Xr.123 Xr.124
>
> StdDev:
> 12.49623 12.49623 12.49623 12.49623 12.49623 12.49623 12.49623 12.49623
>
> Xr.125 Xr.126
> Xr.127
>
> StdDev:
> 12.49623 12.49623 12.49623
>
>
>
> Formula: ~1 | code_tripnr %in% g.1
>
> (Intercept) Residual
>
> StdDev: 0.8132693 5.077804
>
>
>
> Variance
> function:
>
> Structure: fixed weights
>
> Formula: ~invwt
>
> Fixed
> effects: list(fixed)
>
> Value Std.Error
> DF t-value p-value
>
> XX 1.3137042 0.09863463 923
> 13.318894 0.0000
>
> Xs(lon,lat)Fx1
> -0.4406352 0.23114503 923 -1.906315
> 0.0569
>
> Xs(lon,lat)Fx2
> -0.6217519 0.24918031 923 -2.495189 0.0128
>
> Correlation:
>
> XX X(,)F1
>
> Xs(lon,lat)Fx1 0.015
>
> Xs(lon,lat)Fx2
> -0.009 -0.148
>
>
>
> Standardized
> Within-Group Residuals:
>
> Min Q1 Med Q3 Max
>
> -3.42951750 -0.37448354
> 0.06432438 0.53690322 8.62026552
>
>
>
> Number
> of Observations: 1000
>
> Number
> of Groups:
>
> g.1 code_tripnr %in% g.1
>
> 1 75
>
>
>
>
>
>
>
> ----------------------------------------
>
> > From: s.wood at bath.ac.uk
> > To: r-help at r-project.org
> > Date: Fri, 23 Jan 2009 11:32:21 +0000
> > Subject: Re: [R] convergence problem gamm / lme
> >
> > Geert,
> >
> > Can you get a simpler model with, say, a quadratic dependence on lon, lat
> > to converge, using glmmPQL? The answer might give a clue about whether
> > the issue is related to using a smoother, or is something more basic.
> >
> > How confident are you that the Poisson assumption is reasonable?
> >
> > Can the model be fitted to a random subsample of the data, or does it
> > always fail? PQL can fail to converge, but it's usually not as obstinate
> > as it seems to be in this case, if the model structure is reasonable and
> > identifiable.
> >
> > best,
> > Simon
> >
> > On Thursday 22 January 2009 15:52, geert aarts wrote:
> >> Hope one of you could help with the following question/problem:
> >> We would like to explain the spatial
> >> distribution of juvenile fish. 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 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? And does someone
> >> have an idea what might cause this? Secondly are there some
> >> tricks/solutions. E.g. perhaps we could use the results from the more
> >> complex model (gamm2 above), but I do not know exactly how. All
> >> help/advice would be greatly appreciated.
> >>
> >>
> >>
> >> Kind regards, Geert
> >>
> >>
> >>
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),
> >> random=list(code_tripnr=~1),family="poisson", correlation = corExp(1,
> >> form=~X + Y),nite
> >>
> >> rPQL=200)
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in recalc.corSpatial(object[[i]],
> >> conLin) :
> >>
> >>
> >> NA/NaN/Inf in foreign function call (arg 1)
> >>
> >>> gamm3<-gamm(count~offset(offsetter)+s(lon,lat,k=c(1,1)),random=list(cod
> >>>e_ tripnr=~1),family="poisson",
> >>
> >> niterPQL=200)
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in lme.formula(fixed = fixed, random
> >> = random, data = data, correlation = correlation, :
> >>
> >> nlminb
> >> problem, convergence error code = 1
> >>
> >>
> >> message = false convergence (8)
> >>
> >> In addition: Warning messages:
> >>
> >> 1: In if (k < M + 1) { :
> >>
> >> the
> >> condition has length> 1 and only the first element will be used
> >>
> >>
> >>
> >>
> >>
> >> .Options$mgcv.vc.logrange=0.001 # we also
> >> tried higher settings
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~
> >>1), family="poisson", niterPQL=200, control=lmeControl(opt="optim"))
> >>
> >>
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in optim(c(coef(lmeSt)),
> >> function(lmePars) -logLik(lmeSt, lmePars),
> >>
> >>
> >>
> >> initial value in 'vmmin' is not finite
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~
> >>1), family="poisson", niterPQL=200,control=lmeControl(minAbsParApV
> >>
> >> ar=0.0000000000001))
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in recalc.corSpatial(object[[i]],
> >> conLin) :
> >>
> >>
> >> NA/NaN/Inf in foreign function call (arg 1)
> >>
> >>
> >>
> >>
> >> 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)
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat,k=c(1,1)),random=list(code
> >>_tr ipnr=~1),family="poisson", niterPQL=200)
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in lme.formula(fixed = fixed, random
> >> = random, data = data, correlation = correlation, :
> >>
> >>
> >> nlminb problem, convergence
> >> error code = 1
> >>
> >>
> >> message = false convergence (8)
> >>
> >> In addition: Warning messages:
> >>
> >> 1: In if (k < M + 1) { :
> >>
> >> the
> >> condition has length> 1 and only the first element will be used
> >>
> >> 2: In smooth.construct.tp.smooth.spec(object,
> >> dk$data, dk$knots) :
> >>
> >>
> >> basis dimension, k, increased to minimum possible
> >>
> >>
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat,k=c(8,8)),random=list(code
> >>_tr ipnr=~1),family="poisson", niterPQL=200)
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in lme.formula(fixed = fixed, random
> >> = random, data = data, correlation = correlation, :
> >>
> >>
> >> nlminb problem, convergence
> >> error code = 1
> >>
> >>
> >> message = false convergence (8)
> >>
> >> In addition: Warning messages:
> >>
> >> 1: In if (k < M + 1) { :
> >>
> >> the
> >> condition has length> 1 and only the first element will be used
> >>
> >> 2: In 1:UZ.len : numerical expression has 2
> >> elements: only the first used
> >>
> >> 3: In if (p.rank> ncol(XZ)) p.rank
> >> <- ncol(XZ) :
> >>
> >> the
> >> condition has length> 1 and only the first element will be used
> >>
> >> 4: In 1:p.rank : numerical expression has 2
> >> elements: only the first used
> >>
> >> 5: In if (p.rank < k - j) Xf <- XZU[,
> >> (p.rank + 1):(k - j), drop = FALSE] else Xf <- matrix(0, :
> >>
> >> the
> >> condition has length> 1 and only the first element will be used
> >>
> >> 6: In (p.rank + 1):(k - j) :
> >>
> >>
> >> numerical expression has 2 elements: only the first used
> >>
> >> 7: In 1:p.rank : numerical expression has 2
> >> elements: only the first used
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat,k=c(4,4),fx=T),random=list
> >>(co de_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)
> >>
> >> In addition: Warning messages:
> >>
> >> 1: In if (k < M + 1) { :
> >>
> >> the
> >> condition has length> 1 and only the first element will be used
> >>
> >> 2: In 1:UZ.len : numerical expression has 2
> >> elements: only the first used
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+te(lon,lat),random=list(code_tripnr=
> >>~1) ,family="poisson", niterPQL=200,control=lmeControl(opt="opti
> >>
> >> m"))
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in optim(c(coef(lmeSt)),
> >> function(lmePars) -logLik(lmeSt, lmePars),
> >>
> >>
> >>
> >> initial value in 'vmmin' is not finite
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~
> >>1), family="poisson", niterPQL=200,control=lmeControl(tolerance=
> >>
> >> 0.00000000000001))
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in MEestimate(lmeSt, grps) :
> >>
> >>
> >> NA/NaN/Inf in foreign function call (arg 1)
> >>
> >>> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=
> >>>~1 ),family="poisson",
> >>
> >> niterPQL=200,control=lmeControl(niterEM=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)
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~
> >>1), family="poisson",
> >> niterPQL=200,control=lmeControl(msTol=0.00000000000000001))
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in MEestimate(lmeSt, grps) :
> >>
> >>
> >> NA/NaN/Inf in foreign function call (arg 1)
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~
> >>1), family="poisson",
> >> niterPQL=200,control=lmeControl(.relStep=0.00000000000000000001))
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in MEestimate(lmeSt, grps) :
> >>
> >>
> >> NA/NaN/Inf in foreign function call (arg 1)
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~
> >>1), family="poisson",
> >> niterPQL=200,control=lmeControl(nlmStepMax=0.00000000000000000001))
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in MEestimate(lmeSt, grps) :
> >>
> >>
> >> NA/NaN/Inf in foreign function call (arg 1)
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~
> >>1), family="poisson",
> >> niterPQL=200,control=lmeControl(minAbsParApVar=0.0000000000001))
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in MEestimate(lmeSt, grps) :
> >>
> >>
> >> NA/NaN/Inf in foreign function call (arg 1)
> >>
> >>
> >>
> >>
> >> gamm3<-gamm(count~offset(offsetter)+s(lon,lat),random=list(code_tripnr=~
> >>1), family="poisson", niterPQL=200, control=lmeControl(returnObject=T))
> >>
> >> Maximum number of PQL iterations: 200
> >>
> >> iteration 1
> >>
> >> iteration 2
> >>
> >> Error in MEestimate(lmeSt, grps) :
> >>
> >>
> >> Singularity in backsolve at level 0, block 1
> >>
> >> In addition: Warning messages:
> >>
> >> 1: In logLik.reStruct(object, conLin) :
> >>
> >>
> >> Singular precision matrix in level -1, block 1
> >>
> >> 2: In logLik.reStruct(object, conLin) :
> >>
> >>
> >> Singular precision matrix in level -1, block 1
> >>
> >> 3: In logLik.reStruct(object, conLin) :
> >>
> >>
> >> Singular precision matrix in level -1, block 1
> >>
> >> 4: In logLik.reStruct(object, conLin) :
> >>
> >>
> >> Singular precision matrix in level -1, block 1
> >>
> >> 5: In logLik.reStruct(object, conLin) :
> >>
> >>
> >> Singular precision matrix in level -1, block 1
> >>
> >> 6: In MEestimate(lmeSt, grps) :
> >>
> >>
> >> Singular precision matrix in level -1, block 1
> >>
> >>
> >> _________________________________________________________________
> >>
> >>
> >> [[alternative HTML version deleted]]
> >
> > --
> >
> >> Simon Wood, Mathematical Sciences, University of Bath, Bath, BA2 7AY UK
> >> +44 1225 386603 www.maths.bath.ac.uk/~sw283
>
> _________________________________________________________________
> De leukste online filmpjes vind je op MSN Video!
> http://video.msn.com/video.aspx?mkt=nl-nl
> ______________________________________________
> R-help at 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 Sciences, University of Bath, Bath, BA2 7AY UK
> +44 1225 386603 www.maths.bath.ac.uk/~sw283
More information about the R-help
mailing list