Set disperssion partameter to 1 in gamma regression
132 views
Skip to first unread message
Gabriela Juri
unread,
Mar 24, 2021, 8:24:23 PM3/24/21
Reply to author
Sign in to reply to author
Forward
Sign in to forward
Delete
You do not have permission to delete messages in this group
Copy link
Report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
to TMB Users
Hi!!
Hello to all.
I have time to an event data as a response variable and two categorical predictors. After intense googling and reading, I got to the conclusion that I should use an exponential distribution. This apparently means a gamma distribution with a dispersion parameter set to 1. I used Gamma(link=log) distribution as it comes with glmmtmb, checked assumptions using dharma package and got the most awful plots in the world... I think the problem is the dispersion parameter... I've read that I can correct the summary but I can't compare mean values in this way.
IS there any way to change this dispersion parameter inside the glmmtmb sentence (or other function78package)? is the exponential distribution supported?
I hope someone can point me in the right direction because I'm feeling very lost..
GaBI
Ben Bolker
unread,
Mar 24, 2021, 8:31:29 PM3/24/21
Reply to author
Sign in to reply to author
Forward
Sign in to forward
Delete
You do not have permission to delete messages in this group
Copy link
Report message
Show original message
Either email addresses are anonymous for this group or you need the view member email addresses permission to view the original message
The short answer is that you can use the 'map' argument to fix the
value of the dispersion parameter, something like map(betad=NA) . That
fixes the dispersion parameter to its starting value. The dispersion
parameter is modeled on the log scale and I believe its default start
value is so (so dispersion=exp(0)=1, which is what you want).
If you need more help or that doesn't work for some reason (I haven't
tested it!), please follow up on r-sig-mix...@r-project.org ...