kaks值分析可视化

yy951010

1.整理所需要的KaKs.xls表格；
2.运行相关R包背景，并开始计算。




rm(list=ls())
library(ggplot2)
library(reshape2)
data <- read.table("all.1.ks",header = T,sep="\t")
data = na.omit(data)
head(data)
p2 <- ggplot(data,aes(KaKs,color = Species, fill = Species))  +
  geom_line(stat="density",alpha=0.8)  + xlim(0,3)  + xlab("KaKs") +theme_classic()  +
  theme(axis.title = element_text(size=16),axis.text=element_text(size=16))
pdf(file=paste("newKaks_density.pdf", sep = "."), width = 10, height = 7)
p2
dev.off()