[R-br] Gráfico Pareto

[Elisa Henning]

Boa tarde

A função paretoPlot do pacote pid, produz um gráfico com cores diferentes
para os efeitos positivos e negativos.
Att
Elisa

Em 18 de fevereiro de 2016 13:20, Ari Clecius <ari072000 em yahoo.com.br>
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> 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
> 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
> 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/fad7a1df/attachment.html>