自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+
  • 博客(3)
  • 收藏
  • 关注

原创 只有一列的数据框删除特定行后如何还能保持原有数据框的行列结构

只有一列的数据框在删除某些行后,会自动转换成numeric类型,可能只剩一行向量了。如何保持原有的数据框行列名和结构呢?加一列重复值,一劳永逸以下展示一个简单的数据框,我们删除尾号为01A的数据。可以看到删除改行后变成了类型为numeric的向量countscounts_normal <- counts[-grep("01A", rownames(counts)),]counts_normal结果如下我们可以添加一模一样的一列,再进行操作。问题可能是以后一直要带着这列“repeat”

2020-12-17 14:52:59 687

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

对特殊字符“.”、“+”这些要特别注意,很有可能影响你的指令!用gatk来call体细胞突变时bam文件样本名始终报错我在用gatk mutect2 call somatic variants,报错bam样本名。用samtools view -h xxx.bam看过了bam header没有问题。但还是一直报错。找了很久,最后把bam文件名的“.”去掉就可以了。小数点在代码里常常有特殊意义,最好不要用到参数上!用clusterProfiler包的bitr转换gene id时识别不了输入的ENSEMB

2020-12-09 21:00:03 8065

原创 根据新旧文件名一一对应的表格批量循环修改文件名

Linux根据新旧文件名一一对应的表格批量循环修改文件名从TCGA上下载的①sample sheet文件包含file name,file ID 和对应sample type②gdc download文件含有counts和FPKM的两个文件③clinical文件:含有case submitter id和一些临床数据。这里的case id和file id不同根据gdc download文件中manifest文件中的file ID(如0096f6b3-94f5-49d9-aa43-66d9ab6b2b5

2020-12-08 17:19:14 1039

空空如也

空空如也

TA创建的收藏夹 TA关注的收藏夹

TA关注的人

提示
确定要删除当前文章?
取消 删除