## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)

## ----message=FALSE, warning=FALSE---------------------------------------------
library(magrittr)
library(dplyr)

## ----message=FALSE, warning=FALSE---------------------------------------------
library(ConfidenceEllipse)

## -----------------------------------------------------------------------------
data(glass, package = "ConfidenceEllipse")

## ----message=FALSE, warning=FALSE---------------------------------------------
ellipsoid <- glass %>%
  confidence_ellipsoid(x = SiO2, y = Na2O, z = Fe2O3, conf_level = 0.95)

## -----------------------------------------------------------------------------
ellipsoid %>% glimpse()

## ----message=FALSE, warning=FALSE---------------------------------------------
rgl::setupKnitr(autoprint = TRUE)
rgl::plot3d(
  x = ellipsoid$x,
  y = ellipsoid$y,
  z = ellipsoid$z,
  xlab = "SiO2 (wt.%)",
  ylab = "Na2O (wt.%)",
  zlab = "Fe2O3 (wt.%)",
  type = "l",
  radius = .05,
  col = "darkgrey"
)
rgl::points3d(
  x = glass$SiO2,
  y = glass$Na2O,
  z = glass$Fe2O3,
  col = "darkred",
  size = 5
)
rgl::view3d(zoom = .8)

## ----message=FALSE, warning=FALSE---------------------------------------------
ellipsoid_grp <- glass %>%
  confidence_ellipsoid(x = SiO2, y = Na2O, z = Fe2O3, .group_by = glassType, conf_level = 0.95)

## -----------------------------------------------------------------------------
ellipsoid_grp %>% glimpse()

## ----message=FALSE, warning=FALSE---------------------------------------------
rgl::setupKnitr(autoprint = TRUE)
rgl::plot3d(
  x = ellipsoid_grp$x,
  y = ellipsoid_grp$y,
  z = ellipsoid_grp$z,
  xlab = "SiO2 (wt.%)",
  ylab = "Na2O (wt.%)",
  zlab = "Fe2O3 (wt.%)",
  type = "s",
  radius = .03,
  col = as.numeric(ellipsoid_grp$glassType)
)
rgl::points3d(
  x = glass$SiO2,
  y = glass$Na2O,
  z = glass$Fe2O3,
  col = as.numeric(glass$glassType),
  size = 5
)
rgl::view3d(zoom = .8)