Skip to content

Instantly share code, notes, and snippets.

@benilton
Created August 15, 2011 10:55
Show Gist options
  • Save benilton/1146031 to your computer and use it in GitHub Desktop.
Save benilton/1146031 to your computer and use it in GitHub Desktop.
Sobreposição de Gráficos Usando polygon()
## Por Benilton Carvalho
##
## Determinando os parâmetros da Normal
media <- 170
stder <- 8
## Criando o grid no qual deseja-se visualizar a densidade
x <- seq(media-3*stder, media+3*stder, .01)
## Determinando a densidade em cada ponto do grid
y <- dnorm(x, media, stder)
## Gráfico da densidade no grid de interesse
plot(x, y, type='l', xlab='Eixo X', ylab='Eixo Y')
abline(h=0, v=media, lty=2)
## Determinando o sub-grid de interesse
## e a densidade em cada ponto da regiao
## Região inferior apenas.
hx <- seq(media-3*stder, media-2*stder, .01)
hy <- dnorm(hx, media, stder)
## A funcao polygon() necessita do caminho inteiro
## da regiao a ser colorida. Então, para cada hx[i],
## existem 2 pontos determinando o polígono: a) na
## altura da densidade e b) em zero (para fechar o
## polígono).
n <- length(hy)
polygon(c(hx, rev(hx)), c(hy, rep(0, n)), col=2)
## Repetindo para a região superior
hx <- seq(media+2*stder, media+3*stder, .01)
hy <- dnorm(hx, media, stder)
n <- length(hy)
polygon(c(hx, rev(hx)), c(hy, rep(0, n)), col=2)
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment