Error with KissDE exploreResults()
0
0
Entering edit mode
12 weeks ago

Hello,

For some time now, the function exploreResults() of kissDE doesn't seem to work, throwing an error :

Error in names(x) <- value : 
  'names' attribute [16] must be the same length as the vector [10]

I know there is a warning on the github page : WARNING: exploreResults will not work correctly on rds file from previous release.

So I ran the whole analysis starting from kissplice (v2.6.7) on a small dataset of 4 samples. Default parameter for kissplice, except for --type1-only. Default parameter for kissplice2refgenome (v2.0.7) as well, except for --pairedEnd and --readLength.

For kissDE, I did the analysis either on the 1.23.0 version from github, or the 1.24.0 version from Bioconductor. In both cases, there is the error. I simply used :

myCounts <- kissplice2counts(k2rg_file, counts = 2, pairedEnd = TRUE, k2rg = TRUE, 
                             remove = c("deletion", "insertion", "-", "indel"))
myConditions <- c("MUT", "MUT", "WT", "WT")
diffSplicing <- diffExpressedVariants(myCounts, myConditions, technicalReplicates = FALSE,
                                      filterLowCountsVariants = 10, nbCore = 32)
writeOutputKissDE(diffSplicing, output = "kissDE_test.tab", adjPvalMax = 0.05, dPSImin = 0.1)
exploreResults(rdsFile = "kissDE_test.tab.rds")

I don't know if this come from my end or if it is a known problem, but just in case, I wanted to point it out.

Thanks !

kissplice kissDE • 186 views
ADD COMMENT

Login before adding your answer.

Traffic: 1445 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