Hi,
I was going to try DEXSeq for the splicing analysis, but I got stuck at the beginning.
I got 'subscript out of bounds' errors when tried 'DEXSeqDataSetFromHTSeq'. The error seems to tell me that sample or column numbers were not matching with something else.
I don't really understand it. Please, somebody tell me what did I do wrong?
I was going to try DEXSeq for the splicing analysis, but I got stuck at the beginning.
I got 'subscript out of bounds' errors when tried 'DEXSeqDataSetFromHTSeq'. The error seems to tell me that sample or column numbers were not matching with something else.
I don't really understand it. Please, somebody tell me what did I do wrong?
Code:
library("DEXSeq")
sampleTable <- data.frame(
rownames = c( "IVF1","IVF2","IVF3",
"cSCNT1","cSCNT2","cSCNT3","cSCNT4",
"fSCNT1","fSCNT2","fSCNT3","fSCNT4"),
countFiles = c( "~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/IVF1.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/IVF2.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/IVF3.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/cSCNT1.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/cSCNT2.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/cSCNT3.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/cSCNT4.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/fSCNT1.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/fSCNT2.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/fSCNT3.counts",
"~/00-NGS/RNAseq/bov/Old_Pico/Ensembl/DEXseq/fSCNT4.counts"),
condition = c( rep("IVF", 3), rep("cSCNT", 4),rep("fSCNT", 4) ),
type = c("single","single","single","single","single","single","single","single","single","single","single"))
sampleTable
flattenedfile="/home/bio1/00-NGS/RNAseq/bov/Old_Pico/Annotation/DEXseq_Bovine_Ensembl.gtf"
dxd = DEXSeqDataSetFromHTSeq(
countFiles,
sampleData=sampleTable,
design= ~ sample + exon + condition:exon,
flattenedfile=flattenedfile)