twoThetaTest: Test for shift in theta at branch

In Jmendo12/evemodel: Expression Variance Evolution Model

Description

Test for shift in theta at branch

Usage

twoThetaTest(tree, gene.data, isTheta2edge,
  colSpecies = colnames(gene.data), ...)

Arguments

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 fitOneTheta and fitTwoTheta

Value

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

Jmendo12/evemodel documentation built on June 17, 2020, 8:56 p.m.