Java error while running HiCDC overview code
1
0
Entering edit mode
IsMM • 0
@76f106c3
Last seen 11 weeks ago
México

Hey folks!

I am trying to run the code provided in the HiCDC overview (https://github.com/mervesa/HiCDCPlus#diff_int).

Everything runs as expected (as shown in the Github code) but at the moment of trying to save the output of HiCDCPlus_parallel() into a .hic file with hicdc2hic() I run into an error I am not sure how to interpret/solve. Also, the .hic file is saved but it cannot be opened properly in Juicebox, the screen shows only white.

I was wondering if someone had run into something similar. Thanks :)

The error reads:

Writing body
java.lang.NumberFormatException: For input string: "score"         
at java.base/jdk.internal.math.FloatingDecimal.readJavaFormatString(FloatingDecimal.java:2054)
        at java.base/jdk.internal.math.FloatingDecimal.parseFloat(FloatingDecimal.java:122)
        at java.base/java.lang.Float.parseFloat(Float.java:455)
        at juicebox.tools.utils.original.mnditerator.MNDLineParser.updatePairScoreIfNeeded(MNDLineParser.java:88)
        at juicebox.tools.utils.original.mnditerator.MNDFileParser.parseShortFormat(MNDFileParser.java:105)
        at juicebox.tools.utils.original.mnditerator.MNDFileParser.parse(MNDFileParser.java:87)
        at juicebox.tools.utils.original.mnditerator.GenericPairIterator.advance(GenericPairIterator.java:56)
        at juicebox.tools.utils.original.mnditerator.AsciiPairIterator.<init>(AsciiPairIterator.java:55)
        at juicebox.tools.utils.original.mnditerator.PairIterator.getIterator(PairIterator.java:54)
        at juicebox.tools.utils.original.Preprocessor.computeWholeGenomeMatrix(Preprocessor.java:583)
        at juicebox.tools.utils.original.Preprocessor.writeBody(Preprocessor.java:674)
        at juicebox.tools.utils.original.Preprocessor.preprocess(Preprocessor.java:436)
        at juicebox.tools.clt.old.PreProcessing.run(PreProcessing.java:165)
        at juicebox.tools.HiCTools.main(HiCTools.java:94)

and the code is:

library(HiCDCPlus) #I am using R 4.3.2 

hicfile_path<-system.file("extdata", "GSE63525_HMEC_combined_example.hic", package = "HiCDCPlus")
outdir<-tempdir(check=TRUE)

#generate features
construct_features(output_path=paste0(outdir,"/hg19_50kb_GATC"),
                   gen="Hsapiens",gen_ver="hg19",
                   sig="GATC",
                   bin_type="Bins-uniform",
                   binsize=50000,
                   chrs=c("chr21","chr22"))

gi_list<-generate_bintolen_gi_list(bintolen_path=paste0(outdir,"/hg19_50kb_GATC_bintolen.txt.gz"))
gi_list<-add_hic_counts(gi_list,hic_path = hicfile_path)
gi_list<-expand_1D_features(gi_list)

#run HiC-DC+ on 2 cores
set.seed(1010) #HiC-DC downsamples rows for modeling
gi_list<-HiCDCPlus_parallel(gi_list,ncore=2)
head(gi_list)
#~~~~~~Up to here everything works fine! The next line throws the error 
hicdc2hic(gi_list,hicfile=paste0(outdir,'/GSE63525_HMEC_combined_result.hic'), mode='normcounts',gen_ver='hg19')
JuicerTools HiC HiCDCPlus Java Juicer • 207 views
ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 1 hour ago
United States

This is a Bioconductor support site. You might try over at biostars.org, or better yet, ask on the GitHub where you got that package.

Login before adding your answer.

Traffic: 836 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6