[R-br] Gráfico de Pareto
Elisa Henning
elisa.henning em gmail.com
Segunda Fevereiro 10 21:23:42 BRST 2014
Boa noite Ari
Talvez ajude olhar a função paretoPlot do pacote qualityTools
Elisa
2014-02-10 21:10 GMT-02:00 Ari Clecius <ari072000 em yahoo.com.br>:
> Alguém poderia me dizer uma forma mais fácil de fazer um gráfico de pareto
> para los efeitos significativos de um experimento fatorial? Segue código
> abaixo.
>
> ATtt,
>
> Ari Clecius Alves de Lima
> Engenheiro Químico
> Me. Engenharia Civil
> (085)88412345
> (085)33669042
>
> #Código
> *resp=scan(nmax=8)*
> *82.8*
> *55.1*
> *74.3*
> *72.2*
> *64.7*
> *51.3*
> *69.9*
> *53.3*
>
> *#replicated*
> *rad=rep(c(+1,-1),time=4,each=1)*
> *rad*
> *pH=rep(c(+1,-1),time=2,each=2)*
> *pH*
> *area=rep(c(+1,-1),each=4)*
> *area*
>
> *order=(1:8)*
> *order*
>
>
> *fat=data.frame(rad,pH,area,resp,order)*
> *fat*
> *aov.2=aov(resp~rad+pH+area, data=fat)*
> *summary(aov.2)*
>
> *#Gráfico de pareto*
>
> *MSE=41.7 #Obtido a partir de *
> *coef1=c(aov.2$coef[2:4])*
> *coef1*
> *n=8*
> *f1=t(abs(coef1))/(sqrt(MSE/n));f1*
> *z=data.frame(f1)*
> *names(z)*
> *v1=qt(0.975,8)*
> *coef3=sort(abs(coef1))*
> *coef3*
> *g=barplot(sort(f1),horiz=TRUE,names.arg=c(
> names(coef3)),col=ifelse(sort(f1)>v1,"blue","gray"))##xlim, ylim*
> *box("plot", col="blue") *
>
> *abline(v=qt(0.975,8),col="red",lty=2)*
> *s=format(qt(0.975,8),digits=2)*
> *text(2.5,1,label=s,cex=1)*
> *mtext("Efeito padronizado", side=3, line=2, cex=1.0, col="blue",
> outer=FALSE)*
> *p=format(f1,scientific = TRUE)*
> *legend(x=6,y=4,cex=0.6,bg="yellow",box.col="blue", legend=c(
> "Rad=7.475","area=5.650","pH=-1.975" ), col="blue", lwd=1:2,
> lty=c(NA,NA,NA))*
>
>
> _______________________________________________
> 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.
>
--
Profa. Elisa Henning
DMAT/CCT/UDESC
(47) 4009 7978
http://www.joinville.udesc.br/portal/professores/elisa/
*Um abismo. *
*Uma ponte a ser construída.*
*Não somente como um monumento, mas para permitir que outros prossigam com
menos dificuldade.*
*(Phillip Ross)*
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140210/83b45371/attachment.html>
Mais detalhes sobre a lista de discussão R-br