Question: Shape parameter of the gamma distribution of rate variation
0
12 months ago by
dieter0
dieter0 wrote:

Hello all, I have another question for the experts where I did not find any answer, yet. It's about the shape parameter of the gamma distribution of rate variation. Known as Alpha. I know what the parameter is, and I know that if when Alpha = 1 the gamma distribution matches the exponential distribution. My question is: What does the setting "exponential" for example at MrBayes mean? At MrBayes I can set Shapepr=exponential(value). And I can set Shapepr=fixed(value). What is the difference? All the best Dieter

sequence • 465 views
modified 12 months ago by natasha.sernova3.1k • written 12 months ago by dieter0

See this pdf-file for details of gamma-distribution:

This Biostar posts on MrBayes and gamma distibution:

The Effect Of The "Invgamma" Setting On Computing The Site Likelihood In Mrbayes?

and ask stackoverflow if the above ones are not enough:

https://stackoverflow.com/questions/29902770/how-to-specify-gamma-distribution-using-shape-and-rate-in-python

Thanks Natasha, But I'm sorry, the pdf only describes the gamma distribution (which I already know) - all is clear there - but it does not answr my questions. The second link leads to a discussion about the Invar Model - this has nothing to do with my question. But thanks anyway. All the best Dieter

0
12 months ago by
Kevin Blighe33k
Republic of Ireland
Kevin Blighe33k wrote:

Hi Dieter, from what I understand, the use of exponential means that the tested gamma distributions become virtually infinite, as it would then be based on an exponential range of α (alpha). I think that the value that's passed to `Exponential` may relate to λ (lambda).

The other options for `Shapepr` are:

• uniform (the default, I believe), α values between 0 and 1.0 are tested
• fixed, fixed at a specific value of α

You may wish to confirm with the author of MrBayes.

Hello Kevin, I checked it again and: You are right! It is exactly how you said. I found the proof here: http://hydrodictyon.eeb.uconn.edu/eebedia/index.php/Phylogenetics:_MrBayes_Lab It is the exponetial distribution and tha parameter is no longer alpha but lambda. The interesting thing is that the new standard setting at MrBayes is the exponetial distribution with Lambda = 1 [setting Shapepr=exponential(1)]. So thanks for all. All the best Dieter