[R-sig-ME] glmm with a tweedie distribution
Zhang,Yanwei
Yanwei.Zhang at cna.com
Mon Mar 19 23:16:35 CET 2012
One problem with the "glmmPQL" is that the variance function can not be estimated - you need to pre-specify it in an ad hoc way. The "cpglmm" function in the "cplm" package estimates it directly from the data along with other parameters using MLE. But of course, you can use glmmPQL to generate starting values that are fed to cpglmm.
Regards,
Wayne
-----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 Bill Pikounis
Sent: Monday, March 19, 2012 4:56 PM
To: r-sig-mixed-models at r-project.org
Subject: Re: [R-sig-ME] glmm with a tweedie distribution
Bryan,
You might also wish to try the glmmPQL function in Venables and
Ripley's MASS package. Someone reported success with it on the SIG-ECO
R list nearly a year ago:
https://stat.ethz.ch/pipermail/r-sig-ecology/2011-May/002158.html
Hope that helps.
Bill
On Mon, Mar 19, 2012 at 17:18, Ben Bolker <bbolker at gmail.com> wrote:
> Danson, Bryan <Bryan.Danson at ...> writes:
>
>>
>> Is there a way to run a GLMM with a tweedie distribution?
>>
>
> Yes, using the 'cplm' package. I expanded the section on
> ZIGLMMs on http://glmm.wikidot.com/faq to include a bit more
> stuff on ZIGLMMs, including a bit on ZIGLMMs with a continuous
> response for the non-zero data.
>
> Ben
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
_______________________________________________
R-sig-mixed-models at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
NOTICE: This e-mail message, including any attachments and appended messages, is for the sole use of the intended recipients and may contain confidential and legally privileged information.
If you are not the intended recipient, any review, dissemination, distribution, copying, storage or other use of all or any portion of this message is strictly prohibited.
If you received this message in error, please immediately notify the sender by reply e-mail and delete this message in its entirety.
More information about the R-sig-mixed-models
mailing list