plot.CZ_MANOVA: 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 multiple parameters: 1 - MPSRF1 and MPSRF2 against batch; 2 - Largest eigenvalue of V and largest eigenvalue of Wc against batch; 3 - Largest eigenvalue of Wm and largest eigenvalue of WmWc against batch.

Usage

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

Arguments

x	Object of class CZ_MANOVA created by CZ_MANOVA 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(500)

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)

{
theta1a<-rnorm(100000,mean=0,sd=1)
theta2a<-rnorm(100000,mean=1,sd=2)
theta3a<-rnorm(100000,mean=1,sd=3)
thetaa<-cbind(theta1a,theta2a,theta3a)

theta1b<-rnorm(100000,mean=0,sd=1)
theta2b<-rnorm(100000,mean=1,sd=2)
theta3b<-rnorm(100000,mean=1,sd=3)
thetab<-cbind(theta1b,theta2b,theta3b)

theta1c<-rnorm(100000,mean=0,sd=1)
theta2c<-rnorm(100000,mean=1,sd=2)
theta3c<-rnorm(100000,mean=1,sd=3)
thetac<-cbind(theta1c,theta2c,theta3c)

theta<-rbind(thetaa,thetab,thetac)
}

#CZ_MANOVA object is saved on test_2
teste_2<-CZ_MANOVA(theta,chains,models,mcmciterations)
plot(teste_2)

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