#' Variance Explained
#'
#' @param x a \code{prcomp} object
#' @return a numeric vector of percentage variance explained for each principal componenet (PC)
#'
#' @export
variance_exp <- function(x)
{
if (class(x) != "prcomp") {
stop(deparse(substitute(x)), " must be a 'prcomp' object", call. = FALSE)
}
variance_explained <-
round(((x$sdev ^ 2) / sum(x$sdev ^ 2)) * 100, digits = 2)
return(as.numeric(variance_explained))
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.