Nothing
ggbootvalid_variables <- function(resmca, axes = c(1,2), type = "partial", K = 30,
ellipse = "norm", level = 0.95, col = NULL, legend = "right") {
classe <- attr(resmca,'class')[1]
if(classe %in% c("MCA","speMCA","csMCA")) {
rate1 <- modif.rate(resmca)$modif$mrate[axes[1]]
rate2 <- modif.rate(resmca)$modif$mrate[axes[2]]
}
if(classe %in% c("stMCA","multiMCA","PCA")) {
rate1 <- modif.rate(resmca)$raw$rate[axes[1]]
rate2 <- modif.rate(resmca)$raw$rate[axes[2]]
}
noms <- getvarnames(resmca)
boot <- bootvalid_variables(resmca, axes = axes, type = type, K = K)
names(boot)[3:4] <- c("dimX", "dimY")
boot <- merge(boot, noms, by = "varcat")
coord <- data.frame(resmca$var$coord[, axes])
names(coord) <- c("dimX", "dimY")
coord$varcat <- rownames(coord)
coord <- merge(coord, noms, by = "varcat")
p <- ggplot2::ggplot()
if(is.null(col)) {
p <- p +
ggrepel::geom_text_repel(data = coord, ggplot2::aes(x = .data$dimX, y = .data$dimY, label = .data$varcat, col = .data$var), key_glyph='blank', parse = FALSE, max.overlaps = Inf) +
ggplot2::stat_ellipse(data = boot, ggplot2::aes(x = .data$dimX, y = .data$dimY, group = .data$varcat, col = .data$var), type = ellipse, level = level, lty = 2) +
ggplot2::geom_point(data = boot, ggplot2::aes(x = .data$dimX, y = .data$dimY, col = .data$var), alpha = 0.2)
} else if(is.character(col)) {
p <- p +
ggrepel::geom_text_repel(data = coord, ggplot2::aes(x = .data$dimX, y = .data$dimY, label = .data$varcat), color = col, key_glyph='blank', parse = FALSE, max.overlaps = Inf) +
ggplot2::stat_ellipse(data = boot, ggplot2::aes(x = .data$dimX, y = .data$dimY, group = .data$varcat), color = col, type = ellipse, level = level, lty = 2) +
ggplot2::geom_point(data = boot, ggplot2::aes(x = .data$dimX, y = .data$dimY), color = col, alpha = 0.2)
}
p <- p +
ggplot2::geom_hline(yintercept = 0, colour = "darkgrey", linewidth = .1) +
ggplot2::geom_vline(xintercept = 0, colour = "darkgrey", linewidth = .1) +
ggplot2::xlab(paste0("dim ", axes[1], " (", round(rate1,1), " %)")) +
ggplot2::ylab(paste0("dim ", axes[2], " (", round(rate2,1), " %)")) +
ggplot2::theme_bw() +
ggplot2::theme(panel.grid.major = ggplot2::element_blank(),
panel.grid.minor = ggplot2::element_blank()) +
ggplot2::guides(color = ggplot2::guide_legend(title = ""))
p
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.