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<-2grau_L2<-2grau_L3<-2grau_L4<-2ag <- 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<-3sy=sqrt(s2/n) #t(ef)=ef/s(ef)tefA=efeitos[1]/sytefB=efeitos[2]/sytefAB=efeitos[3]/sy #Gráfico de pareto f1<-c(tefA,tefB,tefAB);f1z=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, ylimbox("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@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 DickEstatísticoMestre em Epidemiologia em Saúde PúblicaTel.: (55 21) 99591-2716 Em 17 de fevereiro de 2016 21:59, ana paula coelho madeira <apcmadeira@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@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@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.

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@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@hotmail.com <mailto:apcmadeira@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@listas.c3sl.ufpr.br <mailto: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.

_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br <mailto: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� c�igo m�imo reproduz�el.

_______________________________________________ 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.