Question: Basemean always multiples of 0.5 in deseq2 results
0
gravatar for batmaiden
16 months ago by
batmaiden0
batmaiden0 wrote:

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 • 231 views
ADD COMMENTlink modified 16 months ago by Michael Love25k • written 16 months ago by batmaiden0
Answer: Basemean always multiples of 0.5 in deseq2 results
0
gravatar for Michael Love
16 months ago by
Michael Love25k
United States
Michael Love25k wrote:

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 COMMENTlink modified 16 months ago • written 16 months ago by Michael Love25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 218 users visited in the last hour