aidenlab / juicer

A One-Click System for Analyzing Loop-Resolution Hi-C Experiments
http://aidenlab.org
MIT License
410 stars 181 forks source link

Unknown command: LibraryComplexity #236

Closed edfajardo closed 3 years ago

edfajardo commented 3 years ago

Line 453 of the juicer.sh script (version 1.6 in CPU directory) results in the following error:

Exception in thread "main" java.lang.RuntimeException: Unknown command: LibraryComplexity
    at juicebox.tools.HiCTools.main(HiCTools.java:99)

The offending line in juicer.sh is:

${juiceDir}/scripts/common/juicer_tools LibraryComplexity $outputdir inter.txt >> $outputdir/inter.txt

This error is not dependent on the version of the version_tools.jar installed (from https://github.com/aidenlab/juicer/wiki/Download).

The output of just running juicer_tools (in the common directory) with juicer_tools.1.6.2_jcuda.0.7.0.jar installed is:


$ juicer_tools
Juicer Tools Version 1.6.2
Usage:
    dump <observed/oe> <NONE/VC/VC_SQRT/KR> <hicFile(s)> <chr1>[:x1:x2] <chr2>[:y1:y2] <BP/FRAG> <binsize> [outfile]
    dump <norm/expected> <NONE/VC/VC_SQRT/KR> <hicFile(s)> <chr> <BP/FRAG> <binsize> [outfile]
    dump <loops/domains> <hicFile URL> [outfile]
    pre [options] <infile> <outfile> <genomeID>
    addNorm <input_HiC_file>
    pearsons [-p] <NONE/VC/VC_SQRT/KR> <hicFile(s)> <chr> <BP/FRAG> <binsize> [outfile]
    eigenvector -p <NONE/VC/VC_SQRT/KR> <hicFile(s)> <chr> <BP/FRAG> <binsize> [outfile]
    apa <hicFile(s)> <PeaksFile> <SaveFolder>
    arrowhead <hicFile(s)> <output_file>
    hiccups <hicFile> <outputDirectory>
    hiccupsdiff <firstHicFile> <secondHicFile> <firstLoopList> <secondLoopList> <outputDirectory>
    validate <hicFile>
    -h, --help print help
    -v, --verbose verbose mode
    -V, --version print version
Type juicer_tools <commandName> for more detailed usage instructions

So clearly, there is no LibraryComplexity command available. I also tried juicer_tools_1.19.02.jar and juicer_tools_0.7.0.jar; neither has such command.

Following immediately the above error line, I found the following (unrelated?) error line:

sort: cannot create temporary file in ‘/home/project/fastq/sample1/HIC_tmp’: No such file or directory

Which is odd because HIC_tmp had been created earlier (I saw the directory soon after the juicer.sh run started) but obviously was deleted at a later step in the script.

In the end, the .hic files were created OK, and I assume these are correct.

Thank you.

nchernia commented 3 years ago

The more recent jars have LibraryComplexity even though it's not in the usage statement

On Fri, Aug 20, 2021 at 4:13 PM edfajardo @.***> wrote:

Line 453 of the juicer.sh script (version 1.6 in CPU directory) results in the following error:

Exception in thread "main" java.lang.RuntimeException: Unknown command: LibraryComplexity

at juicebox.tools.HiCTools.main(HiCTools.java:99)

The offending line in juicer.sh is:

${juiceDir}/scripts/common/juicer_tools LibraryComplexity $outputdir inter.txt >> $outputdir/inter.txt

This error is not dependent on the version of the version_tools.jar installed (from https://github.com/aidenlab/juicer/wiki/Download).

The output of just running juicer_tools (in the common directory) with juicer_tools.1.6.2_jcuda.0.7.0.jar installed is:

$ juicer_tools

Juicer Tools Version 1.6.2

Usage:

dump <observed/oe> <NONE/VC/VC_SQRT/KR> <hicFile(s)> [:x1:x2] [:y1:y2] <BP/FRAG> [outfile]

dump <norm/expected> <NONE/VC/VC_SQRT/KR> <hicFile(s)> <BP/FRAG> [outfile]

dump <loops/domains> [outfile]

pre [options]

addNorm

pearsons [-p] <NONE/VC/VC_SQRT/KR> <hicFile(s)> <BP/FRAG> [outfile]

eigenvector -p <NONE/VC/VC_SQRT/KR> <hicFile(s)> <BP/FRAG> [outfile]

apa <hicFile(s)>

arrowhead <hicFile(s)>

hiccups

hiccupsdiff

validate

-h, --help print help

-v, --verbose verbose mode

-V, --version print version

Type juicer_tools for more detailed usage instructions


So clearly, there is no LibraryComplexity command available. I also tried juicer_tools_1.19.02.jar and juicer_tools_0.7.0.jar; neither has such command.

Following immediately the above error line, I found the following (unrelated?) error line:

sort: cannot create temporary file in ‘/home/project/fastq/sample1/HIC_tmp’: No such file or directory

Which is odd because HIC_tmp had been created earlier (I saw the directory soon after the juicer.sh run started) but obviously was deleted at a later step in the script.

In the end, the .hic files were created OK, and I assume these are correct.

Thank you.

— You are receiving this because you are subscribed to this thread. Reply to this email directly, view it on GitHub https://github.com/aidenlab/juicer/issues/236, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAK2EW6QNAGU66KU7C2UKN3T52ZM7ANCNFSM5CQ6T34Q . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&utm_campaign=notification-email .

-- Neva Cherniavsky Durand, Ph.D. | she, her, hers Assistant Professor | Molecular and Human Genetics Aiden Lab | Baylor College of Medicine www.aidenlab.org

edfajardo commented 3 years ago

Thank you for your reply. As mentioned in my message, I tried many versions of the juicer tools jar from https://github.com/aidenlab/juicer/wiki/Download, the latest being 1.19.02. The point I am trying to make is, running the juicer.sh script throws an error. The error happens at line 453 of the CPU version because it calls LibraryComplexity and it thinks LibraryComplexity is not there. Whether it is true of not that the jar has LibraryComplexity, the output message is that LibraryComplexity is not a valid command (as documented in my original post above).

Thank you.

nchernia commented 3 years ago

We are able to run this command without a problem. I would suggest trying to run the command, not just looking at usage. I cannot reproduce your problem.

On Mon, Aug 30, 2021 at 2:25 PM edfajardo @.***> wrote:

Thank you for your reply. As mentioned in my message, I tried many versions of the juicer tools jar from https://github.com/aidenlab/juicer/wiki/Download, the latest being 1.19.02. The point I am trying to make is, running the juicer.sh script throws an error. The error happens at line 453 of the CPU version because it calls LibraryComplexity and it thinks LibraryComplexity is not there. Whether it is true of not that the jar has LibraryComplexity, the output message is that LibraryComplexity is not a valid command (as documented in my original post above).

Thank you.

— You are receiving this because you commented.

Reply to this email directly, view it on GitHub https://github.com/aidenlab/juicer/issues/236#issuecomment-908580263, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAK2EW45UHTRMQSIU7OAWRTT7PEIPANCNFSM5CQ6T34Q . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

-- Neva Cherniavsky Durand, Ph.D. | she, her, hers Assistant Professor | Molecular and Human Genetics Aiden Lab | Baylor College of Medicine www.aidenlab.org

edfajardo commented 3 years ago

OK, thank you. It is working now (tried both version 1.19.02 and 1.22.01 of the jar). I assume that the .hic files generated in runs where the LibraryComplexity command failed are OK. Is this correct?

Thank you very much.

nchernia commented 3 years ago

Yes they are fine.

On Mon, Aug 30, 2021 at 3:47 PM edfajardo @.***> wrote:

OK, thank you. It is working now (tried both version 1.19.02 and 1.22.01 of the jar). I assume that the .hic files generated in runs where the LibraryComplexity command failed are OK. Is this correct?

Thank you very much.

— You are receiving this because you commented.

Reply to this email directly, view it on GitHub https://github.com/aidenlab/juicer/issues/236#issuecomment-908639147, or unsubscribe https://github.com/notifications/unsubscribe-auth/AAK2EWZI34WGMMM472YGLC3T7PN6LANCNFSM5CQ6T34Q . Triage notifications on the go with GitHub Mobile for iOS https://apps.apple.com/app/apple-store/id1477376905?ct=notification-email&mt=8&pt=524675 or Android https://play.google.com/store/apps/details?id=com.github.android&referrer=utm_campaign%3Dnotification-email%26utm_medium%3Demail%26utm_source%3Dgithub.

-- Neva Cherniavsky Durand, Ph.D. | she, her, hers Assistant Professor | Molecular and Human Genetics Aiden Lab | Baylor College of Medicine www.aidenlab.org

edfajardo commented 3 years ago

Thanks!