library(ggplot2)
data <- read.table("C:\\Users\\86150\\Desktop\\fragments\\RD20220208-LPC-12462FND.fragment.length.txt",header=F)
length_cutoff <- 1200
fragment <- data$V1[data$V1 <= length_cutoff]
breaks_sum <- 500
res <- hist(fragment,breaks = breaks_sum,plot = FALSE)
plot(x=c(0,res$breaks),y=c(0,0,res$counts)/10^3,type="l",col="red",xlab = "Fragment length(bp)",ylab = expression(Normalized ~ read ~ density ~ 10^3),main="Sample Fragment sizes")
结果图

网友评论