search

smorabit / hdWGCNA

High dimensional weighted gene co-expression network analysis
https://smorabit.github.io/hdWGCNA/
Other
346 stars 34 forks source link

RNA assay in seurat object for hdWGCNA, should l removeintegrated assay? #309

Open yulchen810 opened 3 weeks ago

yulchen810 commented 3 weeks ago

Hi, I have a seurat object containing 6 samples(seurat v4), I performed integration to correct batch effect-following https://satijalab.org/seurat/articles/integration_introduction.html. Now I have two assays: "RNA" and "integrated".

I want to use this seurat object to do hdWGCNA analysis.

  1. Should I use "RNA" assay? Because I tried to use "integrated" assay, an error occur when SetupForWGCNA, and it worked when I use "RNA" assay.
  2. Should I remove "integrated" assay? I tried to do MetacellsByGroups after the fisrt step in 1, it return me an ERROR again, Error in MetacellsByGroups(seurat_obj = seurat_obj, group.by = c("cell_type", : Invalid reduction (harmony). Reductions in Seurat object: pca, umap

The structure of my seurat object is like below: image

The data in "RNA" assay only be normalized, so I need to perform ScaleData, PCA, Harmony and UMAP to "RNA" assay, and use the result to do hdWGCNA, right? If then, do I need to remove integrated assay from seurat object or I just keep it which will not affect hdWGCNA analysis?

Thank you very much! Best

smorabit commented 1 day ago

I think you should use the RNA assay because it seems like your integrated assay only contains 2000 features.