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
0.166804 0.150406 0.334731 0.651941 0.255857509
0.230704926
0.143072 0.076084 0.28832 0.507476 0.281928604
0.149926302
0.102943 0.102883 0.296767 0.502593 0.204823784
0.204704403
0.076466 0.132665 0.333456 0.542587 0.140928552
0.244504568
0.085353 0.146785 0.282928 0.515066 0.165712744
0.284982895
0.076825 0.156807 0.28215 0.515782 0.148948587
0.304017977
0.088717 0.130882 0.258809 0.478408 0.185442133
0.273578201
0.119562 0.155156 0.254589 0.529307 0.225884033
0.293130452
0.092786 0.210945 0.329216 0.632947 0.146593633
0.33327435
_______________________________________________
R-br mailing list
[email protected]
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código
mínimo reproduzível.
