library(zoo)
m <- factor(months(zoo::as.yearmon(time(datasets::sunspots))),
            levels=month.name)
plot(m, datasets::sunspots, axes=FALSE)
axis(2)
axis(1, at=1:12, labels=FALSE)


plot(m, datasets::sunspots)