[R-br] Gráfico Pareto

[salah]

Caro,

segue sugestão:

## dados
dd = data.frame(valor = 2*abs(f1), grupo = names(f1))

## equação
lb1 = paste("S[ef]%.%t[crit]==", s)

library(ggplot2)

## gráfico
p = ggplot(dd, aes(grupo, valor, fill=grupo)) +
     theme_bw() +
     geom_bar(stat = "identity") +
     xlab(NULL) +
     ylab(NULL) +
     ggtitle("Gráfico de pareto dos efeitos") +
     geom_text(aes(label = round(valor, 1)),
             position = position_dodge(width = 0.9),
             angle = -90,
             vjust = -0.25) +
     annotate("text", x = 2, y = 10, label = lb1, parse=TRUE) +
     coord_flip()

p

ggsave("grafico_pareto.jpg", p, width = 6, height = 4)

saudações

Em 18/02/2016 13:20, Ari Clecius escreveu:
> Prezados colegas, este assunto me interessa, gostaria de saber como 
> automatizar este gráfico, fiz um script, mas é todo manual, os 
> gráficos de pareto dos pacotes não formatam por cor, seria possível 
> visualizar os efeitos negativos e positivos em cores diferentes?
>
> Att,
> Ari Clecius Alves de Lima
> Engenheiro Químico
> Me. Engenharia Civil
> Dr. Engenharia Civil
> (085)84174333
> (085)33669042
>
> */#############3/*
> */#A=temperatura/*
> */#B=vibração/*
> */#A<-c(20,20,60,60,20,20,60,60,20,20,60,60)/*
> */#B<-c(5,15,5,15,5,15,5,15,5,15,5,15)/*
> */#y<-c(74,50,46,34,68,48,44,32,65,51,40,26)/*
> */
> /*
> */A<-c(-1,-1,1,1,-1,-1,1,1,-1,-1,1,1)/*
> */B<-c(-1,1,-1,1,-1,1,-1,1,-1,1,-1,1)/*
> */y<-c(74,50,46,34,68,48,44,32,65,51,40,26)/*
> */group<-c("A","B","C","D","A","B","C","D","A","B","C","D")/*
> */falha<-data.frame(A,B,y,group)/*
> */aov.2=aov(y~A+B+A*B, data=falha)/*
> */summary(aov.2)/*
> */
> /*
> */
> /*
> */efeitos<-2*(aov.2$coef[2:4])/*
> */
> /*
> */#Estimativa da variância do erro experimental (s^2)/*
> */grau_L1<-2/*
> */grau_L2<-2/*
> */grau_L3<-2/*
> */grau_L4<-2/*
> */ag <- data.frame(aggregate(. ~ group, falha, function(x) c(mean = 
> mean(x), sd = sd(x))))/*
> */s2<-(sum(2*(ag$y[,2])^2))/(grau_L1+grau_L2+grau_L3+grau_L4)/*
> */
> /*
> */#syi^2=s2/n onde n é o número de replicatas(3)/*
> */#sy1^2=sy2^2=sy3^2=sy4^2=syq/*
> */#sy=raiz(syq)/*
> */n<-3/*
> */sy=sqrt(s2/n)/*
> */
> /*
> */#t(ef)=ef/s(ef)/*
> */tefA=efeitos[1]/sy/*
> */tefB=efeitos[2]/sy/*
> */tefAB=efeitos[3]/sy/*
> */
> /*
> */
> /*
> */#Gráfico de pareto/*
> */
> /*
> */f1<-c(tefA,tefB,tefAB);f1/*
> */z=data.frame(2*abs(f1))/*
> */z=t(z)/*
> */v1=sy*(qt(0.975,8))/*
> */coef3=sort(z)/*
> */coef3/*
> */
> /*
> */#png(file = "pareto_ceto.png",res=300, width = 4, height = 4, units 
> = 'in')/*
> */g=barplot(sort(z),horiz=TRUE,names.arg=c( 
> rev(names(f1))),col=ifelse(sort(abs(f1))>v1,"blue","gray"),xlim=c(0,25))##xlim, 
> ylim/*
> */box("plot", col="blue") /*
> */abline(v=v1,col="red",lty=2)/*
> */s=format(v1,digits=2)/*
> */text(6,1,label=expression("S"[ef]*" x t"[crit]*"="),cex=1)/*
> */text(9,1,label=s,cex=1)/*
> */mtext("Gráfico de pareto dos efeitos", side=3, line=2, cex=1.0, 
> col="blue", outer=FALSE)/*
> */p=format(f1,scientific = TRUE)/*
> */legend(x=16,y=2,cex=0.6,bg="yellow",box.col="blue", legend=c( 
> "A=-22.3","B=-16.0","AB=3.33" ), col="blue", lwd=1:2, lty=c(NA,NA,NA))/*
> */#dev.off()/*
>
>
> Em Quinta-feira, 18 de Fevereiro de 2016 10:18, Paulo Dick 
> <paulopcdick em gmail.com> escreveu:
>
>
> Bom dia Ana Paula,
>
> Seu gráfico não aparece para mim.
>
> Sugiro também que envie a sintaxe que utilizou.
>
> Abraço
>
> *
> *
> *Paulo Dick*
> Estatístico
> Mestre em Epidemiologia em Saúde Pública
> Tel.: (55 21) 99591-2716
>
> Em 17 de fevereiro de 2016 21:59, ana paula coelho madeira 
> <apcmadeira em hotmail.com <mailto:apcmadeira em hotmail.com>> escreveu:
>
>     Prezados,
>
>     boa noite!
>
>     Construí um gráfico de Pareto utilizando a função paretoPlot e a
>     saída é da forma:
>
>     Gostaria de fazer uma das seguintes alterações:
>
>     1) Apresentar as barras na horizontal;
>     ou
>     2) Colocar as estimativas na vertical ( para não ficar tão difícil
>     de visualizar).
>
>     Desde já agradeço.
>
>     Att.,
>
>     Ana Paula
>
>     _______________________________________________
>     R-br mailing list
>     R-br em listas.c3sl.ufpr.br <mailto:R-br em 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.
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br <mailto:R-br em 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� 
> c�igo m�imo reproduz�el.
>
>
>
>
> _______________________________________________
> R-br mailing list
> R-br em 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.

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160218/7ff60e24/attachment.html>