123anjuan / HMA

Human Muscle Ageing Cell Atlas
9 stars 2 forks source link

about:block 1/03.figure1d_LTSR.R #1

Open kaqisekuzi opened 2 months ago

kaqisekuzi commented 2 months ago

@iramirezpardo @123anjuan Hello,
I saw "patient_id_sample" in the code below you. May I ask if "patient_id_sample" is the sample information (such as when I sequenced and measured 6 samples, 3 blanks, and 3 treatments)? Or does it refer to "sample ID+cell ID"? Thank you very much. Looking forward to your reply.

meta = read.csv(args[1],row.names=1) #the matedata of merged snRNA,scRNA and snATAC GeneScore RDS Y = table(meta$patient_id_sample,meta$Final_annotation)

123anjuan commented 2 months ago

Hi,I apologize for any confusion. The "patient id sample" should refer to single-cell libraries for each omics technology, typically corresponding to the "orig.ident" of a seurat object. For example, if you have sequenced and measured 6 samples (with 3 libraries for each individual sample), then the total number of "patient id sample" should be 18... Please let me know if you need further clarification.