estParaGamma: Estimation of parameters of selection strength Gamma...

Description Usage Arguments Value Examples

Description

Estimation of parameters of selection strength Gamma Distribution

Usage

1
estParaGamma(exptable = NULL, corrmatinv = NULL)

Arguments

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

Value

returns a vector of parameters estimated from gamma distribution

Examples

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)

jingwyang/TreeExp documentation built on June 11, 2019, 6:17 p.m.