Closed cristalliao closed 1 year ago
Hi Cristal,
Most information in the AnnData object comes from file "GSE147747_meta_table.tsv", which is provided by the original dataset. To perform the alignment of slices, we added the 2D location information (an N-by-2 matrix) as .obsm['loc_use']. The aligned 2D locations are in .obsm['spatial_aligned'], while 3D locations are in .obsm['3D_coor']. The 3D spatial location neighbourhood graph is .obsm['graph']. We also saved learned representations as .obsm['latent'], which can be treated as a dimension reduction result when using scanpy for analysis. Basically, you can still use scanpy to perform other standard analysis like preprocessing, visualization, finding marker genes, etc..
Dear Geifei, Thanks for your explanations. It is really helpful!! Also, I still have some questions regarding this result:
Thanks a lot! Best regards, Cristal
Hi Cristal,
Best, Gefei
Dear Geifei,
Thanks for your explanations. It is really helpful!! Also, I want to know where I can find the gene expression information. Which variable is gene expression data?
AnnData object with n_obs × n_vars = 1480 × 3817 obs: 'nGene', 'nUMI', 'Sample', 'weeks', 'ChipBatch', 'ChipNr', 'Experiment_date', 'Experiment_procedure', 'Sequencing_date', 'Raw_reads', 'new_x', 'new_y', 'percent.mito', 'res.0.8', 'selected', 'array_row', 'array_col', 'slice', 'batch', 'library_size', 'n_genes' var: 'n_cells' uns: 'log1p' obsm: 'spatial', 'loc_use', 'spatial_aligned', 'count', 'graph', '3D_coor', 'latent'
Does the gene expression data change after training the STitch3D model? Where is the gene expression data before the model training?
Moreover, I want to know whether the x, y, and z axis(3D information) changed due to the STitch3D model training.
Thanks a lot! Best regards, Cristal
Hi Cristal,
You can get gene expression in AnnData objects using "adata_st.X". Check here for more information about the data structure in AnnData. In STitch3D, we saved normalized and log-tansformed gene expression data in "adata_st.X". You can get raw data from the original datasets. Also, gene expression data and location information are not changed after STitch3D model training.
Best, Gefei
Dear Professors, I am interested in performing some "scanpy" analysis on the results using the STitch3D model in Example 2(6PCW human heart dataset), which is "model.adata_st". This object is Anndata, so I want to perform some "scanpy" analysis. Do you have any suggestions for me to do this? Also, could I get some information on the different meanings of variables?
AnnData object with n_obs × n_vars = 1480 × 3817 obs: 'nGene', 'nUMI', 'Sample', 'weeks', 'ChipBatch', 'ChipNr', 'Experiment_date', 'Experiment_procedure', 'Sequencing_date', 'Raw_reads', 'new_x', 'new_y', 'percent.mito', 'res.0.8', 'selected', 'array_row', 'array_col', 'slice', 'batch', 'library_size', 'n_genes' var: 'n_cells' uns: 'log1p' obsm: 'spatial', 'loc_use', 'spatial_aligned', 'count', 'graph', '3D_coor', 'latent'
I am a little confused about each variable in 'obs', 'var', 'uns' and 'obsm'. Could you provide me a dictionary to know with these meanings? Thanks a lot!