createGeneModelFromTfList: Create a model for a target gene using a Trena object

In PriceLab/trena: Fit transcriptional regulatory networks using gene expression, priors, machine learning

Create a model for a target gene using a Trena object

Description

Create a model for a target gene using a Trena object

Usage

## S4 method for signature 'Trena'
createGeneModelFromTfList(obj, targetGene, solverNames, tfList, mtx)

Arguments

obj	An object of class Trena
targetGene	The name of a target gene to use for building a model
solverNames	A character vector containing the solver names to be used for building the model
tfList	A character list, often the gene symbols for known transcription factors
mtx	An assay matrix of expression data

Value

A data frame containing the gene model

Examples

if(interactive()){  # takes too long for the bioconductor build
   # Create a Trena object for human and make a gene model for "MEF2C" using a footprint filter
   trena <- Trena("hg38")
   chromosome <- "chr5"
   mef2c.tss <- 88904257
   loc.start <- mef2c.tss - 1000
   loc.end   <- mef2c.tss + 1000

   database.filename <- system.file(package="trena", "extdata", "mef2c.neigborhood.hg38.footprints.db")
   database.uri <- sprintf("sqlite://%s", database.filename)
   sources <- c(database.uri)
   load(system.file(package="trena", "extdata/ampAD.154genes.mef2cTFs.278samples.RData"))

   model.mef2c <- createGeneModelFromTfList(trena, "MEF2C", c("lasso","ridge","randomforest"),
                                            tfList=c())
   } # if interactive

PriceLab/trena documentation built on March 16, 2023, 10:01 a.m.