In flaviusb/homosapienDEE2CellScore: Example Data Package for CellScore

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

Introduction

This data package contains code to generate various filtered forms of a normalised dataset based on downloading the accessions in inst/hsapiens_colData_transitions_v3.5.csv for species "hsapiens" from DEE2, and it also has a way (for reproducibility and also ease of use) to use AnnotationHub (as the data files were too large to fit inside of the data package itself) to fetch data pregenerated at the initial package development time.

Installation

We assume you already have R installed and configured. If you do not have Bioconductor installed, you can install it like so:

if (!requireNamespace("BiocManager", quietly=TRUE))
    install.packages("BiocManager")
BiocManager::install()

Then, if you have not installed this package and the dependencies from Bioconductor, you can install them like so:

BiocManager::install(c("CellScore", "homosapienDEE2CellScore", "devtools", "getDEE2", "SummarizedExperiment"))

Setup

In order to make the needed libraries accessible, you will need the following setup:

library(DESeq2)
library(S4Vectors)
library(Biobase)
library(SummarizedExperiment)
library(getDEE2)
library(devtools)
library(CellScore)
library(homosapienDEE2CellScore)

Getting and using the data

All of the permutations of filtering and normalising the data can be downloaded in a labeled list like so:

the_data<-downloadAllTheData()

We can then get the data we are after - non-normalised data including samples with quality control warnings - and use it in CellScore to calculate the on/off score for cell transitions from fibroblast to embryonic stem cells:

sm <- the_data$HomosapienDEE2_QC_WARN_Raw
## We could have just run `sm <- homosapienDEE2CellScore::readInSEZip(homosapienDEE2CellScore::HomosapienDEE2_QC_PASS_Raw())`
## instead of downloading all the data.

# Here we want to analyse all of the raw data to calculate the
# on/off score for cell transitions from fibroblast to embryonic stem cells
test1 <- sm[, sm$category == 'test']
standard <- sm[, sm$category == 'standard']
sm1 <- cbind(test1, standard)
cell.change <- data.frame(start=c("FIB"), test=c("nESC"), target=c("ESC"))
group.OnOff <- OnOff(sm1, cell.change, out.put="marker.list")