Skip to content

cmdError in system2(path, "callpeak -h", stdout = TRUE, stderr = TRUE) #2314

@aowenson-imm

Description

@aowenson-imm

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.

Metadata

Metadata

Assignees

No one assigned

    Labels

    bugSomething isn't working

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions