本篇内容介绍了“怎么用r语言统计相同id出现的频率”的有关知识，在实际案例的操作过程中，不少人都会遇到这样的困境，接下来就让小编带领大家学习一下如何处理这些情况吧！希望大家仔细阅读，能够学有所成！
r语言代码，统计向量的频次，其实很简单，用table方法一步就完成，后面就是绘图了：
注意本代码的主题为cowplot主题适合sci文章发表，并设置了柱状图上加数字。
library(reshape2)
local({r <- getoption("repos") ;r["cran"] <- "http://mirrors.tuna.tsinghua.edu.cn/cran/" ;options(repos=r)})
library(ggplot2)
library(cowplot)
library(rcolorbrewer)
pairjoin <- function(x){
ran=x
ran[length(ran)]=">=11"
ran=ran[-1]
ran
}
#统计频率，并把大于11的 归类到一起：
data=table(data)
data[data>10]=11
massstatn <- with(hist(data, breaks=seq(0, 11, by = 1), plot=false),
data.frame(n=counts, mass=pairjoin(breaks), pct=counts/sum(counts)))
massstatn
massstatn$mass=factor(massstatn$mass,levels = massstatn$mass,order=t)
pn=ggplot(data=massstatn, aes(x=mass, y=n))
geom_bar(fill="#4daf4a",alpha = .9, stat="identity",width=0.8)
geom_text(aes(x=mass,y=n 20,label=n))
guides(fill=false)
theme(legend.key = element_blank(),legend.title = element_blank()
) xlab("peptide number") ylab("protein number") ggtitle("peptide number distribution")
pn
“怎么用r语言统计相同id出现的频率”的内容就介绍到这里了，感谢大家的阅读。如果想了解更多行业相关的知识可以关注本站网站，小编将为大家输出更多高质量的实用文章！
