[R-sig-ME] Log likelihood of a glmer() binomial model .

[Phillip Alday]

And speaking of things not getting through, here's the my message again
without digital signing.

Phillip

On 25/4/19 11:26 am, Phillip Alday wrote:
The X.txt got through, but not the demo.R. Care to resend the demo.R as
demo.X and/or put it up as a gist/code snippet on GitHub/GitLab/etc.?

Best,
Phillip

> 
> On 25/4/19 11:23 am, Rolf Turner wrote:
>>
>> On 20/04/19 12:44 PM, Ben Bolker wrote:
>>
>> <SNIP>
>>
>>>   Here is an **inefficient** method for computing the likelihood
>>>
>>>     coefs <- unlist(getME(fit,c("theta","beta"))
>>>     newdev <- update(fit, data=VS, devFunOnly=TRUE)
>>>     newdev(coefs)
>>
>> <SNIP>
>>
>> Well I (finally) tried that, and got an error.  Since I haven't a clue
>> as to what is going on, there is no way that I can do any debugging.
>>
>> I'm hoping that Ben, or someone, can tell me what I doing rong ( :-) ).
>>
>> I have attached a sample data set (simulated; in the file X.txt) and the
>> code that I used in the file demo.R.  I *think* that *.txt and *.R files
>> pass through the list server OK.  This should demonstrate the error that
>> arises.
>>
>> If you source demo.R (having made sure that X.txt is in the working
>> directory) you will the error message:
>>
>>> Error in pp$setTheta(theta) : theta size mismatch
>>
>> The crucial code is essentially copy-and-pasted from Ben's email
>> (referred to above) so I have no idea where I'm messing up.
>>
>> Can anyone enlighten me?  Thanks.
>>
>> cheers,
>>
>> Rolf
>>
>>
>> _______________________________________________
>> R-sig-mixed-models using r-project.org mailing list
>> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
>>
>