b.limour.top/source/_posts/格里森评分数据处理.md

---

title: 格里森评分数据处理 tags: [] id: '1958' categories:

• • uncategorized

    date: 2022-08-09 22:27:03

    f_TCGA_gleason_grade <- function(primary_gleason_grade, secondary_gleason_grade){
    primary_gleason_grade <- as.numeric(unlist(data.frame(strsplit(primary_gleason_grade, ' '))[2,]))
    secondary_gleason_grade <- as.numeric(unlist(data.frame(strsplit(secondary_gleason_grade, ' '))[2,]))
    primary_gleason_grade + secondary_gleason_grade
    }
    clinical <- readRDS('../../../zl_liu_new/Bulk/tcga/PRAD/clinical.rds')
    tmp <- clinical[c('psa', 'age_at_initial_pathologic_diagnosis', 'ajcc_pathologic_t', 'ajcc_pathologic_n', 'ajcc_clinical_m')]
    names(tmp) <- c('PSA', 'Age', 'T', 'N', 'M')
    tmp$GS <- f_TCGA_gleason_grade(clinical$primary_gleason_grade, clinical$secondary_gleason_grade)
    tmp$T <- factor(tmp$T, ordered = T)
    tmp$N <- factor(tmp$N, ordered = T)
    tmp$M <- factor(tmp$M, ordered = T)