## ----setup, echo = FALSE--------------------------------------------------- library(BiocStyle) knitr::opts_chunk$set(tidy = FALSE, cache = TRUE, autodep = TRUE) ## ----samples--------------------------------------------------------------- dataFilesDir = system.file("extdata", package = "pasilla", mustWork=TRUE) pasillaSampleAnno = read.csv(file.path(dataFilesDir, "pasilla_sample_annotation.csv")) pasillaSampleAnno ## ----dirextData------------------------------------------------------------ dir(system.file("extdata", package = "pasilla", mustWork=TRUE), pattern = ".txt$") ## ----ecs------------------------------------------------------------------- gffFile = file.path(dataFilesDir, "Dmel.BDGP5.25.62.DEXSeq.chr.gff") ## ----read, message = FALSE, warning = FALSE-------------------------------- library("DEXSeq") dxd = DEXSeqDataSetFromHTSeq( countfiles = file.path(dataFilesDir, paste(pasillaSampleAnno$file, "txt", sep=".")), sampleData = pasillaSampleAnno, design= ~ sample + exon + condition:exon, flattenedfile = gffFile) dxd ## ----dxd------------------------------------------------------------------- genesforsubset = readLines(file.path(dataFilesDir, "geneIDsinsubset.txt")) dxd = dxd[geneIDs( dxd ) %in% genesforsubset,] ## ----save, eval = FALSE---------------------------------------------------- # save(dxd, file = file.path("..", "data", "pasillaDEXSeqDataSet.RData")) ## ----sessionInfo----------------------------------------------------------- sessionInfo()