#Codigo
mx=0
my=0
s2x=1
s2y=1
corxy=0
x=seq(-3,3,length=50)
y=x
z=outer(x, y, function(x,y){
(1/(2*pi*sqrt(s2x*s2y*(1-corxy^2))))*
exp(-1/(2*(1-corxy^2))*
((x-mx)^2/s2x +
(y-my)^2/s2y -
2*corxy*((x-mx)*(y-my))/
(sqrt(s2x)*sqrt(s2y))
)
)
}
)
z1=outer(x, y, function(x,y){x*y*sin(x^2+y^2)})
for(i in 1:1439){
persp(x, y, z, theta=i/4, phi=i/8,
expand = 0.5,col = rainbow(i),
ltheta = 120, shade = 0.75,
ticktype = "detailed",
xlab = "X", ylab = "Y", zlab = "Z")
title(sub="Funcion de densiodad de probabilidad")
title(main = "Distribucion normal Bivariada")
}
Si hay problemas en Rstudio, probar directamente en R
###################################################
mx=0
my=0
s2x=1
s2y=1
corxy=0
x=seq(-3,3,length=50)
y=x
z=outer(x, y, function(x,y){
(1/(2*pi*sqrt(s2x*s2y*(1-corxy^2))))*
exp(-1/(2*(1-corxy^2))*
((x-mx)^2/s2x +
(y-my)^2/s2y -
2*corxy*((x-mx)*(y-my))/
(sqrt(s2x)*sqrt(s2y))
)
)
}
)
z1=outer(x, y, function(x,y){x*y*sin(x^2+y^2)})
for(i in 1:1439){
persp(x, y, z, theta=i/4, phi=i/8,
expand = 0.5,col = rainbow(i),
ltheta = 120, shade = 0.75,
ticktype = "detailed",
xlab = "X", ylab = "Y", zlab = "Z")
title(sub="Funcion de densiodad de probabilidad")
title(main = "Distribucion normal Bivariada")
}
###################################################