search

chr1swallace / coloc

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

coloc.abf error #141

Closed maguileraf closed 7 months ago

maguileraf commented 7 months ago

Hi Chris, I am trying to use coloc for the first time but I keep getting this error: Error in coloc.abf(dataset1 = D1, dataset2 = D2) : dataset1 and dataset2 should contain the same snps in the same order, or should contain snp names through which the common snps can be identified

However, the snps in both lists are the same (see attached). Any ideas on what could be the issue? Screenshot 2023-11-21 at 7 56 04 AM Screenshot 2023-11-21 at 7 55 53 AM

jjuhyunkim commented 7 months ago

Did you figure out how to resolve this error? I encountered the same issue. 스크린샷 2023-12-06 오전 1 11 56

chr1swallace commented 7 months ago

I don't see this error. Could you share a minimal dataset that exhibits this issue?

-- https://chr1swallace.github.io

From: Juhyun Kim @.> Sent: Wednesday, December 6, 2023 6:12:14 AM To: chr1swallace/coloc @.> Cc: Subscribed @.***> Subject: Re: [chr1swallace/coloc] coloc.abf error (Issue #141)

Did you figure out how to resolve this error? I encountered the same issue. 2023-12-06.1.11.56.png (view on web)https://github.com/chr1swallace/coloc/assets/131193680/05fdefa7-68d1-4578-aab9-766b5f382f42

— Reply to this email directly, view it on GitHubhttps://github.com/chr1swallace/coloc/issues/141#issuecomment-1842158475, or unsubscribehttps://github.com/notifications/unsubscribe-auth/AAQWR2H3YC3L5RLXH6L6KFLYIAED5AVCNFSM6AAAAAA7UQ2JKGVHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTQNBSGE2TQNBXGU. You are receiving this because you are subscribed to this thread.Message ID: @.***>

jjuhyunkim commented 7 months ago

Thank you for your reply! ☺️

I have identified the issue during the data export process.

The problem arises from using hg19 for dataset1 and hg38 for dataset2. I mistakenly did not replace the hg38 positions after liftovering variants from hg19.

LLH0514 commented 8 hours ago

I encountered the same issue.Error in coloc.abf(dataset1 = list(pvalues = input$pval_gwas, N = as.integer(input$samplesize_gwas), : dataset1 and dataset2 should contain the same snps in the same order, or should contain snp names through which the common snps can be identified.Did you figure out how to resolve this error?l use hg19 for dataset1 and hg19 for dataset2