search

jrs95 / hyprcoloc

Hypothesis Prioritisation in multi-trait Colocalization
https://jrs95.github.io/hyprcoloc/
GNU General Public License v3.0
46 stars 12 forks source link

Error in which.max(labf) : 'list' object cannot be coerced to type 'double' #27

Closed WaheedURAhmed1 closed 1 year ago

WaheedURAhmed1 commented 1 year ago

Hi there, I am trying to run the hyprcoloc package, but getting the below error.

hyprcoloc(loc01_beta_matrix, loc01_se_matrix, trait.names=traits, snp.id=rsid_loc01) Error in which.max(labf) : 'list' object cannot be coerced to type 'double'

here's a header of all files to help you:

head(loc01_beta_matrix) T1 T2 T3 T4 T5 T6 T7 T8 T9 T10 T11 T12 rs12139158 0.0253761 -0.016001100 -0.01320890 -0.000542281 0.0627357 0.01115330 -2.61864e-02 0.07725800 -0.006135420 -0.0145061 0.0362987 -0.0381493 rs12064209 -0.0302395 -0.026537200 -0.04993110 0.049364300 0.0093351 0.06371270 5.03399e-02 -0.02919430 -0.100037000 0.1210630 0.0400980 -0.2859180 rs17034627 0.0295129 0.005371080 -0.01997840 -0.021614600 -0.0829874 -0.00970612 4.87015e-02 0.00238708 0.057445100 -0.0258380 0.0446427 -0.1901770 rs9332414 -0.0203766 -0.000274617 -0.00301425 0.007538970 -0.0287159 0.00426423 6.48232e-05 -0.05930950 -0.000932442 0.0223630 0.0265988 -0.0354991 rs115221625 -0.0175293 0.008864500 0.02718990 0.043303200 -0.1099040 0.00459884 -5.09681e-02 0.01865400 -0.106147000 0.1523170 0.1734120 -0.1494930 rs12141246 -0.0202726 0.000263744 -0.00195736 0.007033050 -0.0291631 0.00463914 1.21057e-03 -0.05901810 -0.001497740 0.0272348 0.0247724 -0.0373881 head(loc01_se_matrix) T1 T2 T3 T4 T5 T6 T7 T8 T9 T10 T11 T12 rs12139158 0.0151063 0.0155211 0.0155490 0.0205578 0.0201483 0.0239369 0.0322175 0.0452013 0.0592372 0.0545188 0.0618186 0.0677466 rs12064209 0.0316663 0.0323276 0.0325916 0.0427048 0.0402632 0.0501843 0.0666161 0.0798421 0.1249230 0.1123670 0.1333200 0.1707720 rs17034627 0.0376127 0.0380625 0.0383871 0.0505247 0.0479451 0.0581393 0.0788080 0.0963338 0.1566580 0.1420700 0.1489930 0.1698140 rs9332414 0.0105585 0.0107960 0.0108623 0.0144101 0.0151486 0.0169013 0.0222721 0.0308595 0.0417838 0.0377977 0.0443587 0.0471825 rs115221625 0.0302246 0.0309781 0.0308677 0.0416579 0.0428905 0.0482738 0.0652175 0.0776992 0.1247510 0.1073410 0.1342290 0.1352530 rs12141246 0.0105535 0.0107901 0.0108567 0.0144007 0.0150632 0.0168901 0.0222639 0.0308681 0.0417423 0.0378083 0.0443200 0.0471482

head(traits) [1] "T1" "T2" "T3" "T4" "T5" "T6"
head(rsid_loc01) [1] "rs12139158" "rs12064209" "rs17034627" "rs9332414" "rs115221625" "rs12141246"

Please may you kindly advise?

Thanks!

WaheedURAhmed1 commented 1 year ago

@cnfoley would appreciate your help! thanks!