Skip to contents

Process ngvb.list when method = "Gibbs". Produces a matrix with the samples of the mixing variables \(\mathbf{V}\).

Usage

# S4 method for ngvb.list
mungeGibbs.mixing(object)

Arguments

object

An ngvb.list object (output of ngvb function)

Examples

if (FALSE) {
 #Here we fit an RW1 latent process to the jumpts time series
 plot(jumpts)

 #Fit LGM with INLA
 LGM     <- inla(y ~ -1 + f(x,  model = "rw1"),
                 data = jumpts)

 #Run 10 iterations from the Gibbs sampler
 LnGM.Gibbs <- ngvb(fit = LGM, selection = list(x=1:100),
                    method = "Gibbs", iter=10)

 Gibbs.V    <- mungeGibbs.mixing(LnGM.Gibbs)
 Gibbs.eta  <- mungeGibbs.ng(LnGM.Gibbs)
 }