estParaWBayesian: Bayesian estimation of expression conservation from previous...

In jingwyang/TreeExp: the R package for analyzing expression evolution based on RNA-seq data

Description

Bayesian estimation of expression conservation from previous data set

Usage

estParaWBayesian(Q_gene = NULL, gammaparas = NULL)

Arguments

Q_gene	a vector specifying parameter Q estimated from estParaQ
gammaparas	a vector specifying parameters estimated from estParaGamma

Value

returns a vector of exptations and confidence intervals of parameter Ws estimated from previous Q estimation

Examples

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)
brain.Q <- estParaQ(brain.exptable, corrmatinv = inv.corr.mat)
brain.post<- estParaWBayesian(brain.Q, gamma.paras)

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