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 !