Tenho esta situação:
 
f=rmultinom(30,3,prob=c(0.5,0.3,0.8))*dmultinom(x=1:3,prob=c(0.5,0.3,0.8))
 n=5
 for(i in 1:n){
 g=f
contour(g)
}
#preciso que gere 5 gráficos diferentes (mas não parece que gera gráficos #diferentes)
jpeg(filename = "MeusGraficos%03d.jpg", width = 480, height = 480)
for (i in 1:5) contour(g)
dev.off()
 
E como adapto o caso da normal para o multinomial para fazer gráfico 3D, já tentei mas acusa medidas diferentes.
z <- matrix(rnorm(50*10),50,10)
image(z)
persp(z)