I am trying to make plots in a loop. But how do I put different titles on each plot? In this example, I want different names for my 8 density plots, such as beta[Treatment], beta[Time Dummy], etc. Thanks!
par(mfrow=c(4,2)
for (i in 2:8) {
plot(density(beta[,i]))
title(main=substitute(paste('Density of ', beta[Treatment]))))
}