# Plot posterior distributions

Source: `R/src_summaries.R`
`plot_hdpglm.Rd`

this function creates a plot with two grids. One is the grid with posterior expectation of betas as function of context-level covariates. The other is the posterior distribution of tau

    plot_hdpglm(
      samples,
      X = NULL,
      W = NULL,
      ncol.taus = 1,
      ncol.betas = NULL,
      ncol.w = NULL,
      nrow.w = NULL,
      smooth.line = FALSE,
      pred.pexp.beta = FALSE,
      title.tau = NULL,
      true.tau = NULL,
      title.beta = NULL,
      tau.x.axis.size = 1.1,
      tau.y.axis.size = 1.1,
      tau.title.size = 1.2,
      tau.panel.title.size = 1.4,
      tau.legend.size = 1,
      beta.x.axis.size = 1.1,
      beta.y.axis.size = 1.1,
      beta.title.size = 1.2,
      beta.panel.title.size = 1.4,
      beta.legend.size = 1,
      tau.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.

ncol.taus

  integer with the number of columns of the grid containing the posterior distribution of tau

ncol.betas

  integer with the number of columns of the posterior expectation of betas as function of context-level features

ncol.w

  integer with the number of columns to use to display the different context-level covariates

nrow.w

  integer with the number of rows to use to display the different context-level covariates

smooth.line

  boolean, if `TRUE` the plot will display a regression line representing the regression of the posterior expectation of the linear coefficients betas on the context-level covariates. Default `FALSE`

pred.pexp.beta

  boolean, if `TRUE` the plots will display a line with the predicted posterior expectation of betas obtained using the posterior expectation of taus, the linear coefficients of the expectation of beta

title.tau

  string, the title for the posterior distribution of the context effects

true.tau

  a `data.frame` with four columns. The first must be named `w` and it indicates the index of each context-level covariate, starting with 0 for the intercept term. The second column named `beta` must contain the indexes of the betas of individual-level covariates, starting with 0 for the intercept term. The third column named `Parameter` must be named `tau<w><beta>`, where `w` and `beta` must be the actual values displayed in the columns `w` and `beta`. Finally, it must have a column named `True` with the true value of the parameter.

title.beta

  string, the title for the posterior expectation of beta as function of context-level covariate

tau.x.axis.size

  numeric, relative size of the x-axis of the plot with tau

tau.y.axis.size

  numeric, relative size of the y-axis of the plot with tau

tau.title.size

  numeric, relative size of the title of the plot with tau

tau.panel.title.size

  numeric, relative size of the title of the panels of the plot with tau

tau.legend.size

  numeric, relative size of the legend of the plot with tau

beta.x.axis.size

  numeric, relative size of the x-axis of the plot with beta

beta.y.axis.size

  numeric, relative size of the y-axis of the plot with beta

beta.title.size

  numeric, relative size of the title of the plot with beta

beta.panel.title.size

  numeric, relative size of the title of the panels of the plot with beta

beta.legend.size

  numeric, relative size of the legend of the plot with beta

tau.xlab

  string, the label of the x-axis for the plot with tau

## Examples

    library(magrittr)
    # 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_hdpglm(samples)
    #>
    #>
    #> Plot being generated ...
    #>
    #>
    #>
    #> Generating plots ...
    #>
    #> Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.

    plot_hdpglm(samples, ncol.taus=2, ncol.betas=2, X='x1')
    #>
    #>
    #> Plot being generated ...
    #>
    #>
    #>
    #> Generating plots ...
    #>
    #> Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.

    plot_hdpglm(samples, ncol.taus=2, ncol.betas=2, X='x1', ncol.w=2, nrow.w=1,
                pred.pexp.beta=TRUE,smooth.line=TRUE )
    #>
    #>
    #> Plot being generated ...
    #>
    #>
    #>
    #> Generating plots ...
    #>
    #> Warning: `guides(<scale> = FALSE)` is deprecated. Please use `guides(<scale> = "none")` instead.
    #> `geom_smooth()` using formula 'y ~ x'
    #> `geom_smooth()` using formula 'y ~ x'