Last data update: 2014.03.03

R: Accessors for the 'counts' slot of a DESeqDataSet object.
countsR Documentation

Accessors for the 'counts' slot of a DESeqDataSet object.

Description

The counts slot holds the count data as a matrix of non-negative integer count values, one row for each observational unit (gene or the like), and one column for each sample.

Usage

## S4 method for signature 'DESeqDataSet'
counts(object, normalized = FALSE,
  replaced = FALSE)

## S4 replacement method for signature 'DESeqDataSet,matrix'
counts(object) <- value

Arguments

object

a DESeqDataSet object.

normalized

logical indicating whether or not to divide the counts by the size factors or normalization factors before returning (normalization factors always preempt size factors)

replaced

after a DESeq call, this argument will return the counts with outliers replaced instead of the original counts, and optionally normalized. The replaced counts are stored by DESeq in assays(object)[['replaceCounts']].

value

an integer matrix

Author(s)

Simon Anders

See Also

sizeFactors, normalizationFactors

Examples


dds <- makeExampleDESeqDataSet(m=4)
head(counts(dds))

dds <- estimateSizeFactors(dds) # run this or DESeq() first
head(counts(dds, normalized=TRUE))

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(DESeq2)
Loading required package: S4Vectors
Loading required package: stats4
Loading required package: BiocGenerics
Loading required package: parallel

Attaching package: 'BiocGenerics'

The following objects are masked from 'package:parallel':

    clusterApply, clusterApplyLB, clusterCall, clusterEvalQ,
    clusterExport, clusterMap, parApply, parCapply, parLapply,
    parLapplyLB, parRapply, parSapply, parSapplyLB

The following objects are masked from 'package:stats':

    IQR, mad, xtabs

The following objects are masked from 'package:base':

    Filter, Find, Map, Position, Reduce, anyDuplicated, append,
    as.data.frame, cbind, colnames, do.call, duplicated, eval, evalq,
    get, grep, grepl, intersect, is.unsorted, lapply, lengths, mapply,
    match, mget, order, paste, pmax, pmax.int, pmin, pmin.int, rank,
    rbind, rownames, sapply, setdiff, sort, table, tapply, union,
    unique, unsplit


Attaching package: 'S4Vectors'

The following objects are masked from 'package:base':

    colMeans, colSums, expand.grid, rowMeans, rowSums

Loading required package: IRanges
Loading required package: GenomicRanges
Loading required package: GenomeInfoDb
Loading required package: SummarizedExperiment
Loading required package: Biobase
Welcome to Bioconductor

    Vignettes contain introductory material; view with
    'browseVignettes()'. To cite Bioconductor, see
    'citation("Biobase")', and for packages 'citation("pkgname")'.

> png(filename="/home/ddbj/snapshot/RGM3/R_BC/result/DESeq2/counts.Rd_%03d_medium.png", width=480, height=480)
> ### Name: counts
> ### Title: Accessors for the 'counts' slot of a DESeqDataSet object.
> ### Aliases: counts counts,DESeqDataSet-method
> ###   counts<-,DESeqDataSet,matrix-method
> 
> ### ** Examples
> 
> 
> dds <- makeExampleDESeqDataSet(m=4)
> head(counts(dds))
      sample1 sample2 sample3 sample4
gene1      61     110      34      88
gene2       8      18      13       9
gene3      23      32      27       8
gene4       6       1       1       2
gene5       1       0       9       2
gene6      30       0      12       4
> 
> dds <- estimateSizeFactors(dds) # run this or DESeq() first
> head(counts(dds, normalized=TRUE))
         sample1     sample2    sample3   sample4
gene1 59.9058084 106.2842100 33.9103838 87.421140
gene2  7.8564995  17.3919616 12.9657350  8.940798
gene3 22.5874359  30.9190429 26.9288342  7.947376
gene4  5.8923746   0.9662201  0.9973642  1.986844
gene5  0.9820624   0.0000000  8.9762781  1.986844
gene6 29.4618730   0.0000000 11.9683707  3.973688
> 
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>