Description
工作流程发布,不足交流探讨及改进
工作流程发布,不足交流探讨及改进
准备文件:geneID
三个碱基翻译成一个蛋白质,但是如何检查CDS文件中几万个基因中,哪些CDS的序列是三的倍数?
存在问题:计算红豆杉中基因对之间的Ka.Ks的值全都是NaN
#柱状图
#=========================================================================
# 加载包
#=========================================================================
library(ggplot2)
library(dplyr)
library(reshape2)
#=========================================================================
# 1.library(ggplot2)
library(dplyr)
library(reshape2)
# ===================== 数据 =====================
eggnog <- data.frame(
Category = c("A","B","C","D","E","F","G","H","I","J","K","L","M","N","O","P","Q","R","S","T","U","V","W","Y","Z"),
disease = c(1331.92,583.42,59792.40,12647.23,72533.57,17147.12,
58640.01,31475.03,36986.05,37400.36,58142.47,67299.68,
55828.29,9119.68,32528.88,45813.37,34376.89,0,
146233.98,50492.20,17621.21,17423.12,279.45,28.24,766.37),
healthy = c(1979.61,1166.96,61561.81,11042.在操作时发现,由于个人习惯喜欢在tmux中操作,所以有时候停止进行,但是会产生临时文件
临时文件不可以直接rm删除
首先,查看文件的进程
lsof .nfs000000003074ec5a0f349c35