Function to plot posterior distribution of tau
plot_tau(
samples,
X = NULL,
W = NULL,
title = NULL,
true.tau = NULL,
show.all.taus = FALSE,
show.all.betas = FALSE,
ncol = NULL,
legend.position = "top",
x.axis.size = 1.1,
y.axis.size = 1.1,
title.size = 1.2,
panel.title.size = 1.4,
legend.size = 1,
xlab = NULL
)Arguments
- samples
an output of the function
hdpGLM- X
a string vector with the name of the first-level covariates whose associated tau should be displayed
- W
a string vector with the name of the context-level covariate(s) whose linear effect will be displayed. If
NULL, the linear effect tau of all context-level covariates are displayed. Note: the context-level covariate must have been included in the estimation of the model.- title
string, title of the plot
- true.tau
a
data.framewith four columns. The first must be namedwand it indicates the index of each context-level covariate, starting with 0 for the intercept term. The second column namedbetamust contain the indexes of the betas of individual-level covariates, starting with 0 for the intercept term. The third column namedParametermust be namedtau<w><beta>, wherewandbetamust be the actual values displayed in the columnswandbeta. Finally, it must have a column namedTruewith the true value of the parameter.- show.all.taus
boolean, if
FALSE(default) the posterior distribution of taus representing the intercept of the expectation of beta are omitted- show.all.betas
boolean, if
FALSE(default) the taus affecting only the intercept terms of the outcome variable are omitted- ncol
number of columns of the grid. If
NULL, one column is used- legend.position
one of four options: "bottom" (default), "top", "left", or "right". It indicates the position of the legend
- x.axis.size
numeric, the relative size of the label in the x-axis
- y.axis.size
numeric, the relative size of the label in the y-axis
- title.size
numeric, the relative size of the title of the plot
- panel.title.size
numeric, the relative size of the titles in the panel of the plot
- legend.size
numeric, the relative size of the legend
- xlab
string, the label of the x-axis
Examples
library(magrittr)
set.seed(66)
# Note: this example is just for illustration. MCMC iterations are very reduced
set.seed(10)
n = 20
data.context1 = tibble::tibble(x1 = rnorm(n, -3),
x2 = rnorm(n, 3),
z = sample(1:3, n, replace=TRUE),
y =I(z==1) * (3 + 4*x1 - x2 + rnorm(n)) +
I(z==2) * (3 + 2*x1 + x2 + rnorm(n)) +
I(z==3) * (3 - 4*x1 - x2 + rnorm(n)) ,
w = 20
)
data.context2 = tibble::tibble(x1 = rnorm(n, -3),
x2 = rnorm(n, 3),
z = sample(1:2, n, replace=TRUE),
y =I(z==1) * (1 + 3*x1 - 2*x2 + rnorm(n)) +
I(z==2) * (1 - 2*x1 + x2 + rnorm(n)),
w = 10
)
data = data.context1 %>%
dplyr::bind_rows(data.context2)
## estimation
mcmc = list(burn.in=1, n.iter=50)
samples = hdpGLM(y ~ x1 + x2, y ~ w, data=data, mcmc=mcmc, n.display=1)
#>
#>
#> Preparing for estimation ...
#>
#>
#>
#> Estimation in progress ...
#>
#>
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 1
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 1
#> Maximum Number of Clusters Active in the Current Iteration : 1
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000e+00
#> 1.0000e+02
#>
#> Clusters in Context 2
#> 1.0000
#> 95.0000
#>
#> [=== ] 3 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 2
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 2
#> Maximum Number of Clusters Active in the Current Iteration : 2
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000
#> 25.0000 75.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 4.0000
#> 10.0000 55.0000 35.0000
#>
#> [===== ] 5 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 3
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 4
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 5.0000
#> 15.0000 10.0000 40.0000 35.0000
#>
#> Clusters in Context 2
#> 1.0000 2.0000 3.0000 4.0000 7.0000
#> 10.0000 30.0000 30.0000 10.0000 15.0000
#>
#> [====== ] 7 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 4
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 6
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 2.0000 4.0000 5.0000 7.0000
#> 10.0000 35.0000 10.0000 35.0000
#>
#> Clusters in Context 2
#> 3.0000 4.0000 5.0000 7.0000
#> 35.0000 15.0000 20.0000 15.0000
#>
#> [======= ] 9 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 5
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 8
#> Maximum Number of Clusters Active in the Current Iteration : 8
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 4.0000 5.0000 6.0000
#> 10.0000 35.0000 10.0000 10.0000 20.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 40.0000
#>
#> [========= ] 11 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 6
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 9
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 2.0000 3.0000 5.0000 6.0000 7.0000 8.0000
#> 15.0000 10.0000 10.0000 15.0000 10.0000 25.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 40.0000
#>
#> [========== ] 13 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 7
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 8
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 2.0000 3.0000 4.0000 5.0000 6.0000 11.0000
#> 20.0000 10.0000 15.0000 15.0000 20.0000 10.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 40.0000
#>
#> [=========== ] 15 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 8
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 7
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 4.0000 5.0000 6.0000
#> 15.0000 10.0000 15.0000 25.0000 25.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 45.0000
#>
#> [============= ] 17 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 9
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 8
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 2.0000 4.0000 5.0000 6.0000
#> 10.0000 20.0000 10.0000 40.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 45.0000
#>
#> [============== ] 19 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 10
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 9
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 2.0000 4.0000 5.0000 6.0000
#> 15.0000 15.0000 10.0000 35.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 40.0000
#>
#> [================ ] 21 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 11
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 8
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 4.0000 6.0000 7.0000 10.0000
#> 15.0000 15.0000 35.0000 10.0000 10.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 50.0000
#>
#> [================= ] 23 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 12
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 8
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 4.0000 6.0000 7.0000 9.0000
#> 15.0000 40.0000 15.0000 10.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 50.0000
#>
#> [================== ] 25 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 13
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 8
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 4.0000 6.0000 7.0000
#> 10.0000 15.0000 45.0000 10.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 50.0000
#>
#> [==================== ] 27 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 14
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 8
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 4.0000 6.0000 7.0000
#> 10.0000 10.0000 15.0000 40.0000 10.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 45.0000 45.0000
#>
#> [===================== ] 29 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 15
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 5
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 4.0000 6.0000
#> 20.0000 10.0000 15.0000 50.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 45.0000
#>
#> [====================== ] 31 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 16
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 4.0000 6.0000
#> 20.0000 15.0000 50.0000
#>
#> Clusters in Context 2
#> 3.0000 7.0000
#> 50.0000 45.0000
#>
#> [======================== ] 33 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 17
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 4.0000 5.0000 6.0000
#> 10.0000 10.0000 10.0000 15.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 7.0000
#> 10.0000 50.0000 40.0000
#>
#> [========================= ] 35 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 18
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 6.0000
#> 10.0000 20.0000 10.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 7.0000
#> 10.0000 50.0000 40.0000
#>
#> [=========================== ] 37 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 19
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 7
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 6.0000 7.0000
#> 10.0000 15.0000 10.0000 45.0000 10.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 6.0000 7.0000
#> 10.0000 50.0000 10.0000 30.0000
#>
#> [============================ ] 39 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 20
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 7
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 5.0000 6.0000
#> 15.0000 10.0000 15.0000 45.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 7.0000
#> 10.0000 50.0000 35.0000
#>
#> [============================= ] 41 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 21
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 4.0000 6.0000
#> 20.0000 10.0000 10.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 7.0000
#> 10.0000 45.0000 40.0000
#>
#> [=============================== ] 43 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 22
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 6.0000
#> 20.0000 15.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 7.0000
#> 10.0000 50.0000 40.0000
#>
#> [================================ ] 45 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 23
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 5
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 6.0000
#> 25.0000 10.0000 10.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 7.0000
#> 15.0000 50.0000 35.0000
#>
#> [================================= ] 47 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 24
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 5
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000 6.0000
#> 20.0000 10.0000 10.0000 10.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 7.0000
#> 15.0000 50.0000 35.0000
#>
#> [=================================== ] 49 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 25
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 5
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 6.0000
#> 25.0000 15.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 7.0000
#> 25.0000 40.0000 25.0000
#>
#> [==================================== ] 50 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 26
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 6.0000
#> 25.0000 10.0000 15.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000 7.0000
#> 35.0000 45.0000 15.0000
#>
#> [====================================== ] 52 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 27
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 6.0000
#> 25.0000 10.0000 15.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 45.0000 45.0000
#>
#> [======================================= ] 54 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 28
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 4.0000 6.0000
#> 20.0000 15.0000 10.0000 45.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [======================================== ] 56 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 29
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 6.0000
#> 20.0000 10.0000 10.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 45.0000 50.0000
#>
#> [========================================== ] 58 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 30
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 6.0000
#> 15.0000 10.0000 15.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [=========================================== ] 60 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 31
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 5
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 6.0000
#> 25.0000 15.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [============================================ ] 62 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 32
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 5
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 6.0000
#> 20.0000 20.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [============================================== ] 64 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 33
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 6.0000
#> 15.0000 20.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [=============================================== ] 66 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 34
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 6.0000
#> 25.0000 20.0000 50.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 45.0000
#>
#> [================================================= ] 68 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 35
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 6.0000
#> 25.0000 25.0000 45.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [================================================== ] 70 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 36
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 6.0000
#> 25.0000 25.0000 45.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 45.0000 50.0000
#>
#> [=================================================== ] 72 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 37
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 5
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 6.0000
#> 20.0000 25.0000 45.0000
#>
#> Clusters in Context 2
#> 1.0000 2.0000 3.0000
#> 40.0000 10.0000 50.0000
#>
#> [===================================================== ] 74 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 38
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 3.0000 5.0000 6.0000
#> 20.0000 30.0000 10.0000 40.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [====================================================== ] 76 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 39
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 5
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 6.0000
#> 20.0000 10.0000 25.0000 40.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [======================================================= ] 78 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 40
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 7
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000 6.0000
#> 10.0000 20.0000 15.0000 15.0000 30.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 45.0000 50.0000
#>
#> [========================================================= ] 80 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 41
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 5
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000
#> 15.0000 45.0000 20.0000 15.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [========================================================== ] 82 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 42
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000 6.0000
#> 15.0000 35.0000 20.0000 15.0000 10.0000
#>
#> Clusters in Context 2
#> 1.0000 2.0000 3.0000
#> 45.0000 10.0000 45.0000
#>
#> [============================================================ ] 84 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 43
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000
#> 15.0000 40.0000 20.0000 15.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [============================================================= ] 86 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 44
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 6
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000
#> 10.0000 45.0000 20.0000 15.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [============================================================== ] 88 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 45
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000
#> 10.0000 50.0000 20.0000 20.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [================================================================ ] 90 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 46
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000
#> 15.0000 50.0000 20.0000 15.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 45.0000
#>
#> [================================================================= ] 92 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 47
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000
#> 10.0000 50.0000 20.0000 20.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [================================================================== ] 94 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 48
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000 4.0000
#> 15.0000 50.0000 25.0000 10.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 45.0000 50.0000
#>
#> [==================================================================== ] 96 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 49
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000
#> 25.0000 50.0000 20.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [===================================================================== ] 98 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 50
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000
#> 20.0000 50.0000 25.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 50.0000 50.0000
#>
#> [=======================================================================] 100 %
#>
#> -----------------------------------------------------
#> MCMC in progress ....
#>
#> Family of the distribution of the outcome variable of the mixture components: gaussian
#>
#> Burn-in: 1
#> Number of MCMC Iterations : 50
#>
#> Iteration: 51
#>
#> Acceptance Rate for beta : 1
#> Average Acceptance Rate for beta : 1
#>
#> Maximum Number of Clusters Allowed (K) : 100
#> Maximum Number of Clusters Activated Among All Contexts : 9
#> Maximum Number of Clusters Active in the Current Iteration : 4
#>
#> (displaying only clusters with more than 5% of the data)
#> Clusters in Context 1
#> 1.0000 2.0000 3.0000
#> 25.0000 50.0000 20.0000
#>
#> Clusters in Context 2
#> 1.0000 3.0000
#> 45.0000 50.0000
#>
#> [=======================================================================] 100 %
plot_tau(samples)
plot_tau(samples, ncol=2)
plot_tau(samples, X='x1', W='w')
plot_tau(samples, show.all.taus=TRUE, show.all.betas=TRUE, ncol=2)
