Ester

ash does not return logLR i.e. if res=ash(...), there is no res$logLR I have data where this happens: `load("/mnt/lustre/home/epantaleo/src/multiseq/local/testing.pipeline.RData")` `library(multiseq) #pull the most recent version` `zdat.ash = ash(zdat[3,ind], zdat[4,ind], prior=prior,...