Bom dia senhores! Estou fazendo uma ilustração animada para demonstrar que a média é influenciada por valores extremos com o seguinte código: png(file='exemplos/med0%1d.png',width=500,height=250) x1 <- list() medias <- list() for(i in 1:10){ x <- c(7,10,13.4,15.3,17.6,19.5,20.6,21.1) y <- rep(0,8) x1[[i]] <- c(x,i+21.1) medias[[i]] <- mean(x1[[i]]) plot(x,y,ylim=c(0,5),xlim=c(0,35),ylab='',xlab='',axes=F) axis(1) par(xpd=TRUE) points(mean(x),-0.2,col='blue',lwd=10,pch=2) points(x1[[i]][9],0,col='red',lwd=7,cex=0.1) points(medias[[i]],-0.2,col='blue',lwd=10,pch=2) } dev.off() png(file='exemplos/med1%1d.png',width=500,height=250) x2 <- list() medias2 <- list() for(j in seq(10,1,-1)){ x2[[j]] <- c(12.3-j,x) medias2[[j]] <- mean(x2[[j]]) plot(x,y,ylim=c(0,5),xlim=c(0,35),ylab='',xlab='',axes=F) axis(1) par(xpd=TRUE) points(mean(x),-0.2,col='blue',lwd=10,pch=2) points(x1[[j]][9],0,col='red',lwd=7,cex=0.1) points(medias[[j]],-0.2,col='blue',lwd=10,pch=2) } dev.off() png(file='exemplos/med2%1d.png',width=500,height=250) for(j in 1:10){ x2[[j]] <- c(12.3-j,x) medias2[[j]] <- mean(x2[[j]]) plot(x,y,ylim=c(0,5),xlim=c(0,35),ylab='',xlab='',axes=F) axis(1) par(xpd=TRUE) points(mean(x),-0.2,col='blue',lwd=10,pch=2) points(x2[[j]][1],0,col='red',lwd=7,cex=0.1) points(medias2[[j]],-0.2,col='blue',lwd=10,pch=2) } dev.off() É possível juntar isso em um único código? Desde já grato por qualquer ajuda. (S,f,P) Allaman \begin{signature} <<>>= Prof. Dr. Ivan Bezerra Allaman Universidade Estadual de Santa Cruz Departamento de Ciências Exatas e Tecnológicas Ilhéus/BA - Brasil Fone: +55 73 3680-5596 E-mail: ivanalaman@yahoo.com.br/ivanalaman@gmail.com @ \end{signature}