phylosim (version 3.0.5)

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

Description

Plot a histogram of omega values from a range.

Usage

# S3 method for 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).

See Also

For more information see CodonSequence.

Examples

# NOT RUN {
	# 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)
 
# }