Objects from the Class. A virtual Class: No objects may be created from it. Class 'GARCHfit, directly.Class 'rGARCH, by class 'GARCHfit, distance 2. Object of class 'vector' Holds data on the fitted model. Object of class 'vector' The model specification common to all. Working with R studio is highly recommended and will be more clearly outlined in this post. R does not support working with spatial data straight out of the box so there are a couple of packages that need to be downloaded to get R working with spatial data. The two packages required are ‘sp' and ‘rgdal'. Depends R (= 3.2.5), lme4 (= 1.1-10), stats, methods Imports numDeriv, MASS, ggplot2. SlotNames(bm) carrots Consumer Preference Mapping of Carrots Description In a consumer study 103 consumers scored their preference of 12 danish carrot types on a scale from 1 to 7. Moreover the consumers scored the degree of sweetness, bitterness. Objects from the Class. A virtual Class: No objects may be created from it. Class 'GARCHfit, directly.Class 'rGARCH, by class 'GARCHfit, distance 2. Object of class 'vector' Holds data on the fitted model.
Seurat is a popular R package that is designed for QC, analysis, and exploration of single cell RNA-seq data. Seurat aims to enable users to identify and interpret sources of heterogeneity from single cell transcriptomic measurements, and to integrate diverse types of single cell data. Further, the authors provide several tutorials, on their website.
Download and expand the expression_tables_cellrangerV3.zip file to extract the single cell matrix files for the three samples. These are isolated mouse cells ran on the 10X genomics platform (3' expression V2) for single cell RNA sequencing, sequenced with UC Davis on 1 HiSeq 4000 lane. The Experiment contains 3 samples, each merged from 2 original samples and then randomly subsamples to 1000 cells each.
The three samples are, Dorsal root ganglion neurons :At weaning, Ttpa+/+ mice were fed a normal diet (35 mg of dl-α-tocopheryl acetate/kg, vitE+) while and Ttpa-/- mice were fed either an α-TOH-deficient diet (<10 mg of dl-α-tocopheryl acetate/kg, vitE-), or α-TOH-supplemented diet (600 mg of dl-α-tocopheryl acetate/kg, vitE+++) diet.
- UCD_Adj_VitE - normal + Vitamin E
- UCD_Supp_VitE - Vitamin E supplimented by diet.
- UCD_VitE_Def - Vitamin E deficient animals
We start with loading needed libraries for R, at this time all we need is the package Seurat.
Load the Cell Ranger Matrix Data and create the base Seurat object.
R Slotnames
Cell Ranger provides a function cellranger aggr that will combine multiple samples into a single matrix file. However, when processing data in R and Seurat this is unnecessary and we can aggregate them in R.
Seurat provides a function Read10X to read in 10X data folder. First we read in data from each individual sample folder. First, we initialize the Seurat object (CreateSeuratObject) with the raw (non-normalized data). Keep all genes expressed in >= 10 cells. Keep all cells with at least 200 detected genes. Also extracting sample names, calculating and adding in the metadata mitochondrial percentage of each cell. Adding in the metadata batchid and cell cycle. Finally, saving the raw Seurat object.
The percentage of reads that map to the mitochondrial genome
- Low-quality / dying cells often exhibit extensive mitochondrial contamination.
- We calculate mitochondrial QC metrics with the PercentageFeatureSet function, which calculates the percentage of counts originating from a set of features.
- We use the set of all genes, in mouse these genes can be identified as those that begin with ‘mt', in human data they begin with MT.
Calculate cell cycle, add to meta data
Using the package scran, get the mouse cell cycle markers and a mapping of m
Lets create a fake batch metadata (used in part 3), Here we determine UCD_Adj_VitE is from one batch and UCD_Adj_VitE/UCD_Adj_VitE are from a second battch
Here we build a new metadata variable ‘batchid' which can be used to specify treatment groups.
Lets spend a little time getting to know the Seurat object.
The Seurat object is the center of each single cell analysis. It stores all information associated with the dataset, including data, annotations, analyses, etc. The R function slotNames can be used to view the slot names within an object.
We can then view the data within a slot with the @ operator.
Question(s)
Setnames In R
- What slots are empty, what slots have data?
- What columns are available in meta.data?
- Look up the help documentation for subset?
Finally, save the original object, write out a tab-delimited table that could be read into excel, and view the object.
Get the next Rmd file
Session Information
