calc_gloglik: Compute loglikelihood for data usign the prior g (usually from an ash fit)
Description
Return the log-likelihood of the data betahat, with standard errors betahatsd,
under the fitted distribution in the ash object.
Usage
calc_gloglik(g, betahat, betahatsd, df, model = c("EE", "ES"))
Arguments
g
the prior for effects or standardized effects
betahatsd
the observed standard errors
df
appropriate degrees of freedom for (t) distribution of betahat/sebetahat
model
indicates whether you want the likelihood under the EE or ES model
Details
See example in CompareBetahatvsZscoreAnalysis.rmd