
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@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@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)*