RNA
An object of class Seurat
29875 features across 4224 samples within 2 assays
Active assay: integrated (2000 features, 2000 variable features)
2 layers present: data, scale.data
1 other assay present: RNA
2 dimensional reductions calculated: pca, umap
ATAC
An object of class Seurat
149666 features across 5309 samples within 2 assays
Active assay: peaks (110042 features, 110041 variable features)
2 layers present: counts, data
1 other assay present: RNA
2 dimensional reductions calculated: lsi, umap
but I get the following error messages:
Performing data integration using Seurat...
Performing log-normalization
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Calculating gene variances
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Calculating feature variances of standardized and clipped values
0% 10 20 30 40 50 60 70 80 90 100%
[----|----|----|----|----|----|----|----|----|----|
**************************************************|
Centering and scaling data matrix
|====================================================================================================| 100%
Running CCA
Merging objects
Finding neighborhoods
Finding anchors
Found 10350 anchors
Filtering anchors
Retained 2825 anchors
Warning: Please provide a matrix that has the same number of columns as the number of reference cells used in anchor finding.
Number of columns in provided matrix : 2976
Number of columns required : 4224
Skipping element 1.
Error: None of the provided refdata elements are valid.
In addition: Warning messages:
1: In LayerData.Assay5(object = assays[[i]], layer = lyr, fast = TRUE) :
multiple layers are identified by counts.1 counts.2
only the first layer is used
2: In LayerData.Assay5(object = object[[assay]], layer = layer, ...) :
multiple layers are identified by data.1 data.2
only the first layer is used
I think this error likely stems from some missing metadata processing step shown in the pre-processing script. However, I am unable to follow along with the data processing script because I am processing my ATAC data with Signac, not ArchR (I am working with plant datasets and have been unsuccessful generating the input genome objects required by ArchR). Do you think this is the most likely cause of this error - or am I missing some other crucial component. Thank you very much for your help, please let me know if you need additional info, I'll include my sessionInfo below.
Hi,
Thanks for maintaining a really nice R package. I am trying to use the CoembedData function as follows:
with the following Seurat objects:
but I get the following error messages:
I think this error likely stems from some missing metadata processing step shown in the pre-processing script. However, I am unable to follow along with the data processing script because I am processing my ATAC data with Signac, not ArchR (I am working with plant datasets and have been unsuccessful generating the input genome objects required by ArchR). Do you think this is the most likely cause of this error - or am I missing some other crucial component. Thank you very much for your help, please let me know if you need additional info, I'll include my sessionInfo below.