This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Shape parameter of the gamma distribution of rate variation

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
• 2,984 views • link
updated by natasha.sernova natasha.sernova • written by dieter • dieter

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 views • link
written by dieter • dieter

Hello dieter!

We believe that this post does not fit the main topic of this site.

Not strictly bioinformatics, more appropriate at https://stats.stackexchange.com/

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

• 0 views • link
written by WouterDeCoster WouterDeCoster

@Wouter: OP is asking about MrBayes package and the options in it. So it is a valid bioinformatics question.

• 0 views • link
written by GenoMax GenoMax

1 answer

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.

• 0 views • link
written by Kevin Blighe Kevin Blighe

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

• 0 views • link
written by dieter • dieter

Log in to answer this question.

More posts like this