Location Dispersion Ellipsoids

> TwoDimEllipsoid <- function(sample.mean, sample.cov, Scale) {
+     x1 <- sample.mean[1] - sqrt(sample.cov[1, 1]) * Scale
+     x2 <- sample.mean[1] + sqrt(sample.cov[1, 1]) * Scale
+     y1 <- sample.mean[2] + sqrt(sample.cov[2, 2]) * Scale
+     y2 <- sample.mean[2] - sqrt(sample.cov[2, 2]) * Scale
+     xrange <- c(range(x[, 1])[1] - sample.mean[1], range(x[,
+         1])[2] + sample.mean[1])
+     yrange <- c(range(x[, 2])[1] - sample.mean[2], range(x[,
+         2])[2] + sample.mean[2])
+     par(mfrow = c(1, 1))
+     plot(x, pch = 19, col = "blue", xlab = "X", ylab = "Y", main = "",
+         xlim = xrange, ylim = yrange)
+     vecs <- eigen(sample.cov)$vectors
+     evals <- eigen(sample.cov)$values
+     length1 <- sqrt(evals[1])
+     length2 <- sqrt(evals[2])
+     segments(sample.mean[1], sample.mean[2], sample.mean[1] +
+         vecs[1, 1] * length1 * Scale, sample.mean[2] + vecs[2,
+         1] * length1 * Scale, lwd = 3, col = "green")
+     segments(sample.mean[1], sample.mean[2], sample.mean[1] +
+         vecs[1, 2] * length2 * Scale, sample.mean[2] + vecs[2,
+         2] * length2 * Scale, lwd = 3, col = "green")
+     par(new = T)
+     angles <- seq(from = 0, to = 2 * pi, by = pi/500)
+     ang.points <- matrix(data = NA, nrow = length(angles), ncol = 2)
+     for (i in seq_along(angles)) {
+         ang.points[i, 1] <- (cos(angles[i])) * sqrt(evals[1]) *
+             Scale
+         ang.points[i, 2] <- (sin(angles[i])) * sqrt(evals[2]) *
+             Scale
+     }
+     ellipsoid <- matrix(data = NA, nrow = n, ncol = 2)
+     ellipsoid <- (ang.points %*% t(vecs))
+     ellipsoid[, 1] <- (ellipsoid[, 1]) + sample.mean[1]
+     ellipsoid[, 2] <- (ellipsoid[, 2]) + sample.mean[2]
+     plot(ellipsoid, xlab = "", ylab = "", xlim = xrange, ylim = yrange,
+         col = "red")
+     abline(v = sample.mean[1], lty = "dotted", col = "grey")
+     abline(h = sample.mean[2], lty = "dotted", col = "grey")
+ }

> n <- 10000
> sample.mean <- c(1, 1)
> sample.cov <- matrix(data = c(1, 0, 0, 1), nrow = 2, ncol = 2)
> x <- rmnorm(n = n, mean = sample.mean, varcov = sample.cov)
> Scale <- 3
> TwoDimEllipsoid(sample.mean, sample.cov, Scale)

> n <- 10000
> sample.mean <- c(1, 1)
> sample.cov <- matrix(data = c(1, 0.9, 0.9, 1), nrow = 2, ncol = 2)
> x <- rmnorm(n = n, mean = sample.mean, varcov = sample.cov)
> Scale <- 3
> TwoDimEllipsoid(sample.mean, sample.cov, Scale)

> n <- 10000
> sample.mean <- c(2, 2)
> sample.cov <- matrix(data = c(1, -0.9, -0.9, 1), nrow = 2, ncol = 2)
> x <- rmnorm(n = n, mean = sample.mean, varcov = sample.cov)
> Scale <- 3
> TwoDimEllipsoid(sample.mean, sample.cov, Scale)

> n <- 10000
> sample.mean <- c(2, 2)
> sample.cov <- matrix(data = c(1, 0.5, 0.5, 1), nrow = 2, ncol = 2)
> x <- rmnorm(n = n, mean = sample.mean, varcov = sample.cov)
> Scale <- 3
> TwoDimEllipsoid(sample.mean, sample.cov, Scale)

> n <- 10000
> sample.mean <- c(2, 2)
> sample.cov <- matrix(data = c(1, -0.5, -0.5, 1), nrow = 2, ncol = 2)
> x <- rmnorm(n = n, mean = sample.mean, varcov = sample.cov)
> Scale <- 3
> TwoDimEllipsoid(sample.mean, sample.cov, Scale)