• 最近做RNA-seq项目的时候准备用R的boxplot()工具画一个各个样品的箱线统计图。
  • 然而,在运行脚本后报错
    • Error in `[.data.frame`(fpkm, , c("MB7409-A", "MB7409-B", "MB7409-C", :
      undefined columns selected
      Calls: [ -> [.data.frame
      Execution halted
  • 抓取数据列的命令哪里错了呢?
  • 打开输入的数据框和脚本做对比
    • 数据是这样的
    • 获取数据列的脚本是这样的
  • 对需要比较的数据列输入应该是正确的,但是为什么会报错?别急,我们用R读取文件后再看看
    • 如图所示
  • 表头的样品名称里的符号“-”居然变成了“.”
  • 原来,特殊符号(如:-、(、?、%等)不能直接被read.table()识别,两个方法:
    • 一个是根据读取文本修改读取列
    • 第二个是可以在read.table()增加check.names=F,避免符号被修改
  • 于是,将脚本重新修改,脚本正常运行,得到箱线图
  • 技巧总结:
    • 1.遇到undefined columns selected 的问题一定是获取列的文本与数据列不对应
    • 2.眼见也不一定为实

      • 特殊符号需要敏感,尽量不要使用特殊符号命名样品;
      • 受兼容性影响,在windows系统生成的文本,在记事本中很多异常符号是显示不出来,可能在liunx系统里不支持(如空格,回车符等);
      • 2021.10.19补充:对于列名,纯数字输入默认会加X,如果设置check.names=F,尽管只显示数字,但仍然是字符串形式,而不是数字。因此,在调取列名的时候需要以字符串形式调用。
  • 欢迎进群讨论RNA-seq流程相关问题,遇见二维码过期可添加VX:bbplayer2021 ,备注 申请加入生信交流群。

2020.10.21【R语言】丨 undefined columns selected 问题解决办法相关推荐

  1. 关于STC8A8K64S4A12(STC15系列升级版)例程中程序出现 'xxx' undefined identifier的问题解决办法

    关于STC8A8K64S4A12(STC15系列升级版)例程中程序出现 'xxx' undefined identifier的问题解决办法 最近上网买了一款STC8A8K64S4A12最小系统版学习, ...

  2. 2020.10.21【转载】丨GWAS全基因组关联分析流程

    感谢CSDN用户 追梦生信人 梳理了GWAS全基因组关联分析的整个流程,并提供了基本的命令,用到的软件包括BWA.samtools.gatk.Plink.Admixture.Tassel等,在此分享出 ...

  3. R极简教程-10:R语言绘图基础

    R语言最强大的一点就是画图,那个无比强大的画图系统,再加上各种各样的神级R包,让我写R语言中最喜欢的一点就是写绘图代码. 简单来说,R语言原生的绘图系统已经非常强大了,根本不需要其他东西的辅助,就可以 ...

  4. 测试开发工作者日记:2020.10.21

    最近忙着搬家,西北旺那边村委会和某开发商 联合下达了腾退指令,所有人一周内滚蛋.作为每次抵挡的靠山-房东 也终于妥协了.周围很大一片区域的几乎半个海淀和半个昌平的公寓和房东们,都一瞬间知道了这个大新闻 ...

  5. R语言丨根据VCF文件自动填充对其变异位点并生成序列fa文件

    根据VCF文件自动填充对其变异位点并生成序列fa文件 首先提出一个问题: 假如有一个重测序结果VCF文件,里面包含了很多个样本在几百个突变位点(snp和iad)的基因型数据,现在想根据这份原始数据,得 ...

  6. R语言作图不显示中文解决办法,如何使用中文字体

    最近我们被客户要求撰写关于R语言作图的研究报告,包括一些图形和统计输出. 预期 R语言绘图中的中文默认字体宋体,希望用其他字体绘图,保存成PDF格式不能识别中文. 实际 使用"雅黑字体&qu ...

  7. R语言丨根据VCF文件设计引物,自动识别两样本差异SNP位点,调用samtools获取上下游参考序列

    根据变异位点设计引物序列 今天碰到一个新问题:假如有一个vcf文件储存了两个样品的变异位点基因型数据,每行代表一个位点,我现在想找出两样本差异的SNP位点,再把差异位点用[REF/ALT]的形式表示, ...

  8. Z变换(2020.10.21)

    1.Z变换 1.1 Z变换公式 1.2 matlab函数(ztrans) 1.2.1 Syntax(语法) ztrans(f)ztrans(f,transVar)ztrans(f,var,transV ...

  9. mdkstc系列器件支持包下载_WPS Office 2019 For Linux(2020/10/21)更新-支持PDF编辑

    WPS公共: 优化字体列表加载规则,解决同系列不同字形的东亚字体在Linux字体列表中加载不全的问题 优化二次开发接口机制,提升接口执行效率.绘制效率 完善JSAPI功能,支持更多接口调用 支持插入图 ...

最新文章

  1. Windows中配置java变量环境
  2. 转:Ubuntu中安装和配置 Java JDK,并卸载自带OpenJDK(以Ubuntu 14.04为例)
  3. volatile关键字(三)
  4. linux下的共享库(动态库)和静态库
  5. 使用async await 封装 axios
  6. Android应用程序安装过程源代码分析(1)
  7. ASP.NET提取多层嵌套json数据的方法
  8. mysql 5.5.18下载_MySQL5.7.18下载和安装过程图文详解
  9. python 列表索引位置,python – 在列表中查找值的第一个位置
  10. Mongodb 与sql 语句对照
  11. java的oracle事务回滚_Oracle事务处理
  12. 智能实验室-全能优化(Guardio) 4.0.0.691 beta 11
  13. 计算机辅助制造讲义翻译,计算机辅助制造讲义-2007-2演示文稿.PPT
  14. ArcGIS 字段计算器取前几位和替换操作
  15. 地理编码涉及的专项技术
  16. 【bzoj1614】[Usaco2007 Jan]Telephone Lines架设电话线 二分+SPFA
  17. java canvas 画线_java – 在扩展Canvas时使用drawLine()绘制粗线
  18. 指纹匹配算法matlab,指纹识别算法的matlab实现..doc
  19. 职场泥潭,这样的IT公司绝对不宜久留
  20. Apple Pencil平替哪个好?Apple Pencil平替笔推荐

热门文章

  1. MYSQL的操作语句(二)
  2. QQ农场启示录之:坐忘峰的盈利模式
  3. 更好的冠状病毒明信片
  4. 2022以后云服务器的发展趋势
  5. 【视觉SLAM十四讲】第一讲 概述与预备知识
  6. 二元函数对xy同时求导_高等数学期末总复习 DAY4. 利用莱布尼茨定理求高阶导 隐函数求导 对数求导法 参数函数求导等...
  7. 高品质的音乐蓝牙耳机有哪些推荐?音质好的蓝牙耳机排行榜
  8. Java SE 流程控制 Struct
  9. 手机扫码登录电脑QQ出现无法登录问题
  10. 第五课 田中さんは会社へ行きます。