This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Monocle3/cicero - data loading in R

I am trying to run cicero on scATAC data, and while loading the data I am getting this error:

Error: CellDataSet 'phenoData' is class 'data.frame' but should be or extend 'AnnotatedDataFrame'

where the command was:

input_cds <- suppressWarnings(newCellDataSet(indata, + featureData = peakinfo, + phenoData = cellinfo, lowerDetectionLimit = 0.1, expressionFamily = VGAM::negbinomial.size()))

cellinfo- barcodes.tsv, peakinfo - peak.bed, indata - matrix.mtx

Please suggest what could be done.

atac sc atac r sequencing next-gen
• 1,303 views • link
written by sidrah.maryam • sidrah.maryam

Hi,

Did you run?

pd <- methods::new("AnnotatedDataFrame", data = cellinfo)

And then give this pd variable to the phenoData option. The new function will coerce the cellinfo data.frame to AnnotatedDataFrame class.

António

• 0 views • link
written by antonioggsousa antonioggsousa

0 answers

No answers yet.

Log in to answer this question.

More posts like this