0

我正在运行 DEXSeq 分析，但是当我到达 DEXSeqDataSetFromSE 时，我收到以下错误。我已经多次使用这个脚本，它运行良好，但这次我添加了另外两个生物复制，我得到了这个错误。我将不胜感激。

DEXSeq：DESeqDataSet 中的错误（rse，设计，ignoreRank = TRUE

library(GenomicRanges)
library(GenomicFeatures)
library(GenomicAlignments)
library(DEXSeq)
 #Download the human genome we aligned to hg19        
library(TxDb.Hsapiens.UCSC.hg19.knownGene)
hse <- TxDb.Hsapiens.UCSC.hg19.knownGene

#Collapse Gene models into counting bins. The exonicParts object contains a GRanges object with our counting bins. We use it to count the number of read fragments that overlap with the bins by means of the function
exonicParts = exonicParts( hse, linked.to.single.gene.only = TRUE )

# show the BAM files 
fls <- list.files (pattern='bam$', full=TRUE )
bamlst <- BamFileList(fls, index=character(), yieldSize=50000, obeyQname=TRUE )


#count reads overlapping bins 

SE <- summarizeOverlaps (exonicParts, bamlst, mode='Union', singleEnd=FALSE, ignore.strand=TRUE, inter.feature=FALSE, fragments=TRUE )
colData (SE)$condition <- c("VP-35", "WT", "VP-35", "WT")
#Build DEXSeq Dataset
DEXSeqDataSetFromSE( SE, design= ~ sample + exon + condition:exon )
DEXSeq : Error in DESeqDataSet(rse, design, ignoreRank = TRUE
dxd <- DEXSeqDataSetFromSE( SE, design= ~ sample + exon + condition:exon )
dxd <- estimateSizeFactors (dxd )
dxd <- estimateDispersions (dxd )
4

0 回答 0