faltou uma leitura mais atenta do polygon()... media <- 170 stder <- 8 x <- seq(media-3*stder, media+3*stder, .01) y <- dnorm(x, media, stder) plot(x, y, type='l') abline(h=0, v=media, lty=2)
## parte inf hx <- seq(media-3*stder, media-2*stder, .01) hy <- dnorm(hx, media, stder) n <- length(hy) polygon(c(hx, rev(hx)), c(hy, rep(0, n)), col=2) ##parte sup 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) 2011/4/1 Augusto Ribas <[email protected]>: > Olá pessoal. > Bem eu estava lendo um livro aqui, o The R Book do M. J Crawley(2007) > E na pagina 212 tem uma figura que to tentando reproduzir mas num to > conseguindo. > A ideia é mais ou menos o seguinte. > > #Primeiro temos la um plot de uma curva de distribuição normal com 3 > desvios: > > x<-seq(-3,3,0.01) > y<-exp(-abs(x)^2) > plot(x,y,type="l") > > #Dai se no exemplo dele ele fala que se trata de dados de alturas de > pessoas, onde a media é 170cm e desvio padrão de 8cm > #dai eu mudei o eixo X da seguinte forma > > altura<-170+(8*x) > plot(altura,y,type="l") > > #dai o grafico fica com a a distribuição baseado no dados em questão altura > e nao em unidades de desvio. > #então chega o meu problema, ele nao deixo codigo de como faze o pintadinho > de ver a chande de alguma coisa. > #por exemplo, alturas menores de 162, qual a chance de ver alguem menor de > 162. > #eu vi naquelas listas de graficos que tem codigo que da pra marcar isso com > o comando polygon > #mas eu nao to conseguindo usalo direito pra isso, e no help dele os > exemplos nao me ajudaram muito > #algo assim > > plot(altura,y,type="l") > polygon(c(146,162,162),pnorm(c(-3,-1,-3)),col="black") > > #mas eu nao quero fazer um triangulo, tinha que ser aquela curvinha seguindo > os valores de y, mas ai como ter varios > #valores de y ao longo dos de x pra fica bonitinho que eu nao to > conseguindo, acho que olhando a figura da pra entender > #minha duvida > > Obrigado a todos. > -- > Grato > Augusto C. A. Ribas > > Site Pessoal: http://augustoribas.heliohost.org > Lattes: http://lattes.cnpq.br/7355685961127056 > > > _______________________________________________ > R-br mailing list > [email protected] > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > > _______________________________________________ R-br mailing list [email protected] https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
