[R-sig-ME] Avoid errors in pwrssUpdate ?

Pierre Morel pier.morel at gmail.com
Wed Jun 13 13:09:13 CEST 2012


OK, thanks for the suggestion of glmmADMB, I will look that up.
I also found out that I could get things working with nAGQ=0.
However the doc says it's less exact, and profile() gives a "pp$setTheta(theta) : theta size mismatch" error in this case.

Pierre Morel


Le 13 juin 2012 à 10:54, Ben Bolker a écrit :

> Pierre Morel <pier.morel at ...> writes:
> 
>> I am getting a lot of pwrssUpdate errors when trying to model my
>> data with gmler (using the most recent version from svn... I don't
>> know if the previous versions were affected).
> 
>> These errors are "PIRLS step failed" or "pwrssUpdate did not
>> converge in 30 iterations". I understand that these means that the
>> algorithm does not manage to work with my data, but there is a
>> peculiar behavior, and my data doesn't seem too unreasonable to fit
>> with the model I want to use, so I am wondering if the problem is on
>> my side !
> 
> [snip to make gmane happier]
> 
>> Here is the model I want to fit, which doesn't seem unreasonable
>> given the figure (random slopes and intercepts for subjects) :
> 
>> model<-glmer(cbind(RuleReach,NTrials-RuleReach)~RuleWeight+
>> (RuleWeight|Subject),data=rewardalldirsub,family=binomial)
> 
>> However this gives me the "pwrssUpdate did not converge in 30
>> iterations" error.  What is surprinsing, is that if I do not use the
>> rightmost points (RuleWeight of 1), the model converges, even though
>> there are less datapoints and the remaining points are the noisiest
>> (subjects follow the rule quite reliably when it has a weight of 1
>> as you can see).
> 
>> Removing the correlation in the random effects works sometimes (but
>> not on all my sub data sets), and having a random intercept only
>> (which is obviously not correct) is the only thing that seems to
>> work in all cases.
> 
>> Centering RuleWeight (ie having it between -1 and 1 instead of 0 and
>> 1) doesn't work.  Any ideas on why this doesn't work / how to make
>> it work ?
> 
>  Thanks for the report: this is an issue the developers are
> (painfully) aware of, and working on.  The issue arises mostly when
> the predictions for some observations are very close to 0 or 1 (which
> explains why using the rightmost points helps ...)  You have tried all
> the obvious things I know of.  I would additionally try (1) setting
> starting values by hand and/or (2) trying out glmmADMB ...
> 
>  Ben Bolker
> 
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models



More information about the R-sig-mixed-models mailing list