Basemean always multiples of 0.5 in deseq2 results
1
0
Entering edit mode
batmaiden • 0
@batmaiden-16484
Last seen 3.2 years ago

Hi,

I analyzing RNA-seq with two samples: control and treatment. I use the following sequence of software: sortmerna -> hisat2 -> DESeq2. This is my R code for DESeq2:

library( "GenomicFeatures" )
hse <- makeTxDbFromGFF( "~/Homo_sapiens.GRCh38.92.gtf", format="gtf" )
trans <- transcriptsBy( hse, by="gene" )

fls <- list.files( "./", pattern="bam$", full=TRUE )

library( "Rsamtools" )
bamLst <- BamFileList( fls, yieldSize=100000 )

library( "GenomicAlignments" )
se <- summarizeOverlaps( trans, bamLst,
mode="Union",
singleEnd=FALSE,
ignore.strand=TRUE,
fragments=TRUE )

library( "DESeq2" )

sampleInfo <- read.csv( "sample.csv" )
sampleInfo <- DataFrame( sampleInfo )
seIdx <- match(colnames(se), sampleInfo$run)

colData(se) <- cbind( colData(se), sampleInfo)

ddsFull <- DESeqDataSet( se, design= ~run )

countdata <- assay(se)
coldata <- colData (se)


ddsFullCountTable <- DESeqDataSetFromMatrix(
countData = countdata,
colData = coldata,
design = ~run)

dds <- ddsFullCountTable

dds <- DESeq(dds)
res <- results( dds )

write.csv( as.data.frame(res), file="results.csv" )
Everything works, but there is a problem with the results: all values ​​in the column "Basemean" are multiples of 0.5 (apparently some setting). What parameters need to be changed? Thank you so much!
deseq results baseMean • 376 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 2 hours ago
United States

I wouldn’t recommend counting in transcripts unless you want intronic and incompatibly spliced reads to count for a gene. I’d recommend using the tximport pipeline, htseq-count or featureCounts with exon-by-gene features.

Otherwise, there isn’t a reason that baseMean would be rounded to 0.5 — try using makeExampleDESeqDataSet()

ADD COMMENT

Login before adding your answer.

Traffic: 283 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6