mash_compute_vloglik: Compute vector of loglikelihood for fitted mash object on new...

Description Usage Arguments Details Value Examples

View source: R/likelihoods_origdata.R

Description

Compute vector of loglikelihood for fitted mash object on new data

Usage

1
mash_compute_vloglik(g, data, algorithm.version = c("Rcpp", "R"))

Arguments

g

A mash object.

data

A set of data on which to compute the loglikelihood.

algorithm.version

Indicate R or Rcpp version

Details

The log-likelihood for each element is p(Bhat_j | Shat_j,g,α) where Bhat_j | B_j, Shat_j \sim N(B_j, Shat_j) and B_j/Shat_j^α | Shat_j \sim g Here the value of α is set when setting up the data object in 'mash_set_data'. If g is a mash object (safest!) then the function will check that this value matches the α used when fitting 'mash'. Note: as a convenience, this function can also be called with g a mixture distribution with same structure as the fitted_g from a mash object. This is mostly useful when doing simulations, where you might want to compute the likelihood under the "true" g. When used in this way the user is responsible for making sure that the g makes sense with the alpha set in data.

Value

The vector of log-likelihoods for each data point computed using g.

Examples

1
2
3
4
simdata = simple_sims(50,5,1)
data = mash_set_data(simdata$Bhat, simdata$Shat)
m = mash(data, cov_canonical(data))
mash_compute_vloglik(m,data)

mashr documentation built on May 24, 2021, 1:06 a.m.