对特殊字符“.”、“+”这些要特别注意,很有可能影响你的指令!

用gatk来call体细胞突变时bam文件样本名始终报错

我在用gatk mutect2 call somatic variants,报错bam样本名。用samtools view -h xxx.bam看过了bam header没有问题。但还是一直报错。找了很久,最后把bam文件名的“.”去掉就可以了。小数点在代码里常常有特殊意义,最好不要用到参数上!

用clusterProfiler包的bitr转换gene id时识别不了输入的ENSEMBL

从TCGA下载的数据中,样本的"ENSEMBL"带有版本号,如ENSG00000000003.13,表示经过多次修正可能现在是第13个版本。而基因集中的‘ENSEMBL’肯定是以ENSG00000000003存在的。可以用字符串替代函数gsub将小数点及其后的数字们都替换成空,即删掉。注意到gsub()函数中一些特殊符号,如小数点"."、加号"+“等要加上两个反斜杠”\"来识别。

library(org.Hs.eg.db)
library(clusterProfiler)
gene <- COAD_counts$ECSG_ID
tmp=gsub("\\..*","",gene)
gene_id <- bitr(tmp, fromType = 'ENSEMBL', toType = c("SYMBOL", "GENENAME"), OrgDb = 'org.Hs.eg.db')

特殊符号“.”对命令识别的影响:bam样本名报错 和 转换gene id时‘ENSEMBL’ keys无法识别相关推荐

  1. 【错误记录】记录 Android 命令行执行 Java 程序中出现的错误 ( dx 打包 PC 可执行文件报错 | dalvik 命令执行 kotlin 编译的 dex 文件报错 )

    文章目录 前言 一.Android 命令行与 PC 可执行 JAR 文件不兼容 二.Android 命令行使用 dalvik 命令不能直接执行 Kotlin 编译的 dex 文件 前言 尝试在 And ...

  2. 执行redis命令redis-trib.rb查看集群信息报错cannot load such file -- redis (LoadError)

    问题描述: 在执行redis-trib.rb命令查看集群状态的时候,报错: [aiprd@hadoop1 ~]$ redis-trib.rb check 192.168.30.10:7000 Trac ...

  3. 怎么解决在打开Visual Studio时的报错“无法找到一个或多个组件”以及关于Visual Studio Installer无法识别已安装的VS并报错找不到路径的问题

    第一个问题:打开Visual Studio时的报错"无法找到一个或多个组件". 情况1:如果是点击安装后桌面生产的Visual Studio的快捷方式而导致出错,建议定位到Visu ...

  4. tar命令解压jdk.tar.gz包 报错 gzip: stdin: not in gzip format

    转自:https://blog.csdn.net/LL_zhuo/article/details/44173355 遇到和这篇博文一样的问题了.用wget 从oracle官网下载jdk, http:/ ...

  5. [debug] 命令行窗口运行.py文件,报错No module named XXX

    问题描述: 在PyCharm中运行.py文件正常,在命令行运行.py文件报错"No module named XXX" 原因分析 说到底是搜索路径的问题. PyCharm中运行的时 ...

  6. 用pytorch官网命令 安装pytorch1.10.1+CUDA11.1报错

    因为服务器的cuda是11.1的所以需要安装历史版本 复制官网的命令 CUDA 11.1 pip install torch==1.10.1+cu111 torchvision==0.11.2+cu1 ...

  7. git 命令行配置及配置文件 解决clone报错 LibreSSL SSL_connect: SSL_ERROR_SYSCALL in connection to github.com:443

    Intro clone 某个来自于 github 的仓库的时候,报错如下: leung@wuyujin coderepo % git clone https://github.com/spring-g ...

  8. 香橙派 ubuntu c# .net 使用ThoughtWorks.QRCodCore库生成 或识别 二维码,及报错处理方法

    在linux下生成和识别二维码,试了大半天,找到个库测试通过可以使用 ThoughtWorks.QRCodCore 其它的库有些仅在window下可用,有些使用困难 这里如果你在linux下跑的话,一 ...

  9. ROS使用tensorflow识别物体、以及出现的报错。

    物体识别使用了tensorflow_object_detector的功能包,具体实现可以参考以下文章,在此不再赘述. https://blog.csdn.net/qq_37464350/article ...

  10. windows系统,使用命令行启动MySQL服务失败,报错:服务名无效。

    方式一:命令行启动 在windows系统中,使用net start mysql启动MySQL服务,报错:服务名无效. 问题原因: 服务器名错了 解决办法 :应该使用命令net start mysql8 ...

最新文章

  1. 网络工程制图论文计算机,学习系统与工程制图论文
  2. php记录网站访问,PHP简单实现记录网站访问量的功能
  3. c语言无法打开源文件xx.h,VS2015 + Qt5.9.2开发中无法打开源文件“ui_*.h” 和 error MSB6006: “cmd.exe”已退出,代码为 3之解决办法。...
  4. jQuery数据转换与提交
  5. java数组的实现_Java数组HashCode实现
  6. java迭代器输入的结果_Java学习之Iterator(迭代器)的一般用法 (转)
  7. python 参数收集_Python参数传递及收集机制原理解析
  8. 【实践】CTR预估在动态样式建模和特征表达学习方面的进展
  9. java word2007_Java解析word2007、Excel2003和Excel2007
  10. Android手机证书安装
  11. 迪士尼超级IP版图日趋完整
  12. Alexa 开发指南
  13. 玩~成语接龙c++代码
  14. HP 8440P 4G内存蓝屏--内存颗粒容量支持问题
  15. 笔记本无线自动配置服务器,笔记本无线网络连接IP地址怎么设置
  16. 戏精程序员凌晨三点半发周报,真的很做作!
  17. 【在CentOS中使用Crow-基于C++的Web服务】
  18. 大学计算思维c语言答案,【程序设计论文】计算思维为下的C语言程序设计教学(共2755字)...
  19. Odoo开发应该怎么学习?
  20. 微信小程序textArea输入框随着输入字数自适应高度

热门文章

  1. 网易微专业——Java Web开发工程师学习笔记(2):Tomcat
  2. ping服务器地址显示TTL传输中过期,ping公网固定IP出现: TTL 传输中过期,怎么处理?...
  3. 计算机病毒中的后门病毒,后门病毒是什么?
  4. party_bid_core三种数据结构总结
  5. linux输入文件后clustalw,合并提取后的domain序列之后,linux系统中的clustalw不能读出蛋白信息...
  6. Vulkan 多线程渲染
  7. 多传感器数据融合简介(转)
  8. Steam 最新VR 游戏畅销榜——这几款游戏绝对值得一玩~不来看看你可就错过啦!
  9. 【RL】同策略(on-policy)与异策略(off-policy)
  10. 详解:分治算法【Java实现】——汉诺塔问题