coloc.susie returns NA results without error message
zhilizheng opened this issue · 3 comments
zhilizheng commented
Hi Dr. Wallace,
Thanks for the great coloc tool. I meet a problem when running the coloc.susie, that I don't understand where the problem comes from.
Relicate of the problem:
Data: https://www.dropbox.com/scl/fo/f0akgky4bbnglgiejykf9/h?rlkey=4zu9ki42hdu3c7323s6l2hwjd&dl=0
# finemap from susie for summary data 1
f1 = readRDS("f1.rds")
# finemap from susie for summary data 2
f2 = readRDS("f2.rds")
# coloc
res=coloc.susie(f1, f2)
As shown in the example, the coloc results is NA, without any errors. Could you figure me what happened here?
BTW, the suma.rds, is the orginal summary data. Columns with name BETA, SE, P, MAF, N were from summary data1, BETA2, SE2, P2, MAF2, N2 were from summary data2.
Thanks for your help.
Regards,
Zhili
htqqdd commented
I have the same problem. Have you solve it?
chr1swallace commented
Sorry, I can't reproduce this. Could either of you send me a small example dataset that produces this error?
…________________________________
From: htqqdd ***@***.***>
Sent: Tuesday, April 9, 2024 2:29 PM
To: chr1swallace/coloc ***@***.***>
Cc: Subscribed ***@***.***>
Subject: Re: [chr1swallace/coloc] coloc.susie returns NA results without error message (Issue #147)
I have the same problem. Have you solve it?
—
Reply to this email directly, view it on GitHub<#147 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AAQWR2CHBH5UFNLHRXHJN2TY4PUMFAVCNFSM6AAAAABDZBOR3CVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDANBVGE4DGOBZGU>.
You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>
chr1swallace commented
I just saw the data - thanks. Easy to check. The issue is that susie doesn't find any credible sets in your second dataset (shown below). That means there is nothing to colocalise. I will make a more informative warning message for this situation.
f2$sets
$cs
NULL
$coverage
NULL
$requested_coverage
[1] 0.95
…________________________________
From: Chris Wallace ***@***.***>
Sent: Wednesday, April 10, 2024 3:34 PM
To: chr1swallace/coloc ***@***.***>; chr1swallace/coloc ***@***.***>
Cc: Subscribed ***@***.***>
Subject: Re: [chr1swallace/coloc] coloc.susie returns NA results without error message (Issue #147)
Sorry, I can't reproduce this. Could either of you send me a small example dataset that produces this error?
________________________________
From: htqqdd ***@***.***>
Sent: Tuesday, April 9, 2024 2:29 PM
To: chr1swallace/coloc ***@***.***>
Cc: Subscribed ***@***.***>
Subject: Re: [chr1swallace/coloc] coloc.susie returns NA results without error message (Issue #147)
I have the same problem. Have you solve it?
—
Reply to this email directly, view it on GitHub<#147 (comment)>, or unsubscribe<https://github.com/notifications/unsubscribe-auth/AAQWR2CHBH5UFNLHRXHJN2TY4PUMFAVCNFSM6AAAAABDZBOR3CVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMZDANBVGE4DGOBZGU>.
You are receiving this because you are subscribed to this thread.Message ID: ***@***.***>