I'm getting an error when running runSeqsumm to read a file. Here's my code and the error message:
> library(seqCNA) > runSeqsumm(file='C1.sam') Reading records from C1.sam Summarizing in windows of 50000 bps. sh: c®Ä: command not found Writing results to disk... Done! >
Naturally I've checked that the file in question (C1.sam) exists and is a valid sam file. The program creates a file with a garbage name (eg, '%E7%D4I???') containing just a header line:
chrom win.start reads.gc reads.mapq counts
The guts of runSeqsumm is a C++ function. It kind of feels like the R-C++ linkage is not right. I'm not an R package developer and don't know how to "remake" the package (other than by reinstalling it from scratch), but if someone can point me in that direction I'm happy to give it a try.
Here's the sessionInfo
> sessionInfo() R version 3.2.3 (2015-12-10) Platform: x86_64-apple-darwin13.4.0 (64-bit) Running under: OS X 10.10.5 (Yosemite) locale: [1] C attached base packages: [1] stats graphics grDevices utils datasets methods base other attached packages: [1] seqCNA_1.16.0 seqCNA.annot_1.6.0 adehabitatLT_0.3.20 [4] CircStats_0.2-4 boot_1.3-17 MASS_7.3-45 [7] adehabitatMA_0.3.10 ade4_1.7-3 sp_1.2-2 [10] doSNOW_1.0.14 snow_0.4-1 iterators_1.0.8 [13] foreach_1.4.3 GLAD_2.34.0 loaded via a namespace (and not attached): [1] codetools_0.2-14 grid_3.2.3 lattice_0.20-33
For some reason the author of the package appears to have decided to require that the user have samtools on their PATH. Hopefully they will update their package to use Rsamtools or bamsignals instead.