View source: R/map_pi_rho_nonnull.R
Currently this function is not exported and is mostly for development purposes. It is just like map_pi_rho except that one can specify gamma, eta, and q rather than having them fixed at zero.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 | map_pi_rho_nonnull(
X,
mix_grid,
sigma_g,
qalpha,
qbeta,
q_start = qbeta(0.5, qalpha, qbeta),
gamma_start = 0,
eta_start = 0,
rho_start = 0,
tol = 1e-05,
n.iter = 20,
null_wt = 10,
z_prior_func = function(z) { dnorm(z, 0, 0.5, log = TRUE) },
q_prior_func = function(q) { dbeta(q, qalpha, qbeta, log = TRUE) },
gammaeta_prior_func = function(g) { dnorm(g, 0, sigma_g, log = TRUE) },
optmethod = c("mixSQP", "mixIP")
)
|
X |
An object of class cause_data containing data for the two traits. |
mix_grid |
An object of class cause_grid containing variance pair candidates |
rho_start |
Starting value for rho |
null_wt |
Specifies the prior weight on the first entry of grid |
z_prior_func |
Prior function for z = arctanh(rho) |
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.