chr1swallace / coloc

Repo for the R package coloc
139 stars 44 forks source link

How to find H0 to H3 from the coloc.abf output for each snp #85

Closed zillurbmb51 closed 2 years ago

zillurbmb51 commented 2 years ago

Hi, The coloc.abf gives us summary (PP for H0-4) and dataframe with the following columns: "snp" "position" "V.df1" "z.df1" "r.df1" "lABF.df1" "V.df2" "z.df2" "r.df2" "lABF.df2" "internal.sum.lABF" "SNP.PP.H4"

How could we find the PP for other hypothesis (H0-3) for each snp?

chr1swallace commented 2 years ago

H1, h2, use finemap.abf. hp doesn't correspond to any individual snp. If you believe h3 is true, then I would use the per trait per snp results from finemap.abf to infer the most likely causal variants for each.

https://chr1swallace.github.io


From: zillurbmb51 @.> Sent: Sunday, May 8, 2022 3:29:59 PM To: chr1swallace/coloc @.> Cc: Subscribed @.***> Subject: [chr1swallace/coloc] How to find H0 to H3 from the coloc.abf output for each snp (Issue #85)

Hi, The coloc.abf gives us summary (PP for H0-4) and dataframe with the following columns: "snp" "position" "V.df1" "z.df1" "r.df1" "lABF.df1" "V.df2" "z.df2" "r.df2" "lABF.df2" "internal.sum.lABF" "SNP.PP.H4"

How could we find the PP for other hypothesis (H0-3) for each snp?

— Reply to this email directly, view it on GitHubhttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fchr1swallace%2Fcoloc%2Fissues%2F85&data=05%7C01%7Ccew54%40universityofcambridgecloud.onmicrosoft.com%7C6921146a08e44da046a108da30ff3c44%7C49a50445bdfa4b79ade3547b4f3986e9%7C0%7C0%7C637876170026337385%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=45fyF3VUzBRO30ElbcmmbhZp6ArDdZnQUguYHzXY8bc%3D&reserved=0, or unsubscribehttps://eur03.safelinks.protection.outlook.com/?url=https%3A%2F%2Fgithub.com%2Fnotifications%2Funsubscribe-auth%2FAAQWR2FSOHTDROXCZ23RIGTVI7FWPANCNFSM5VMAKHGQ&data=05%7C01%7Ccew54%40universityofcambridgecloud.onmicrosoft.com%7C6921146a08e44da046a108da30ff3c44%7C49a50445bdfa4b79ade3547b4f3986e9%7C0%7C0%7C637876170026337385%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=uOzaCostLWQO5ueFcQVutRAayCq10EJTSunPDbQBjd4%3D&reserved=0. You are receiving this because you are subscribed to this thread.Message ID: @.***>