solve.Pearson: Run the Pearson Solver

In TReNA: Fit transcriptional regulatory networks using gene expression, priors, machine learning

Description

Given a TReNA object with Pearson as the solver, use the cor function to estimate coefficients for each transcription factor as a perdictor of the target gene's expression level. This method should be called using the solve method on an appropriate TReNA object.

Usage

## S4 method for signature 'PearsonSolver'
run(obj, target.gene, tfs, tf.weights = rep(1,
  length(tfs)), extraArgs = list())

Arguments

obj	An object of class Solver with "pearson" as the solver string
target.gene	A designated target gene that should be part of the mtx.assay data
tfs	The designated set of transcription factors that could be associated with the target gene.
tf.weights	A set of weights on the transcription factors (default = rep(1, length(tfs)))
extraArgs	Modifiers to the Pearson solver

Value

The set of Pearson Correlation Coefficients between each transcription factor and the target gene.

See Also

cor, PearsonSolver

Other solver methods: run,BayesSpikeSolver-method, run,EnsembleSolver-method, run,LassoPVSolver-method, run,LassoSolver-method, run,RandomForestSolver-method, run,RidgeSolver-method, run,SpearmanSolver-method, run,SqrtLassoSolver-method, solve,TReNA-method

Examples

# Load included Alzheimer's data, create a TReNA object with Bayes Spike as solver, and solve
load(system.file(package="TReNA", "extdata/ampAD.154genes.mef2cTFs.278samples.RData"))
trena <- TReNA(mtx.assay = mtx.sub, solver = "pearson")
target.gene <- "MEF2C"
tfs <- setdiff(rownames(mtx.sub), target.gene)
tbl <- solve(trena, target.gene, tfs)

TReNA documentation built on Nov. 17, 2017, 12:35 p.m.