Skip to content

Commit

Permalink
updated sign test
Browse files Browse the repository at this point in the history
  • Loading branch information
explodecomputer committed Aug 4, 2018
1 parent c787765 commit bec1ea7
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions tests/test_mr_sign.R
Expand Up @@ -21,12 +21,12 @@ library(dplyr)
group_by(param, n) %>% summarise(minp = min(pval))

# Looks like 6 is the smallest number of SNPs that can achieve p < 0.05

library(TwoSampleMR)
a <- extract_instruments(2)
b <- extract_outcome_data(a$SNP,7)
dat <- harmonise_data(a, b)
mr(dat, method_list=c("mr_ivw", "mr_sign"))

mr_scatter_plot(mr(dat, method_list=c("mr_ivw", "mr_sign")), dat)

# Could consider doing a parametric bootstrap type analysis

0 comments on commit bec1ea7

Please sign in to comment.