Open gagelo01 opened 3 years ago
Hi, in the function gwasvcf_to_coloc, the code chunk : else if (type2 == "cc") { out2$s <- mean(tab1$NC/tab1$SS, na.rm = TRUE) }
should be replaced by if (type2 == "cc") { out2$s <- mean(tab2$NC/tab2$SS, na.rm = TRUE) }
Hi, in the function gwasvcf_to_coloc, the code chunk : else if (type2 == "cc") { out2$s <- mean(tab1$NC/tab1$SS, na.rm = TRUE) }
should be replaced by if (type2 == "cc") { out2$s <- mean(tab2$NC/tab2$SS, na.rm = TRUE) }
1 the supression of the else statement allow for colocalisation between two case-control traits.
2 the replacement of tab1 by tab2 allow to get information from the good object.