Изменение названия отдельных графиков на панели графиков в ggplot2 / rstan на греческие буквы с индексом

Question

Я создаю трассировки моих объектов stanfit через mcmc_trace.Я хочу переименовать заголовки трасс.Мне уже удалось изменить заголовок, но я не знаю, как переименовать графики в греческие буквы с индексами, подобными функции выражения.

array <- as.array(fit)
array[1,1,1:3]
dimnames(array)[[3]][1:3] <- c("alpha1", "alpha2", "alpha3")

trace <- mcmc_trace(array, pars = c("alpha1", "alpha2", "alpha3")

Я хочу заменить alpha1 выражением (gamma [0]), но это не работает.

Answer 1

Хорошо, поэтому после небольшого копания просто кажется, что bayesplott::mcmc_trace не имеет возможности для греческого письма.

Однако, у очень похожего пакета mcmcplots::traplot есть опция greek = TRUE.Я сделал случайный пример для вас, чтобы увидеть:

library(mcmcplots)

nc <- 3; nr <- 1000
pnames <- c(paste('alpha[', 1:2, ']', sep = ''), paste('gamma[1]', sep = ''))
means <- rpois(10, 20)
fakemcmc <- coda::as.mcmc.list(
  lapply(1:3, function(i) coda::mcmc(matrix(rnorm(nc*nr, rep(means, each=nr)),
                                            nrow=nr, dimnames=list(NULL,pnames)))))
traplot(fakemcmc, greek = TRUE)

, который производит этот вывод: