用maftools分析GISTIC文件报错

好吃懒做的大白猫 2020-05-28 01:37:59
我用了两种方法做GISTIC,分别在不同的地方报错,有没有大神可以帮忙解决一下

方法一:从TCGA官网下载文件后在GISTIC在线平台分析
在读取GISTIC文件处报错

代码如下:
rm(list = ls())
options(stringsAsFactors = F)
library(maftools)

laml.gistic = readGistic(gisticAllLesionsFile = 'C:/Users/Meredith/Desktop/maftools/all_lesions.conf_90.txt',
gisticAmpGenesFile = 'C:/Users/Meredith/Desktop/maftools/amp_genes.conf_90.txt',
gisticDelGenesFile = 'C:/Users/Meredith/Desktop/maftools/del_genes.conf_90.txt',
gisticScoresFile = 'C:/Users/Meredith/Desktop/maftools/scores.gistic',
isTCGA = TRUE)

#结果报错
>Error in `levels<-`(`*tmp*`, value = as.character(levels)) :
factor level [61] is duplicated

方法二:在firehose直接下载GISTIC分析结果

在生成图片处报错

代码如下:
plotCBSsegments(cbsFile = 'C:/Users/Meredith/Desktop/maftools/BLCA.snp__genome_wide_snp_6__broad_mit_edu__Level_3__segmented_scna_minus_germline_cnv_hg19__seg.seg.txt' )

#结果报错
>Error in rep(chr.colors, seg.spl.transformed[, .N, Chromosome][, N]) :
'times'参数不对
...全文
635 2 打赏 收藏 转发到动态 举报
写回复
用AI写文章
2 条回复
切换为时间正序
请发表友善的回复…
发表回复
  • 打赏
  • 举报
回复
下载最新的R包就可以了
makoudada 2021-07-23
  • 举报
回复
@好吃懒做的大白猫 是因为下载的TCGA数据的版本和maftools与冲突么

3,423

社区成员

发帖
与我相关
我的任务
社区描述
其他开发语言 其他开发语言
社区管理员
  • 其他开发语言社区
加入社区
  • 近7日
  • 近30日
  • 至今
社区公告
暂无公告

试试用AI创作助手写篇文章吧