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 σ2a, σ2p, σ2ep, σ2e. 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