[R-sig-ME] Extract/Calculate score vector from "mer" objects
Krishnaraj Ajay Gharpure
Krishnaraj.Gharpure at mu-sigma.com
Wed Dec 7 11:42:49 CET 2011
Brute force method: if you wish to compute derivatives, you can declare your log-likelihood as an R function and use the grad function (found in the numDeriv package) to evaluate that at convergence.
--
Sagely, he observed: one man's playstation is another man bakery.
-----Original Message-----
From: r-sig-mixed-models-bounces at r-project.org [mailto:r-sig-mixed-models-bounces at r-project.org] On Behalf Of Dimitris Rizopoulos
Sent: Wednesday, December 07, 2011 1:26 AM
To: r-sig-mixed-models at r-project.org
Subject: [R-sig-ME] Extract/Calculate score vector from "mer" objects
Dear All,
I'm interested in calculating the score vector (first order derivative of the log-likelihood) evaluated at the final parameter values (i.e., at convergence), for a mixed effects logistic regression model, returned by glmer(). Has anybody ever done something like this?
Thanks in advance for any hints.
Best,
Dimitris
--
Dimitris Rizopoulos
Assistant Professor
Department of Biostatistics
Erasmus University Medical Center
Address: PO Box 2040, 3000 CA Rotterdam, the Netherlands
Tel: +31/(0)10/7043478
Fax: +31/(0)10/7043014
Web: http://www.erasmusmc.nl/biostatistiek/
_______________________________________________
R-sig-mixed-models at r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
This email message may contain proprietary, private and confidential information. The information transmitted is intended only for the person(s) or entities to which it is addressed. Any review, retransmission, dissemination or other use of, or taking of any action in reliance upon, this information by persons or entities other than the intended recipient is prohibited and may be illegal. If you received this in error, please contact the sender and delete the message from your system.
Mu Sigma takes all reasonable steps to ensure that its electronic communications are free from viruses. However, given Internet accessibility, the Company cannot accept liability for any virus introduced by this e-mail or any attachment and you are advised to use up-to-date virus checking software.
More information about the R-sig-mixed-models
mailing list