Attach your log file
ArchR.log
Describe the bug
ArchR function addReproduciblePeakSet tries to call program callpeak here:
|
cmd <- sprintf("callpeak -g %s --name %s --treatment %s --outdir %s --format BED --call-summits --keep-dup all %s", |
Problem is that is not a program. It is a command of program macs2, so you actually want macs2 callpeak.
$ callpeak
callpeak: command not found
$ macs2 callpeak
usage: macs2 callpeak [-h] -t TFILE [TFILE ...] ...
Similarly for other macs2 commands:
usage: macs2 [-h] [--version] {callpeak,bdgpeakcall,bdgbroadcall,bdgcmp,bdgopt,cmbreps,bdgdiff,filterdup,predictd,pileup,randsample,refinepeak} ...
I'm confused current ArchR can work for other users. Did an old version of Macs2 create aliases but not anymore? Are they manually creating aliases?
To Reproduce
python -m venv macs2
source macs2/bin/activate
pip install macs2
callpeak
# callpeak: command not found
macs2 callpeak
# usage: macs2 callpeak [-h] ...
Additional information
I'm aware of closed issue #1483, but replies are being ignored. Those replies corroborate my experience.
Attach your log file
ArchR.log
Describe the bug
ArchR function
addReproduciblePeakSettries to call programcallpeakhere:ArchR/R/ReproduciblePeakSet.R
Line 788 in 6feec35
Problem is that is not a program. It is a command of program
macs2, so you actually wantmacs2 callpeak.Similarly for other macs2 commands:
I'm confused current ArchR can work for other users. Did an old version of Macs2 create aliases but not anymore? Are they manually creating aliases?
To Reproduce
Additional information
I'm aware of closed issue #1483, but replies are being ignored. Those replies corroborate my experience.