This is a section

for (i in 1:3) {
  dat.sim <- data.frame(x = rnorm(1000))
  pandoc.header(sprintf("Subsection %i", i), level = 2)
  plot(ggplot(dat.sim, aes(x)) + geom_density())
  pandoc.p("")
}

Subsection 1

Subsection 2

Subsection 3