This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Error with KissDE exploreResults()

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

0 answers

No answers yet.

Log in to answer this question.