摘要: 参考地址: https://www.biostars.org/p/56246/ -q INT only include reads with mapping quality >= INT [0] Like for getting the unique reads (a single read mapping at one best position) 命令如下,也有文献报告使用的-q... 阅读全文
posted @ 2019-06-13 15:36 raisok 阅读(1555) 评论(0) 推荐(0) 编辑
摘要: 使用miniconda3进行安装 输入下面的代码进行测试 关闭warning信息使用下面的代码,不然会出现一些warning信息 阅读全文
posted @ 2019-06-13 08:40 raisok 阅读(2620) 评论(0) 推荐(0) 编辑
摘要: 输入数据格式 pathway = read.table("kegg.result",header=T,sep="\t") pp = ggplot(pathway,aes(richFactor,Pathway)) Pathwy是ID,richFactor是富集的基因数目除以背景的基因数目 改变点的大小 阅读全文
posted @ 2019-06-13 08:37 raisok 阅读(3624) 评论(0) 推荐(0) 编辑
摘要: 去掉某一列重复的行 index <- duplicated(a[,2]) a[!index,] 去掉重复的行用unique unique(a) 阅读全文
posted @ 2019-06-13 08:12 raisok 阅读(159) 评论(0) 推荐(0) 编辑