## ----------------------------------------------------------------------------- 库(ExperimentHub)图书馆(SummarizedExperiment)库(“ggplot2 ") ## ----------------------------------------------------------------------------- 嗯= ExperimentHub()查询(呃,“CLLmethylation”)冰毒=呃[[“EH1071”]]#提取甲基化数据 ## ----------------------------------------------------------------------------- methData = t(化验(冰毒))#过滤器只包括大多数变量网站5000强ntop = 5000 methData = methData[、订单(应用(methData 2 var, na.rm = TRUE),减少= TRUE) [1: ntop ]] ## ----------------------------------------------------------------------------- # 主成分分析pcaMeth = prcomp (methData中心= TRUE,规模。= FALSE ) ## ----------------------------------------------------------------------------- 总结(pcaMeth ) ## ----------------------------------------------------------------------------- tmp = data.frame (pcaMeth $ x) ggplot (data = tmp, aes (x = PC1, y = PC2)) + geom_point () + theme_bw () ## ----------------------------------------------------------------------------- sessionInfo ()