filter_genes_without_1to1_homolog: filter_genes_without_1to1_homolog

Description Usage Arguments Value Examples

View source: R/filter_genes_without_1to1_homology.r

Description

filter_genes_without_1to1_homolog Takes the filenames of celltype_data files, loads them, and drops any genes which don't have a 1:1 homolog based on biomart. The new files are saved to disc, appending '_1to1only' to the file tag of the previous file.

Usage

1

Arguments

filenames

Array of filenames for sct_data saved as .Rda files

Value

Array of filenames included the ones with only 1:1 homologs

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
library(ewceData)
# Load the single cell data
cortex_mrna <- cortex_mrna()
expData <- cortex_mrna$exp
expData <- expData[1:100, ] # Use only a subset to keep the example quick
l1 <- cortex_mrna$annot$level1class
l2 <- cortex_mrna$annot$level2class
annotLevels <- list(level1class = l1, level2class = l2)
fNames_ALLCELLS <- generate_celltype_data(exp = expData,
    annotLevels = annotLevels,
    groupName = "allKImouse")
fNames_ALLCELLS <- filter_genes_without_1to1_homolog(fNames_ALLCELLS)

NathanSkene/EWCE documentation built on June 19, 2021, 5:40 a.m.