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))
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)*
participantes (2)
-
Ari Clecius -
Elisa Henning