get_map: Try to find the MAP using coordinate descent

Description Usage Arguments

Description

v4_li_mat_func and v4_prior_func use working model 4 described here: https://jean997.github.io/sh2Ash/mh_test.html

Usage

1
2
3
get_map(pars_start, data, grid, par_bounds, tol = 1e-07, n.iter = Inf,
  type = c("null", "causal", "shared"), fix_pi = FALSE,
  fix_pars = c(), null_wt = 10, ...)

Arguments

pars_start

Starting parameters: eg rho, b, q. Should be appropriate for v4_li_mat_func

data

data

grid

Should be a data.frame with at least two columns, S1 and S2.

par_bounds

should be data.frame or matrix number of pars by 2

type

One of the models null, causal or shared. If null, pars should be length 1 (rho). If causal, pars should be length 2 (rho, b), If shared pars should be length 3 (rho, b, q) The prior on pi is hard coded as dirichlet((10, 1, 1, ...))

fix_pi

if FALSE, only estimate pars, leave mixing parameters fixed as given in grid

fix_pars

vector of indices of parameters to fix

null_wt

Specifies the prior weight on the first entry of grid

...

Additional parameters (z_prior_func, b_prior_func, q_prior_func)


jean997/sherlockAsh documentation built on May 18, 2019, 11:45 p.m.