Description Usage Arguments Value
Test for shift in theta at branch
1 2 | twoThetaTest(tree, gene.data, isTheta2edge,
colSpecies = colnames(gene.data), ...)
|
tree |
Species phylogeny (phylo object) |
gene.data |
A matrix of expression values with samples in columns and genes in rows |
isTheta2edge |
Logical vector with same length as number of edges in tree specifying whether the corresponding edge theta parameter should be theta2 (TRUE) or theta1 (FALSE) |
colSpecies |
A character vector with same length as columns in gene.data, specifying the species for the corresponding column. |
... |
parameters passed to |
List with: indivBetaRes (results from fitOneTheta
)
oneThetaRes: results from fitOneTheta
twoThetaRes: results from fitTwoTheta
LRT: log likelihood ratio test statistic between the two-theta and one-theta model
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.