毕设正文和图片

毕设正文和图片（草稿）

计数

构建Seurat对象和数据预处理

构建Seurat对象：

> seu
An object of class Seurat 
93660 features across 181653 samples within 2 assays 
Active assay: RNA (78691 features, 0 variable features)
 1 layer present: counts
 1 other assay present: HERV
> colnames(seu@meta.data)
 [1] "orig.ident"               "nCount_RNA"               "nFeature_RNA"             "nCount_HERV"             
 [5] "nFeature_HERV"            "GSM"                      "group"                    "tissue"                  
 [9] "diagnosis"                "age"                      "sex"                      "PMD"                     
[13] "hist.diagnosis"           "diag.1"                   "diag.2"                   "diag.3"                  
[17] "age.related.plaque.score" "braak.tangle.stage"       "APOE"                     "percent_mito"            
[21] "percent_ribo"             "HERV_fraction" 
> table(seu$orig.ident)
  AD1  AD10  AD13  AD19   AD2  AD20  AD21   AD4   AD5   AD6   AD8   AD9  NC11  NC12  NC14  NC15  NC16  NC17  NC18 
 6165 16437  1841  3898 15607  9430  8736  5895 12623  3204  2679 10269  3560 15727 10951  8492  6899 14572 11821 
  NC3   NC7 
 6439  6408 

质控：

> scRNA_cells
An object of class Seurat 
35970 features across 164219 samples within 2 assays 
Active assay: RNA (21002 features, 0 variable features)
 1 layer present: counts
 1 other assay present: HERV
> table(scRNA_cells$group)
   NC    AD 
76567 87652 
> table(scRNA_cells$orig.ident)
  AD1  AD10  AD13  AD19   AD2  AD20  AD21   AD4   AD5   AD6   AD8   AD9  NC11  NC12  NC14  NC15  NC16  NC17  NC18 
 5865 15703  1517  3549 15080  9037  7691  5435 12140   230  2167  9238  1520 15185 10596  8175  3554 13900 11451 
  NC3   NC7 
 6056  6130 

过滤前：

过滤后：