plot.CZ_ANOVA: Diagnostic graphics proposed by...

In TPMarques/RJCA: Convergence assessment of RJMCMC algorithms

Diagnostic graphics proposed by \insertCiteCZ;textualconvRJMCMC for a single parameter

Description

This function plots the 3 diagnostic graphics proposed by \insertCiteCZ;textualconvRJMCMC for a single parameter: 1 - PSRF1 and PSRF2 plus their upper bounds against batch; 2 - V and Wc against batch; 3 - Wm and WmWc against batch.

Usage

## S3 method for class 'CZ_ANOVA'
plot(x,...)

Arguments

x	Object of class CZ_ANOVA created by CZ_ANOVA function.
...	Extra arguments to be parsed in the function.

Details

The plots are created in the currently active device. Users must press enter to see other plots.

Note

Axis and legend labels are in Brazilian Portuguese in this preliminary version. Future versions will have English labels as standard and allow user-defined labels.

Author(s)

Tiago Pereira Marques and Miriam Harumi Tsunemi

References

\insertAllCited

Examples

#Example 1

set.seed(634963)
theta1<-rnorm(100000,mean=20,sd=1.5)
theta2<-rnorm(100000,mean=20,sd=1.5)
theta3<-rnorm(100000,mean=20,sd=1.5)
theta<-c(theta1,theta2,theta3)
chains<-c(rep(1,100000),rep(2,100000),rep(3,100000))
models<-factor(sample(c(1,2,3),300000,TRUE,c(0.25,0.5,0.25)))
mcmciterations<-rep(seq(1:100000),3)

#CZ_ANOVA object is saved on test_1
teste_1<-CZ_ANOVA(theta,chains,models,mcmciterations)
plot(teste_1)

TPMarques/RJCA documentation built on Feb. 22, 2024, 8:05 p.m.