simTwoTheta: Simulate expression data for the two-theta EVE model, i.e....

In Jmendo12/evemodel: Expression Variance Evolution Model

Description

Simulate expression data for the two-theta EVE model, i.e. two different thetas assigned to specific edges

Usage

simTwoTheta(n, tree, colSpecies, isTheta2edge, theta1, theta2, sigma2,
  alpha, beta)

Arguments

n	Number of "genes" to simulate
tree	Phylogeny
colSpecies	A character vector with same length as columns in returned expression matrix, specifying the species, i.e. tip labels in the phylogeny, for the corresponding column.
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)
theta1	Value of the theta1 parameter
theta2	Value of the theta2 parameter
sigma2	Value of the sigma2 parameter
alpha	Value of the alpha parameter
beta	Value of the beta parameter

Value

Matrix of simulated gene expression values with samples in columns and genes in rows

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