[R-br] Rodar um exemplo com distribuição normal

[Mauro Sznelwar]

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()
}
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120120/4fcf8758/attachment.html>