Caros colegas, mais uma vez eu recorro a vcs para resolver um problema com gráfico.

Gostaria que aparecessem os valores de média e intervalos de credibilidade em cada caixa;

segundo, no eixo x, aparecesse em forma de letra grega sigma^2 para representar os componentes

de variancia σ²a, σ²p, σ²ep, σ²e.

Abaixo tem um exemplo da minha tentativa mas não deu certo. Aceito sugestões!!!

#exemplo

require(coda)
mat.pn=read.table("cadeia.txt", header=F)
mat.pn=as.matrix(mat.pn) # falar que é uma matriz

y=mcmc(mat.pn)
vara= y[,1]
varp= y[,3]
varep= y[,4]
vare= y[,2]
library(grDevices)
x.1=cbind(vara,varp,varep,vare)
y.1=boxplot(x.1,range=0,ylab='Observações',font.main=3, font.lab=1,
main = "PN", cex.main = 1.2, font.main = 3)
axis(1, at = c(sigma^2(a), sigma^2(p),sigma^2(ep),sigma^2(e)),
labels = expression(sigma^2(a), sigma^2(p),sigma^2(ep),sigma^2(e)))

Abraços

Daniele Portela de Oliveira

Mestranda Produção Animal - UEM

Cel: 44 9949 3525

67 9617 0965