相关性分析作图(热图+数值)/(直方图+数值)R语言
摘要:library('corrplot') library(xts) library(zoo) library(PerformanceAnalytics) setwd("F:/BaiduNetdiskWorkspace/") #读取数据 a<-readxl::read_xlsx("phe.xlsx",s
阅读全文
交互或批量计算区间内基因数目(while read)
摘要:###1.交互状态 #!/bin/bash #count genes in region,need three para #author lee if [ $# -ne 3 ] then echo "Usage chr start end" else num=$(awk -v chr=$1 -v s
阅读全文
曼哈顿图绘制-qqman-CMplot-ggplot(强推)-多个性状-区间放大-批量注释(更改SNP名称为基因名,或添加新列进行自定义注释)
摘要:if(!require("qqman"))install.packages("qqman") if(!require("CMplot"))install.packages("CMplot") library(qqman) #library(CMplot) library(tidyverse) set
阅读全文