a = rnorm(100000,60000,1000)
b = a
c = data.frame(a,b)
color <- c("black", "red")
color_1 <- sample(color, nrow(c), replace=TRUE, prob=c(0.9, 0.1))
c$color_1 = as.factor(color_1)
hist(c$a, col = c$color_1, 100000)
But there are no colors that are showing up. Does anyone know how to fix this?