error in hclust (--pcatype select) #47
Description
Hello,
The first step Rscript dchicf.r --file input_files.txt --pcatype cis is working very well for the compartment calling.
But I am getting an error for the second step Rscript dchicf.r --file input_files.txt --pcatype select:
Error in hclust(as.dist(round(1 - cor(pc.mat), 4))) :
NA/NaN/Inf in foreign function call (arg 10)
Calls: pcselect -> pcselectioncore -> hclust
Execution halted
It seems that it is working for the first two samples as the stdout return:
Running intra 1 in poll_0197 sample
Running intra 1 in poll_3654 sample
Here is my input files:
Bovin-0197.ARS-UCD1.2.mapq_10.50000.txt Bovin-0197.ARS-UCD1.2.mapq_10.50000.bed poll_0197 poll
Bovin-3654.ARS-UCD1.2.mapq_10.50000.txt Bovin-3654.ARS-UCD1.2.mapq_10.50000.bed poll_3654 poll
Bovin-669.ARS-UCD1.2.mapq_10.50000.txt Bovin-669.ARS-UCD1.2.mapq_10.50000.bed unp_669 unp
Bovin-977.ARS-UCD1.2.mapq_10.50000.txt Bovin-977.ARS-UCD1.2.mapq_10.50000.bed unp_977 unp
Any idea?