phylosim (version 3.0.5)

omegaVarM0.CodonSequence: The M0 (one-ratio) model of variable omega ratios among sites

Description

The M0 (one-ratio) model of variable omega ratios among sites.

This method sets the omega site-process specific parameter in the specified range to values sampled from the M0 (one-ratio) model of variable omega ratios among sites.

Distribution of omega values:

	CATEGORY	PROBABILITY
omega		1
       

Usage

# S3 method for CodonSequence
omegaVarM0(this, process, omega, index, ...)

Arguments

this

A CodonSequence object.

process

A process object inheriting from GY94.

omega

The fixed omega value.

index

A vector of positions.

...

Not used.

Value

Invisible TRUE.

References

Yang, Z., Nielsen, R., Goldman, N., Pedersen Krabbe, A-M. (2000) Codon-Substitution Models for Heterogeneous Selection Pressure at Amino Acid Sites - Genetics 155:431-449 http://bit.ly/bvjucn

Goldman, N., Yang, Z. (1994) A codon-based model of nucleotide substitution for protein-coding DNA sequences - Mol Biol Evol 11(5):725-36 http://bit.ly/aSVEoa

See Also

For more information see CodonSequence.

Examples

# NOT RUN {
	# create a GY94 object
	p<-GY94(kappa=2)
	# create a CodonSequence object, attach process p
	s<-CodonSequence(length=20, processes=list(list(p)))
	# sample states
	sampleStates(s)
	# sample omegas in range 1:5 from model M0
	omegaVarM0(s,p,omega=2,1:5)
	# get omega values
	getOmegas(s,p)
	# get a histogram of omega values in range 1:5
	omegaHist(s,p,breaks=50,1:5)
 
# }