library(coda)

Chain1=read.coda("Pars1.out","Pars.ind")
Chain2=read.coda("Pars2.out","Pars.ind")

#Convergence diagnostics
geweke.diag(Chain1) 
geweke.diag(Chain2)
heidel.diag(Chain1) 
heidel.diag(Chain2)
raftery.diag(Chain1) 
raftery.diag(Chain2)
gelman.diag(mcmc.list(Chain1,Chain2))

#Find 95% intervals of highest posterior density
#Separately
HPDinterval(mcmc.list(Chain1,Chain2),prob=0.95) 
#Both chains combined
HPDinterval(mcmc(rbind(Chain1,Chain2)),prob=0.95) 

#The above can also be done using the interactive codamenu() command