tests/ash_error.R

### ash error

ll <- get(load("ash_error_list.rda"))

m <- ll$m
betahat <- ll$betahat
sebetahat <- ll$sebetahat

ashr::ash(betahat, sebetahat)

library(ashr)
comp_postmean(m, betahat, sebetahat)
kkdey/flashr documentation built on May 20, 2019, 10:36 a.m.