Как использовать индексы в легендах ggplot2 [R]

Могу ли я использовать индексы в легендах ggplot2? Я вижу этот вопрос о греческих буквах в легендах и других местах, но я не могу понять, как его адаптировать.

Я подумал, что с помощью выражения () , которое работает в метки оси, сделают свое дело. Но моя попытка ниже терпит неудачу. Спасибо!

library(ggplot2)
temp <- data.frame(a = rep(1:4, each = 100), b = rnorm(4 * 100), c = 1 + rnorm(4 * 100))
names(temp)[2:3] <- c("expression(b[1])", "expression(c[1])")
temp.m <- melt(temp, id.vars = "a")
ggplot(temp.m, aes(x = value, linetype = variable)) + geom_density() + facet_wrap(~ a)