hahahueg 2023-03-01 12:02 采纳率: 87.5%
浏览 118
已结题

tcga预后森林图绘制

使用coxph函数绘制森林图,警告相关性无穷,最后报错不知道是什么东西,不知道是不是数据问题


 risk=read.table(riskFile, header=T, sep="\t", check.names=F, row.names=1)
> cli=read.table(cliFile, header=T, sep="\t", check.names=F, row.names=1)
> cli=cli[apply(cli,1,function(x)any(is.na(match('unknow',x)))),,drop=F]
> cli$Age=as.numeric(cli$Age)
> sameSample=intersect(row.names(cli), row.names(risk))
> risk=risk[sameSample,]
> cli=cli[sameSample,]
> rt=cbind(futime=risk[,1], fustat=risk[,2], cli, riskScore=risk[,(ncol(risk)-1)])
> multiCox=coxph(Surv(futime, fustat) ~ ., data = rt)
Warning message:
In coxph.fit(X, Y, istrat, offset, init, control, weights = weights,  :
  Loglik converged before variable  7 ; coefficient may be infinite. 
> multiCoxSum=summary(multiCox)
> multiTab=data.frame()
> multiTab=cbind(
+                  HR=multiCoxSum$conf.int[,"exp(coef)"],
+                  HR.95L=multiCoxSum$conf.int[,"lower .95"],
+                  HR.95H=multiCoxSum$conf.int[,"upper .95"],
+                  pvalue=multiCoxSum$coefficients[,"Pr(>|z|)"])
> multiTab=cbind(id=row.names(multiTab),multiTab)
> write.table(multiTab, file="cox.result.txt", sep="\t", row.names=F, quote=F)
> #????ɭ??ͼ
> pdf(file="forest.pdf", width=8, height=6, onefile = FALSE)
> ggforest(multiCox,
+          main = "Hazard ratio",
+          cpositions = c(0.02, 0.22, 0.4), 
+          fontsize = 0.8, 
+          refLabel = "reference", 
+          noDigits = 3)
Error in axisTicks(rangeb/2, log = TRUE, nint = 7) : 
  log - axis(), 'at' creation, _LARGE_ range: invalid {xy}axp or par; nint=7
     axp[0:1]=(1e-307,1e+308), usr[0:1]=(0,inf); i=615, ni=87
In addition: Warning message:
In .get_data(model, data = data) :
  The `data` argument is not provided. Data will be extracted from model fit.

  • 写回答

2条回答 默认 最新

  • 社区专家-Monster-XH 2023-03-01 12:09
    关注

    基于Monster 组和GPT的调写:
    看起来是绘图时的一个坐标轴问题。由于 Cox 回归模型中某个变量可能具有无限的系数,因此在绘图过程中出现了这个问题。此外,在警告信息中,coxph 函数也提示您在拟合模型时出现了收敛问题。这些问题可能是由于数据问题、模型选择或参数设置等多种因素导致的。

    本回答被题主选为最佳回答 , 对您是否有帮助呢?
    评论
查看更多回答(1条)

报告相同问题?

问题事件

  • 系统已结题 3月11日
  • 已采纳回答 3月3日
  • 创建了问题 3月1日

悬赏问题

  • ¥50 深度学习运行代码直接中断
  • ¥15 关于#单片机#的问题,请各位专家解答!
  • ¥15 关于#单片机#的问题,请各位专家解答!
  • ¥20 需要完整的共散射点成像代码
  • ¥15 编写vba代码实现数据录入工作
  • ¥15 做过TCL海信电视小米电视相关影视会员软件私我
  • ¥15 Mapreduce是正常的,在运行其他jar包时并没有任何问题,只是在做LogCount.jar 时出的问题。如图所示
  • ¥15 ImportError: DLL load failed while importing _iterative: 找不到指定的模块。
  • ¥15 如何通过交互分析得出某高危患者对放疗获益更多
  • ¥15 相关性分析中,p<0.05, r=0.29,怎么评价相关性呢