omegaHist.CodonSequence: Plot a histogram of omega values from a range

In phylosim: Flexible Simulations of Biological Sequence Evolution

Description

Plot a histogram of omega values from a range.

Usage

## S3 method for class 'CodonSequence'
omegaHist(this, process, breaks, index, ...)

Arguments

this	A CodonSequence object.
process	A process object inheriting from GY94.
breaks	breaks parameter for hist().
index	A vector of positions.
...	Not used.

Value

The CodonSequence object (invisible).

Author(s)

Botond Sipos, Gregory Jordan

See Also

For more information see CodonSequence.

Examples

	# create a GY94 process
	p<-GY94()
	# create a CodonSequence object,
	# attach a process p
	s<-CodonSequence(length=20,processes=list(list(p)))
       # set omega values through omegaVarM2.CodonSequence
       omegaVarM2(s,p,p0=0.5,p1=0.2,omega=1.5)
       # get a histogram of omega values from the range 1:15
       omegaHist(s,p,breaks=10,1:15)
 

phylosim documentation built on Nov. 22, 2019, 1:07 a.m.