Hi everyone,
I am working with 2 datasets of single cell transcriptomics. I normalized and SCT transform the 2 datasets for the percentage of mitochondrial/ribosomal genes. I now want to integrate them and I tried to use PrepSCTIntegration on my seurat object containing my 2 groups with :
split_seurat <- PrepSCTIntegration(object.list = split_seurat, anchor.features = integ_features)
But I got the following error :
Warning: SCT model not present in assay Error in scale.data[anchor.features, ] : no 'dimnames' attribute for array
I don't get it as integ_features contained 2000 genes (features) obtained with this: integ_features <- SelectIntegrationFeatures(object.list = split_seurat, nfeatures = 2000)
Any idea on how to fix this issue ?
Thanks a lot for your help ! Best regards, Ludivine
