Hi,
I have recently started using SeqGSEA to analyse RNAseq data from 12 tumors. (Paired-end, 70-90 mill reads, mapped with tophat, counted with htseq-count).
I would like to compare two of these tumors (samples v6 and v11, see below) to the other ten regarding certain gene sets. However, I run into the following issue when following 6.3 in the manual (DE only):
Genes with read count 0 across all 12 samples have been removed from my counts table.
> nrow(counts)
[1] 45418
> head(counts)
v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12
ENSG00000000003 43888 15516 15307 4287 17635 3342 3904 7235 22298 12464 4476 10353
ENSG00000000005 2238 35 27 23 81 0 12 46 335 39 4 56
ENSG00000000419 2781 3032 1805 2644 2027 2651 3383 2230 2416 2070 2503 1381
ENSG00000000457 1212 727 699 962 707 626 1054 944 950 647 507 493
ENSG00000000460 524 538 396 640 415 1467 852 1073 568 419 409 399
ENSG00000000938 173 1349 616 467 280 1769 987 1159 415 448 1591 679
> label <- as.factor(c(0,0,0,0,0,1,0,0,0,0,1,0))
> DEG <- newCountDataSet(counts,label)
> DEG <- estimateSizeFactors(DEG)
> DEG <- estimateDispersions(DEG, method="pooled", fitType="local")
> DEGres <- DENBStat4GSEA(DEG)
> permuteMat <- genpermuteMat(label, times=perm.times)
> DEpermNBstat <- DENBStatPermut4GSEA(DEG, permuteMat)
Error in { :
task 1 failed - "Parametric dispersion fit failed. Try a local fit and/or a pooled estimation. (See '?estimateDispersions')"
I am using estimateDispersions and estimateSizeFactors instead of runDESeq in order to use fitType=local and method=pooled, but am still experiencing trouble. Any idea what I'm doing wrong? Anyone else having this problem? Any help would be greatly appreciated!
Cheers,
Thale
I have recently started using SeqGSEA to analyse RNAseq data from 12 tumors. (Paired-end, 70-90 mill reads, mapped with tophat, counted with htseq-count).
I would like to compare two of these tumors (samples v6 and v11, see below) to the other ten regarding certain gene sets. However, I run into the following issue when following 6.3 in the manual (DE only):
Genes with read count 0 across all 12 samples have been removed from my counts table.
> nrow(counts)
[1] 45418
> head(counts)
v1 v2 v3 v4 v5 v6 v7 v8 v9 v10 v11 v12
ENSG00000000003 43888 15516 15307 4287 17635 3342 3904 7235 22298 12464 4476 10353
ENSG00000000005 2238 35 27 23 81 0 12 46 335 39 4 56
ENSG00000000419 2781 3032 1805 2644 2027 2651 3383 2230 2416 2070 2503 1381
ENSG00000000457 1212 727 699 962 707 626 1054 944 950 647 507 493
ENSG00000000460 524 538 396 640 415 1467 852 1073 568 419 409 399
ENSG00000000938 173 1349 616 467 280 1769 987 1159 415 448 1591 679
> label <- as.factor(c(0,0,0,0,0,1,0,0,0,0,1,0))
> DEG <- newCountDataSet(counts,label)
> DEG <- estimateSizeFactors(DEG)
> DEG <- estimateDispersions(DEG, method="pooled", fitType="local")
> DEGres <- DENBStat4GSEA(DEG)
> permuteMat <- genpermuteMat(label, times=perm.times)
> DEpermNBstat <- DENBStatPermut4GSEA(DEG, permuteMat)
Error in { :
task 1 failed - "Parametric dispersion fit failed. Try a local fit and/or a pooled estimation. (See '?estimateDispersions')"
I am using estimateDispersions and estimateSizeFactors instead of runDESeq in order to use fitType=local and method=pooled, but am still experiencing trouble. Any idea what I'm doing wrong? Anyone else having this problem? Any help would be greatly appreciated!
Cheers,
Thale
Comment