G R A F I C O S    3 D

# O Generar Gif

###################################################

install.packages("rgl")

install.packages("magick")

install.packages("MASS")

library("rgl")
library("magick")

library(MASS)

# from the fitdistr example

set.seed(123)

x <- rgamma(100, shape = 5, rate = 0.1)

fit <- fitdistr(x, dgamma, list(shape = 1, rate = 0.1), lower = 0.001)

loglik <- function(shape, rate) sum(dgamma(x, shape=shape, rate=rate,

                                           log=TRUE))

loglik <- Vectorize(loglik)

xlim <- fit$estimate[1]+4*fit$sd[1]*c(-1,1)

ylim <- fit$estimate[2]+4*fit$sd[2]*c(-1,1)

mfrow3d(1, 2, sharedMouse = TRUE)

persp3d(loglik,

        xlim = xlim, ylim = ylim,

        n = 30)

zlim <- fit$loglik + c(-qchisq(0.99, 2)/2, 0)

next3d()

persp3d(loglik,

        xlim = xlim, ylim = ylim, zlim = zlim,

        n = 30)

play3d( spin3d( axis = c(0, 0, 1), rpm = 7), duration = 100 )

###################################################