怎么用R語(yǔ)言統(tǒng)計(jì)相同ID出現(xiàn)的頻率

本篇內(nèi)容介紹了“怎么用R語(yǔ)言統(tǒng)計(jì)相同ID出現(xiàn)的頻率”的有關(guān)知識(shí)，在實(shí)際案例的操作過(guò)程中，不少人都會(huì)遇到這樣的困境，接下來(lái)就讓小編帶領(lǐng)大家學(xué)習(xí)一下如何處理這些情況吧！希望大家仔細(xì)閱讀，能夠?qū)W有所成！

R語(yǔ)言代碼，統(tǒng)計(jì)向量的頻次，其實(shí)很簡(jiǎn)單，用table方法一步就完成，后面就是繪圖了：

注意本代碼的主題為cowplot主題適合SCI文章發(fā)表，并設(shè)置了柱狀圖上加數(shù)字。

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
}
#統(tǒng)計(jì)頻率，并把大于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語(yǔ)言統(tǒng)計(jì)相同ID出現(xiàn)的頻率”的內(nèi)容就介紹到這里了，感謝大家的閱讀。如果想了解更多行業(yè)相關(guān)的知識(shí)可以關(guān)注億速云網(wǎng)站，小編將為大家輸出更多高質(zhì)量的實(shí)用文章！