- 博客(9)
- 收藏
- 关注
原创 ggplot+geom_col
library(tidyverse)df<-read.csv('ref_thinTae_wheatH03H04.csv',sep = ',',header = T)head(df)ggplot(df)+geom_col(aes(gsub("Chr","",chr),reas.count.ratio,fill=type),position="dodge",width=0.6,alpha=0.75)+ geom_hline(aes(yintercept=1.5), colour="#99000.
2022-03-17 15:48:28 835
原创 count计数
library(tidyverse)df<-read.table('ID.txt',sep='\t',header=F)head(df)df1<-df%>%count(V1)#df2<-data.frame(df1)head(df2)write.table(df2,file="idcount1.txt")
2022-01-13 16:14:17 226
原创 mosdepth计算depth
##make pos.bedbedtools makewindows -g CS_genome.txt -w 5000000 -s 2500000 > CS_win.bedCS_genome.txtCS_win.bed###calculate the depth of region mosdepth -b CS_win.bed H3_output sample1.Aligned.sortedByCoord.out.bam &###resultPer-..
2021-12-22 10:39:24 646
原创 ggplot+geom_bar
data type:library(tidyverse)df <-read.csv("spikenum.csv",header = T,sep = ',')head(df)ggplot(df,aes(x=x,y=heritability,fill=variable))+geom_bar(position="dodge",stat="identity")+ theme_classic()+ theme( axis.text = element_text(size=12),
2021-12-14 22:31:58 1909
原创 histogram
df <- read.table("kenelnum.txt",header = F,sep="\t")ggplot(df,aes(V2,fill=V1))+geom_histogram()+theme_bw()+theme(axis.text = element_text(size=12),axis.title=element_text(size=12),plot.margin = unit(c(1,1,1,1),"cm"),legend.title = element_blank...
2021-12-14 20:36:19 226
原创 Linux安装anaconda3
1、打开anaconda网址,复制linux版本链接,在Linux系统下用命令行下载:wget -c https://repo.anaconda.com/archive/Anaconda3-2021.05-Linux-x86_64.sh2、bash安装bash Anaconda3-2021.05-Linux-x86_64.sh3、配置环境vim ~/.bashrc##添加环境export PATH=/public-supool/home/zhangzm/anaconda3/bin
2021-07-08 14:20:25 353
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人