On Aug 30, 2011 4:47 AM, "Daniele Portela de oliveira" <dane_portela@yahoo.com.br> wrote:
> 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