Open matteofloris opened 3 years ago
MAF information is missing in several datasets, and colocalisation cannot be performed. An example:
GWAS <- "ebi-a-GCST005063" chrpos <- "11:61396322-61796322" out <- ieugwasr_to_coloc(id1=GWAS, id2=GWAS, chrompos=chrpos) res <- coloc::coloc.abf(out[[1]], out[[2]]) Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 0 (non-NA) cases Inoltre: Warning message: In sdY.est(d$varbeta, d$MAF, d$N) : estimating sdY from maf and varbeta, please directly supply sdY if known
MAF information is missing in several datasets, and colocalisation cannot be performed. An example:
GWAS <- "ebi-a-GCST005063" chrpos <- "11:61396322-61796322" out <- ieugwasr_to_coloc(id1=GWAS, id2=GWAS, chrompos=chrpos) res <- coloc::coloc.abf(out[[1]], out[[2]]) Error in lm.fit(x, y, offset = offset, singular.ok = singular.ok, ...) : 0 (non-NA) cases Inoltre: Warning message: In sdY.est(d$varbeta, d$MAF, d$N) : estimating sdY from maf and varbeta, please directly supply sdY if known