compartment_score - Error in if (any(newwidths > 2 * widths)) { : missing value where TRUE/FALSE needed #373
Description
Hi, all
I am encounting similar error as thread before when calculating compartment_score using a list of hic inputs.
It seems that it might be due to the difference between male and female samples within my project.
data$cool_objects$GW25T103C07$CHRS
[1] "1" "10" "11" "12" "13" "14" "15" "16" "17" "18" "19" "2" "20" "21" "22"
[16] "3" "4" "5" "6" "7" "8" "9" "M" "X"
data$cool_objects$GW25T103C02$CHRS
[1] "1" "10" "11" "12" "13" "14" "15" "16" "17" "18" "19" "2" "20" "21" "22"
[16] "3" "4" "5" "6" "7" "8" "9" "M" "X" "Y"
$CENTROMERES
chrom start end
1 1 243 249
2 10 3431 3437
3 11 3724 3730
4 12 3961 3967
5 13 4159 4161
6 14 4352 4353
7 15 4529 4535
8 16 4769 4775
9 17 4924 4930
10 18 5077 5083
11 19 5257 5263
12 2 682 688
13 20 5378 5384
14 21 5467 5473
15 22 5544 5550
16 3 1164 1170
17 4 1479 1485
18 5 1853 1859
19 6 2241 2247
20 7 2582 2588
21 8 2872 2878
22 9 3170 3176
23 X 5732 5738
24 Y NA NA
Any suggestions will be of great help.
Best,
Junfeng