Question: (Closed) help with start DESeq
0
gravatar for batmaiden
5 weeks ago by
batmaiden0
batmaiden0 wrote:

Hi,

I want to analyze my RNA-seq data in DESeq2. I have the following code:

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

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

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

library( "GenomicAlignments" ) se <- summarizeOverlaps( exonsByGene, 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= ~FALSE )

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

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

dds <- ddsFullCountTable as.data.frame( colData(dds) )

dds <- DESeq(dds, betaPrior = FALSE)

And function DESeq gives an error: Error in solve.default(qr.R(qrx)) : 'a' (1 x 0) must be square

Help me, please!

rna-seq software error deseq • 102 views
ADD COMMENTlink modified 5 weeks ago • written 5 weeks ago by batmaiden0

the question is solved

ADD REPLYlink written 5 weeks ago by batmaiden0

Hello batmaiden!

Questions similar to yours can already be found at:

We have closed your question to allow us to keep similar content in the same thread.

If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.

Cheers!

PS: For line command Deseq2 usage help, please see https://support.bioconductor.org/
ADD REPLYlink written 4 weeks ago by Jennifer Hillman Jackson25k
Please log in to add an answer.
The thread is closed. No new answers may be added.

Help
Access

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