calc.fixed.shared-methods: Fixed and shared polymorphisms
Description
A generic function to calculate the number of fixed and shared polymorphisms.
Usage
# S4 method for GENOME
calc.fixed.shared(object,
subsites=FALSE,
new.populations=FALSE,
fixed.threshold=1,
fixed.threshold.fst=1)
Arguments
object
An object of class "GENOME"
subsites
"transitions": SNPs that are transitions.
"transversions": SNPs that are transversions.
"syn": synonymous sites.
"nonsyn": nonsynonymous sites.
"exon": SNPs in exon regions.
"intron": SNPs in intron regions.
"coding": SNPs in coding regions (CDS).
"utr": SNPs in UTR regions.
"gene": SNPs in genes.
"intergenic" : SNPs in intergenic regions.
new.populations
list of populations. default=FALSE
fixed.threshold
Polymorphisms are considered as fixed >= threshold value
fixed.threshold.fst
Polymorphisms are considered as fixed >= threshold value
Value
Returned value is a modified object of class "GENOME"
---------------------------------------------------------
The following slots will be modified in the "GENOME" object
---------------------------------------------------------
|
|
Slot |
Reference |
Description |
|
1. |
n.fixed.sites |
[x] |
Number of fixed sites |
|
2. |
n.shared.sites |
[x] |
Number of shared sites |
|
3. |
n.monomorphic.sites |
[x] |
Number of monomorphic sites |
Details
The nucleotide diversities have to be devided by GENOME.class@n.sites to give diversities per site.
Examples
Run this code# NOT RUN {
# GENOME.class <- readData("\home\Alignments")
# set.populations
# GENOME.class <- calc.fixed.shared(GENOME.class)
# }
Run the code above in your browser using DataLab