Entering edit mode
6 months ago
carolofharvest
▴
50
Cell_Types <- c("Epi","T Cell","Myeloid","B Cell","Fibroblast","RBC","NK", "Endo","Acinar")
Epi_Markers <- c("KRT7","KRT8","KRT18","KRT19","EPCAM","CDH1")
T_Cell_Markers <- c("CD3E","CD3G","CD3D","CD4","IL7R","CD8A","LEF1")
Myeloid_Markers <- c("CD14","ITGAM","MNDA","MPEG1","ITGAX")
B_Cell_Markers <- c("CD79A","MS4A1","CD19")
Fibroblast_Markers <- c("CDH11","PDGFRA","PDGFRB","ACTA2")
RBC_Markers <- c("HBA1","HBB","HBA2")
NK_Markers <- c("NCR3","FCGR3A","NCAM1","KLRF1","KLRC1","CD38","KLRC1")
Endo_Markers <- c("CDH5","PECAM1")
Acinar_Markers <- c("TRY4","SPINK1","AMY2A")
All_Markers <- list(Epi_Markers,T_Cell_Markers,Myeloid_Markers,B_Cell_Markers,Fibroblast_Markers,RBC_Markers,NK_Markers,Endo_Markers,Acinar_Markers)
Epi_Score <- 0
T_Cell_Score <- 0
Myeloid_Score <- 0
B_Cell_Score <- 0
Fibroblast_Score <- 0
RBC_Score <- 0
NK_Score <- 0
Endo_Score <- 0
Acinar_Score <- 0
All_Scores <-list(Epi_Score,T_Cell_Score,Myeloid_Score,B_Cell_Score,Fibroblast_Score,RBC_Score,NK_Score,Endo_Score,Acinar_Score)
Weighted_Scores <- c()
Score_Weights <- c(1.85,1.85,2.22,3.7,2.78,3.7,1.85,5.56,3.7)
Is this ALL the context you can give us?
I've found this code in this site : . It doesn't seem like they've shared any additional information about this in the previous lines of code. I think these score weights are associated with the number of marker genes defined for each cell type, and they are inversely proportional. However, I couldnt figure out exactly how they are calculated