Description Usage Arguments Details Value Author(s) References See Also Examples
View source: R/genomedistances.R
This function computes weights for gene cluster according to their distribution in a pan-genome.
1 | geneWeights(pan.matrix, type = c("shell", "cloud"))
|
pan.matrix |
A pan-matrix, see |
type |
A text indicating the weighting strategy. |
When computing distances between genomes or a PCA, it is possible to give weights to the different gene clusters, emphasizing certain aspects.
As proposed by Snipen & Ussery (2010), we have implemented two types of weighting: The default "shell" type means gene families occuring frequently in the genomes, denoted shell-genes, are given large weight (close to 1) while those occurring rarely are given small weight (close to 0). The opposite is the "cloud" type of weighting. Genes observed in a minority of the genomes are referred to as cloud-genes. Presumeably, the "shell" weighting will give distances/PCA reflecting a more long-term evolution, since emphasis is put on genes who have just barely diverged away from the core. The "cloud" weighting emphasizes those gene clusters seen rarely. Genomes with similar patterns among these genes may have common recent history. A "cloud" weighting typically gives a more erratic or ‘noisy’ picture than the "shell" weighting.
A vector of weights, one for each column in pan.matrix
.
Lars Snipen and Kristian Hovde Liland.
Snipen, L., Ussery, D.W. (2010). Standard operating procedure for computing pangenome trees. Standards in Genomic Sciences, 2:135-141.
1 | # See examples for distManhattan
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.