我通过命令行安装了“ caTools” R软件包:
$ R
$ install.packages("caTools", lib="~/R/library")
然后,执行以下命令:
INPUT=/home/user/file.bam
OUTPUT=/home/user/file_cor.bam
Rscript run_spp_nodups.R -c=$INPUT -savp -out=$OUTPUT
出现错误:
Error: could not find function "runmean"
Execution halted
函数“ runmean”属于我安装的软件包“ caTools”。 R版本是合适的,因为我机器上的R版本是3.3.2,而'caTools'取决于R(≥2.2.0)。
“ run_spp_nodups.R”的R代码很大,要粘贴在这里。我只显示了带有runmean的部分:
# Smooth the cross-correlation curve if required
cc <- crosscorr$cross.correlation
crosscorr$min.cc <- crosscorr$cross.correlation[ length(crosscorr$cross.correlation$y) , ] # minimum value and shift of cross-correlation
cat("Minimum cross-correlation value", crosscorr$min.cc$y,"\n",file=stdout())
cat("Minimum cross-correlation shift", crosscorr$min.cc$x,"\n",file=stdout())
sbw <- 2*floor(ceiling(5/iparams$sep.range[2]) / 2) + 1 # smoothing bandwidth
cc$y <- runmean(cc$y,sbw,alg="fast")
正在发生什么以及如何解决?