mcmc_Rhat {sirt} | R Documentation |
Computation of the Rhat Statistic from a Single MCMC Chain
Description
Computes the Rhat statistic from a single MCMC chain.
Usage
mcmc_Rhat(mcmc_object, n_splits=3)
Arguments
mcmc_object |
Object of class |
n_splits |
Number of splits for MCMC chain |
Value
Numeric vector
Examples
## Not run:
#############################################################################
# EXAMPLE 1: Computation Rhat statistic for 2PNO model fitting by MCMC
#############################################################################
data(data.read)
# estimate 2PNO with MCMC with 3000 iterations and 500 burn-in iterations
mod <- sirt::mcmc.2pno( dat=data.read, iter=1000, burnin=100 )
# plot MCMC chains
plot( mod$mcmcobj, ask=TRUE )
# compute Rhat statistics
round( sirt::mcmc_Rhat( mod$mcmcobj[[1]] ), 3 )
## End(Not run)
[Package sirt version 4.1-15 Index]