Code taken from the DESeq vignette for my own convenience.


exampleFile = system.file ("extra/",package="DESeq")

countsTable = read.delim(exampleFile, header=TRUE, stringsAsFactors=TRUE)

rownames(countsTable) = countsTable$gene

countsTable = countsTable[ , -1]

conds = c("T","T","T","Tb","N","N")

cds = newCountDataSet (countsTable, conds)

cds = cds[,-1]

cds = estimateSizeFactors(cds)


cds <- estimateVarianceFunctions( cds )

res <- nbinomTest( cds, "N", "T" )

resSig <- res[ res$padj < .1, ]

Print Friendly, PDF & Email

Creative Commons License
This work is licensed under a Creative Commons
Attribution 4.0 International License
Posted in RTagged ,

Leave a Reply

Your email address will not be published. Required fields are marked *

This site uses Akismet to reduce spam. Learn how your comment data is processed.