getFinalState: Extract posterior draws from last iteration

In jacsor/MPK: Mixtures of Perturbed Kernels

Description

Extract posterior draws from last iteration

Usage

getFinalState(object, iter = object$mcmc$nsave)

Arguments

object	An MPK object
iter	Integer indicating which iteration to extract. By default, the last iteration is extracted.

Value

A list with the values of the MCMC at the iter-th iteration.

Examples

n = c(250, 250)
p = 4

Y1 = rbind( matrix( rnorm( n[1]*p), ncol = p), matrix( rnorm(n[2]*p) + 3, ncol = p))
Y2 = rbind( matrix( rnorm( n[1]*p), ncol = p), matrix( rnorm(n[2]*p) + 4, ncol = p))
Y = rbind(Y1, Y2)
C = c( rep(1,sum(n)), rep(2,sum(n)))

ans = mpk(Y, C)
state = getFinalState(ans)
ans2 = mpk(Y, C, state = state)

jacsor/MPK documentation built on May 18, 2019, 9:05 a.m.