对结构变异长度大小图示的R脚本

if (!require(pacman)){
    install.packages("pacman")
    library(pacman)
}
### Load packages
pacman::p_load(ggplot2, ggthemes)
### import sv size file
sizebed=read.delim("XXX.bed",header=FALSE)
#bed文件第4列为结构变异分类（"DEL","INS","DUP","INV"），第5列为变异长度
colnames(sizebed)=c("chr","start","end","type","size","source")
sizebed$type=factor(sizebed$type,levels=c("DEL","INS","DUP","INV"))
pdf(file="XXX_svsize.pdf",width=6.5,height=5)
ggplot(sizebed, aes(x=type, y=size, color=type)) +labs(x="Type",y="Size")+scale_color_manual(name="Type",values=c("#007600","#1f1a4f","#E69F00","#999999"))+facet_grid(.~source)+
  geom_jitter(position=position_jitter(0.2),alpha=I(0.6))+theme_tufte(base_family="Helvetica")+theme(panel.border=element_blank(),panel.grid.minor=element_blank(),panel.grid.major=element_blank(),panel.background=element_blank(),axis.line.x=element_line("black"),axis.line.y=element_line("black"))+scale_y_continuous(limits=c(0,6000),breaks=c(0,500,1000,2000,3000,4000,5000,6000))
dev.off()