A generic function to calculate the number of fixed and shared polymorphisms.
# S4 method for GENOME
calc.fixed.shared(object,
subsites=FALSE,
new.populations=FALSE,
fixed.threshold=1,
fixed.threshold.fst=1)
An object of class "GENOME"
"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.
list of populations. default=FALSE
Polymorphisms are considered as fixed >= threshold value
Polymorphisms are considered as fixed >= threshold 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 |
The nucleotide diversities have to be devided by GENOME.class@n.sites
to give diversities per site.
[x]
# 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