Error in DUTest about sierra

[rj-patrick]

Hi @soobeomlee,

Can you confirm whether you are able to run the example in the vignette without error?

Would you be able to show the command that you are running and any other text printed by DUTest?

Cheers,
Ralph

[soobeomlee]

Hi, thanks for the quick reply.
I ran the example in the vignette and got the same error. I believe the version difference of some package might be a problem. I will try to reinstall and run. please let me know when you get a hint for this issue.

FYI, here's traceback.
traceback()
11: stop(gettextf("assignment of an object of class %s is not valid for slot %s in an object of class %s; is(value, "%s") is not TRUE",
dQuote(valueClass), sQuote(name), dQuote(cl), slotClass),
domain = NA)
10: checkSlotAssignment(object, name, value)
9: slot<-(*tmp*, what, value = slot(value, what))
8: asMethod(object, Class, value)
7: as<-(*tmp*, Classi, value = obj)
6: initialize(value, ...)
5: initialize(value, ...)
4: new("DEXSeqResults", LRTresults, modelFrameBM = object@modelFrameBM,
sampleData = sampleAnnotation(object), dispersionFunction = object@dispersionFunction)
3: DEXSeq::DEXSeqResults(dxd)
2: apply_DEXSeq_test_seurat(apa.seurat.object = peaks.object, population.1 = population.1,
population.2 = population.2, exp.thresh = exp.thresh, fc.thresh = fc.thresh,
adj.pval.thresh = adj.pval.thresh, num.splits = num.splits,
seed.use = seed.use, feature.type = feature.type, include.annotations = include.annotations,
filter.pA.stretch = filter.pA.stretch, verbose = verbose,
do.MAPlot = do.MAPlot, return.dexseq.res = return.dexseq.res,
ncores = ncores)
1: DUTest(peaks.seurat, population.1 = "F-SL", population.2 = "EC1",
exp.thresh = 0.1, feature.type = c("UTR3", "exon"))

[rj-patrick]

Hi @soobeomlee,

Sorry for the slow response. I have tried running Sierra using the same software versions for the main requirements (DEXSeq version 1.36), Seurat v4 and don't have any issues. Only possibility I can think of is whether there is a conflict between your R and Bioconductor versions. I don't know how likely this is to be a problem, but it appears you are using Bioconductor version 3.12, which is supposed to be for R version 4.03, so I wonder if you were to downgrade to version 3.11 (or upgrade your R version) whether this would solve the issue?

[soobeomlee]

Awesome! It works on Bioconductor 3.11!

Thanks a lot!

[rj-patrick]

No worries, let us know if you have any more issues.