M8a — Beta-distributed omega with an additional
neutral category
Usage
M8a(submodel: Double
-> CTMC<Codons<a>>, n: Int, mu: Double, v: Double, posP: Double) → DiscreteDist<CTMC<Codons<a>>>
Arguments
A default beginning with ~ specifies a prior
distribution.
-
submodel: -
Codon model as a function of omega
-
Default:
|w:GTR +> MNM +> MutSelAA +> dNdS(omega=w)| -
n: -
Number of categories approximating the Beta distribution
-
Default:
10 -
mu: -
The mean of the Beta distribution
-
Default:
~Uniform(0,1) -
v: -
The normalized Beta variance Var(omega)/(mu*(1-mu)); 0 < v < 1
-
Default:
~Uniform(0,1) -
posP: -
The fraction of neutral sites
-
Default:
~Beta(1,10)
Description
M8a combines an M7 component of weight 1-posP with an additional neutral category of weight posP and omega=1. M8a is equivalent to M8(posW=1) and serves as the null model in M8a_test.
The default background uses GTR, MNM, and MutSelAA. Here omega scales nonsynonymous rates after accounting for amino-acid fitnesses. Thus omega=1 retains these fitness effects; neutrality and positive selection are assessed relative to this background.
The Beta component uses Gauss-Jacobi quadrature with generally unequal category weights. There are n+1 categories in total.
Citation
Swanson, Willie J.; Nielsen, Rasmus; Yang, Qiaofeng (2003). Pervasive Adaptive Evolution in Mammalian Fertilization Proteins. Molecular Biology and Evolution 20(1): 18--20. DOI: 10.1093/oxfordjournals.molbev.a004233