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)

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