[R-br] Plotar intervalo de confiança

[Éder Comunello]

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.eder em gmail.com>omunello.eder em gmail.com>
Dourados, MS - [22 16.5'S, 54 49'W]
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20131128/b19e0f3c/attachment.html>