Вывод gelman.diag по умолчанию форматирует данные так, чтобы они отображали 2 десятичных знака. Однако вы можете сохранить выходные данные gelman.diag в объекте и получить более непосредственный доступ к данным, как в следующем минимальном и воспроизводимом примере:
library(coda)
data(line)
ret <- gelman.diag(line)
# The output of gelman.diag automatically formats data to show 2 decimal places
ret
#Potential scale reduction factors:
#
# Point est. Upper C.I.
#alpha 1.02 1.02
#beta 1.00 1.00
#sigma 1.04 1.12
#
#Multivariate psrf
#
#1.01
# ret$psrf contains the full data
ret$psrf
# Point est. Upper C.I.
#alpha 1.019377 1.019838
#beta 1.000695 1.002321
#sigma 1.037599 1.115930
В целом, для большей гибкости при форматировании чисел с плавающей запятой, смотрите ссылку, предоставленную @NelsonGon в комментариях выше.