getPosterior | R Documentation |
Either the patient level data or both mu_hat as well as sd_hat must to be provided. If patient level data is provided mu_hat and se_hat are calculated within the function using a linear model. This function calculates the posterior for every dose group independently via the RBesT function postmix().
getPosterior(
prior_list,
data = NULL,
mu_hat = NULL,
se_hat = NULL,
calc_ess = FALSE
)
prior_list |
a prior list with information about the prior to be used for every dose group |
data |
dataframe containing the information of dose and response. Default NULL Also a simulateData object can be provided. |
mu_hat |
vector of estimated mean values (per dose group). |
se_hat |
vector of estimated standard deviations (per dose group). |
calc_ess |
boolean variable, indicating whether effective sample size should be calculated. Default FALSE |
posterior_list, a posterior list object is returned with information about (mixture) posterior distribution per dose group
prior_list <- list(Ctrl = RBesT::mixnorm(comp1 = c(w = 1, m = 0, s = 5), sigma = 2),
DG_1 = RBesT::mixnorm(comp1 = c(w = 1, m = 1, s = 12), sigma = 2),
DG_2 = RBesT::mixnorm(comp1 = c(w = 1, m = 1.2, s = 11), sigma = 2) ,
DG_3 = RBesT::mixnorm(comp1 = c(w = 1, m = 1.3, s = 11), sigma = 2) ,
DG_4 = RBesT::mixnorm(comp1 = c(w = 1, m = 2, s = 13), sigma = 2))
mu <- c(0, 1, 1.5, 2, 2.5)
se <- c(5, 4, 6, 7, 8)
posterior_list <- getPosterior(
prior_list = prior_list,
mu_hat = mu,
se_hat = se)
summary(posterior_list)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.