kaks值分析可视化
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()
页:
[1]