M2a_test
Usage
M2a_test(submodel: Double
-> CTMC<Codons<a>>, omega0: Double, p0: Double, posP: Double, posW: Double, posSelection: Int) →
DiscreteDist<CTMC<Codons<a>>>
Arguments
A default beginning with ~ specifies a prior
distribution.
-
submodel: -
The model as a function of dN/dS
-
Default:
|w:GY94(omega=w)| -
omega0: -
The dN/dS value for conserved sites
-
Default:
~Uniform(0,1) -
p0: -
The fraction of conserved sites among non-positively selected sites.
-
Default:
~Uniform(0,1) -
posP: -
The fraction of positively selected sites
-
Default:
~Beta(1,10) -
posW: -
The dN/dS value for positively selected sites
-
Default:
~LogGamma(4,0.25) -
posSelection: -
The model selector: 1 if positive selection, 0 if not.
-
Default:
~Bernoulli(0.5)
Description
This model performs a Bayesian version of the test of M2a against M1a. The posterior probability that posSelection=1 is the posterior probability of positive selection.
The posterior mean of PrPosSelection can provide a more accurate estimate of this probability than the posterior mean of posSelection. The statreport tool uses LogOddsPosSelection to report the corresponding posterior log odds accurately, even when the probability is extremely close to 0 or 1. (Do not average LogOddsPosSelection directly.)
Examples
|w:GY94(omega=w,pi=F1x4)| +> M2a_test|w:MG94(omega=w)| +> M2a_test|w:FMutSel0(omega=w)| +> M2a_testCitation
Wong, Wendy S. W.; Yang, Ziheng; Goldman, Nick; Nielsen, Rasmus (2004). Accuracy and Power of Statistical Methods for Detecting Adaptive Evolution in Protein Coding Sequences and for Identifying Positively Selected Sites. Genetics 168(2): 1041--1051. DOI: 10.1534/genetics.104.031153