search

stephens999 / ashr

An R package for adaptive shrinkage
GNU General Public License v3.0
80 stars 36 forks source link

Negative PositiveProb Values Returned #23

Closed jhmarcus closed 8 years ago

jhmarcus commented 8 years ago

Data

pgc_scz.sumstat.gz

Code

library(ashr)
df <- read.table("pgc_scz.sumstat.gz", header=TRUE, stringsAsFactors=FALSE)
ash_res <- ash(df$BETA, df$SE)

# negative value returned
min(ash_res$PositiveProb)
stephens999 commented 8 years ago

how do you put the data in the issue like that? ... oh i see it just drags and drops.. very cool.