Comments (10)
Hi @mattcai10x ,
Thanks for sharing your findings! Unfortunately I won't be able to investigate the issue until the week of 04/08. Are you saying that, forget about SpotClean, the raw barcode matrix from SpaceRanger output would be pre-filtered to remove empty barcodes? This seems inconsistent with CellRanger where I remembered the raw barcode matrix usually has a fixed dimension (~737,280 barcodes?) Or are there recent updates in CellRanger and SpaceRanger where the empty barcodes are now filtered out?
from spotclean.
I actually did not know this but it looks like the change happened in CellRanger 6.0 and SpaceRanger 1.3.
from spotclean.
Hey!
I am facing the same issue. Upon reviewing my slide_info file, I noticed a discrepancy of 120 barcodes, which are different from the 14216 barcodes in my raw data. Is there an existing solution to address this problem, or is it sufficient to simply eliminate those inconsistent barcodes and proceed with the analysis?
Thank you!
from spotclean.
Hi @Streisenberg , thanks for reporting your issue. Are you seeing 120 barcodes in slide_info
but not in the raw matrix? And is your data coming from Visium platform? Not aware they are extending the number of spots to >14000.
from spotclean.
Hi @zijianni , thank you for your prompt attention to this matter. Yes, that is correct, I have 120 more barcodes in my slide_info
object.
And yes the data comes from the Visium CytAssist platform. I attached their 11x11mm slides here which contain >14000 barcodes within the capture area.
Thank you!
from spotclean.
Thanks for sharing! Yes, you are good to drop the 120 missing barcodes from slide_info
.
Note that the bigger Visium slide (~ 14000 spots) was not available at the time our paper was published, so you might want to validate the generated results before downstream analysis. Another note that the computation time and memory usage will increase with more spots in your data. The optimization process cannot easily be parallelized due to the iterative nature of EM algorithm.
from spotclean.
Thanks a lot (:
from spotclean.
I have the same issue in one of my data. Discard those not in the raw matrix can solve this.
m.raw <- read10xRaw(file.path(data.home, sample.name,"outs/raw_feature_bc_matrix"))
DEPRECATED.idx <- grep("DEPRECATED", rownames(m.raw))
m.raw <- m.raw[-DEPRECATED.idx,]
m.slideInfo <- read10xSlide(file.path(data.home, sample.name,"outs/spatial/tissue_positions.csv"),
file.path(data.home, sample.name,"outs/spatial/tissue_lowres_image.png"),
file.path(data.home, sample.name,"outs/spatial/scalefactors_json.json"))
if(sum(!m.slideInfo$slide$barcode %in% colnames(m.raw)) > 0){m.slideInfo$slide = m.slideInfo$slide[m.slideInfo$slide$barcode %in% colnames(m.raw),] }
m.obj <- createSlide(count_mat = m.raw, slide_info = m.slideInfo)
from spotclean.
Given that more people are facing the same issue, I will implement a fix in createSlide()
to remove spots in slide metadata but not in raw expression matrix.
from spotclean.
Addressed in #24 . After the PR merge, this feature is available in the dev package. The Bioconductor version needs to wait for the next release.
from spotclean.
Related Issues (20)
- install for R <= 4.1.3 HOT 3
- Adapt to Spaceranger 2.0 output HOT 1
- Error in optim(x_init, .fn_optim, .gr_optim, method = "L-BFGS-B", obs_exp = obs_exp, : non-finite value supplied by optim HOT 15
- SpatialExperiment::spatialData now deprecated HOT 1
- Questions about tissue regions HOT 2
- Questions about the output HOT 1
- How to modify strength of SpotClean? HOT 5
- could spotclean work under situation that all spot 's tissue==1? that my ST data 's tissue position csv 's tissue is all 1,no o,that no backgroud spot HOT 3
- The "max difference of decontaminated expressions" HOT 4
- How to run with large datasets in limited memory? HOT 2
- Error in createSlide(), 'barcodes in count matrix do not match barcodes in slide info' HOT 9
- The tissue photographs were presented in the wrong horizontal and vertical ratio HOT 2
- Error in convertToSeurat function, "Error in FUN(left, right) : non-numeric argument to binary operator" HOT 2
- Prepare for upcoming Seurat v5 release HOT 1
- SpotClean performance on 11x11mm slides and questions regarding to output HOT 1
- Error: non-finite value supplied by optim HOT 7
- Ensembl IDs to Gene Names
- Error in keepHighGene(raw_ts_data, verbose = verbose) : no slot of name "meta.data" for this object of class "Assay" HOT 5
- Spotclean() stalls on testdata HOT 2
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from spotclean.