Alexandre, bom dia!
Adotando a mesma ideia utilizada no caso da Fátima, dá pra fazer algo bem
simples, mas que pode servir.
### <BEGIN>
require(gplots)
# plot means and
data(state)
### pra apagar o que não vai usar:
rm(list=ls(patt='^state')[!(ls() %in% c('state.area', 'state.region'))])
tmp <- split(state.area, state.region)
means <- sapply(tmp, mean)
stdev <- sqrt(sapply(tmp, var))
n <- sapply(tmp,length)
ciw <- qt(0.975, n) * stdev / sqrt(n)
# plain
plotCI(x=means, uiw=ciw)
###
rbind(means,stdev,n,ciw) ### sumário
### envelope simples
plotCI(x=means, uiw=ciw)
eixoX <- 1:4
lines(eixoX, means+ciw, lty='dashed', col=2)
lines(eixoX, means-ciw, lty='dashed', col=2)
### envelope simples hachurado
plotCI(x=means, uiw=ciw)
polX <- c(eixoX,rev(eixoX))
polY <- c(means+ciw,rev(means-ciw))
#colors()[grep("light",colors())] ### pesquisando cores
#polygon(polX, polY, lty='dashed', col='lightgray', bor=NA) ### sem borda
ao redor do polígono
polygon(polX, polY, lty='dashed', col='lightgray', bor=2) ### com borda
vermelha
plotCI(x=means, uiw=ciw, add=T)
### <END>
Éder Comunello <c <comunello.e...@gmail.com>omunello.e...@gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
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.
