美文网首页
计算FRIP值

计算FRIP值

作者: 余绕 | 来源:发表于2021-09-28 17:17 被阅读0次

1.安装软件
软件安装在R4.1版本(激活conda的R4.1环境)

install.packages("remotes")
remotes::install_github("imbforge/encodeChIPqc")

出现如下选择时候最好都升级(升级时候如果个别安装不上先不管)

> remotes::install_github("imbforge/encodeChIPqc")
Downloading GitHub repo imbforge/encodeChIPqc@HEAD
These packages have more recent versions available.
It is recommended to update all of them.
Which would you like to update?

1: All                                 
2: CRAN packages only                  
3: None                                
4: gtools     (3.8.2  -> 3.9.2 ) [CRAN]
5: Rcpp       (1.0.6  -> 1.0.7 ) [CRAN]
6: data.table (1.14.0 -> 1.14.2) [CRAN]

Enter one or more numbers, or an empty line to skip updates: 1
gtools     (3.8.2  -> 3.9.2 ) [CRAN]
Rcpp       (1.0.6  -> 1.0.7 ) [CRAN]
data.table (1.14.0 -> 1.14.2) [CRAN]
spp        (NA     -> 1.16.0) [CRAN]

安装好后
加载包

library(encodeChIPqc)
library(GenomicRanges) #如果没有需要安装

进行计算

 peaks =  read.table('../../H4K5ac/WT/WT_1_peaks.narrowPeak')
 colnames(peaks) = c('chr','start','end','name','score','strand','signal','pval','qval','peak')
 peaks.gr = makeGRangesFromDataFrame(peaks,keep.extra.columns=TRUE)
 bam.file = '../../H4K5ac/WT/WT_1.bam' #sam转化后的bam
frip(bam.file,peaks.gr)

frip值

0.3997

encodeChIPqc好像只能在linux下用,其中有个依赖包只有Linux版本。

相关文章

网友评论

      本文标题:计算FRIP值

      本文链接:https://www.haomeiwen.com/subject/abdrnltx.html