Function to plot a distribution of fitnesses based on codon equilibrium freqs
PlotExpectedFitness(codon.fitnesses.matrices, codon.eq.matrices, values,
optimal.aa = NULL, palette = "Set1", lwd = 2,
include.stop.codon = FALSE, type = "histogram", fitness = TRUE,
numcode = 1, ...)
A 3d array of aa.fitness.matrix returned from ComputeEquilibriumAAFitness (first element in return)
A 3d array of codon equilibrium frequencies
The vector of labels for each matrix (i.e., different Phi values)
Single letter code for the optimal aa. If NULL, integrates across aa.
Color palette to use from RColorBrewer
Line width
Include stop codons
If "histogram", do a histogram plot; if "density", do a density plot
If TRUE, plot W; if FALSE, plot S (= 1 - W)
The genetic code
Other paramters to pass to plot()