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