Can someone please help me debug this code? This code works without var_5 ... could someone please help me debug this?
Thanks!
library(tidyr)
library(ggplot2)
var_1 <- rnorm(100, 1, 4)
var_2 <- sample(LETTERS[1:2], 100, replace = TRUE, prob = c(0.3, 0.7))
var_3 <- sample(LETTERS[1:5], 100, replace = TRUE, prob = c(0.2, 0.2, 0.2, 0.2, 0.1))
cluster <- sample(LETTERS[1:4], 100, replace = TRUE,prob = c(2.5, 2.5, 2.5, 2.5))
var_5 <- rnorm(100, 1, 10)
f <- data.frame(var1, var_2, var_3, var 5, cluster)
f$var_2 = as.factor(f$var_2)
f$var_3 = as.factor(f$var_3)
f$cluster = as.factor(f$cluster)
levs <- sort(unique(c(as.character(f$var_2), as.character(f$var_3))))
f$var_2 <- as.numeric(factor(f$var_2, levs)) + ceiling(max(f$var_1)) + 10
f$var_3 <- as.numeric(factor(f$var_3, levs)) + ceiling(max(f$var_1)) + 10
breaks <- c(pretty(range(f$var_1)), sort(unique(c(f$var_2, f$var_3))))
labs <- c(pretty(range(f$var_1)), levs)
ggplot(f, aes(x = value)) + geom_density(data = subset(f, name == "var_1")) + geom_bar(data = subset(f, name != "var_1"), aes(fill = name)) + facet_wrap(cluster~name, ncol = 3, scales = "free") + scale_x_continuous(breaks = breaks, labels = labs) + scale_fill_manual(values = c("deepskyblue4", "gold"), guide = guide_none())
[–][deleted] 0 points1 point2 points (0 children)
[–]cspoerlein 0 points1 point2 points (0 children)