- 博客(7)
- 收藏
- 关注
原创 用survminer优雅的画生存曲线
生存曲线绘制#Survival Curves#The ggsurvplot() function creates ggplot2 plots from survfit objects.data(lung)head(lung) inst time status age sex ph.ecog ph.karno pat.karno meal.cal wt.loss1 3 3...
2018-09-08 21:24:44 14707 4
原创 单因素cox
outsur <- coxph(Surv(Time,Status)~age, data=hcc_risk_clinical_validate)forest_model(outsur,factor_separate_line=F, format_options = list(colour= "black", shape = 16, ...
2018-08-12 12:34:09 1932
原创 森林图
d &lt;- c(NA,NA,0.3,0.3,0.6,rep(c(NA,0.3,0.3,0.6),8))forestplot(as.matrix(forestplot.OS[,c(1,2,3,7,8)]),mean=forestplot.OS$V4,lower = forestplot.OS$V5, upper = forestplot.OS$V6,graph.pos =...
2018-08-12 12:31:21 1553
原创 TCGA-HCC-差异分析-富集分析-pathview
library(pathview)library(TCGAbiolinks)library(SummarizedExperiment)library(clusterProfiler)#下载HCC的Gene expression文件query &lt;- GDCquery(project = "TCGA-LIHC", data.category = "G...
2018-08-09 16:24:12 3917
原创 ROC曲线
rm(list=ls())library(dplyr)library(survivalROC)library(rbsurv)library(parallel)#detectCores()检查当前电脑可用核数cl.cores &lt;- detectCores()cl &lt;- makeCluster(24)registerDoParallel(cl)dim(cox_hcc_me...
2018-08-08 12:28:05 1412
原创 TCGAbiolinks下载临床数据报错
query <- GDCquery(project = "TCGA-GBM", data.category = "Clinical", barcode = c("TCGA-08-0516","TCGA-02-0317")) GDCdownload(query)clinical <- GDCprepare_cli...
2018-08-03 13:46:53 3038
原创 R包安装报错
ERROR: failed to lock directory ‘J:/software/R-3.5.1/library’ for modifying#添加参数:INSTALL_opts = c('--no-lock')install.packages("C:/Users/Administrator/AppData/Local/Temp/RtmpQJeYze/downloaded_pac...
2018-08-01 11:27:08 1883 1
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人