Oi pessoal, eu consegui rodar este exemplo debaixo da distribuição multinomial frames = 5 for(i in 1:frames){ # x = seq(0, i, 1) f.3 = rmultinom(30,3,prob=c(0.5,0.3,0.8))*dmultinom(x=1:3,prob=c(0.5,0.3,0.8)) if (i <= 5) {name = paste('000',i,'plot.jpeg',sep='')} #saves the plot as a .jpeg file in the working directory jpeg(name) contour(f.3, ylab ='probability',main = paste('Multinomial density with n = ', i), col = 'red') dev.off() } Mas não consegui rodar este da normal, alguém sabe como fazê-lo? frames = 5 for(i in 1:frames){ x = seq(0, i, 1) f.4 = rnorm(30)*dnorm(1) if (i <= 5) {name = paste('000',i,'plot.jpeg',sep='')} #saves the plot as a .jpeg file in the working directory jpeg(name) dens= density(f.4, ylab ='probability',main = paste('Normal density with n = ', i), col = 'red') points(dens,type='l',col=1) # ou de preferência com o gráfico de contorno contour(f.4, ylab ='probability', main = paste('Normal density with n = ', i), col = 'red') dev.off() }