dev/test-vectors3d.R

In friendly/matlib: Matrix Functions for Teaching and Learning Linear Algebra and Multivariate Statistics

---
  title: "Test plot math in vectors3d"
---
  
library(rgl)
library(matlib)


# problem with `col =` in current version when more than one vector is given

open3d()
E <- diag(3)
rownames(E) <- c("x1", "y", "x2")
vectors3d(E, lwd=2)
planes3d(0, 0, 1, 0, col="gray",
         alpha=0.1)
# rotate
V <- matrix(c(1, 1, 1, -1, -1, -1), 2, 3, byrow=TRUE)
rownames(V) <- c("plus", "minus")
vectors3d(V, lwd=2, col = c("blue", "red"))


# want to use labels that are expressions.
# text3d() now allows this, using `plotmath3d()` for expressions
labs <- c(expression(x[1]), "y", expression(x[2]))
is.expression(labs)

#source(here::here("dev", "vectors3d.R"))



open3d()
E <- diag(3)
rownames(E) <- c(expression(x[1]), "y", expression(x[2]))
vectors3d(E, lwd=2)
vectors3d(c(1, 1, 1),
          labels=c(expression(hat(y))), color="red",
          lwd=3)
vectors3d(c(1, 1, 0),
          labels=c("", "x+y"),
          color="green", lwd=2)
planes3d(0, 0, 1, 0, col="gray",
         alpha=0.1)
segments3d(rbind(c(1, 1, 1),
                 c(1, 1, 0)))
arc(c(1, 1, 1), c(0, 0, 0),
    c(1, 1, 0))
corner(c(0, 0, 0), c(1, 1, 0), c(1, 1, 1))