Description Usage Arguments Value Examples
Estimation of parameters of selection strength Gamma Distribution
1 | estParaGamma(exptable = NULL, corrmatinv = NULL)
|
exptable |
an expression level table: column corresponds to median expression value of all biological samples within one taxa_subtaxa group; row corresponds to othologous genes |
corrmatinv |
an inversed correlation matrix between specie |
returns a vector of parameters estimated from gamma distribution
1 2 3 4 5 6 | data('tetraExp')
species.group <- c("Human", "Chimpanzee", "Bonobo", "Gorilla",
"Macaque", "Mouse", "Opossum", "Platypus")
inv.corr.mat <- corrMatInv(tetraExp, taxa = species.group, subtaxa = "Brain")
brain.exptable <- exptabTE(tetraExp, taxa = species.group, subtaxa = "Brain" ,logrithm = TRUE)
gamma.paras <- estParaGamma(exptable =brain.exptable, corrmatinv =inv.corr.mat)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.