Create tmIndMarkReg.R

dtm -> topic proportions + dates -> topic proportion time series -> plot(fitted(lm(tpts_1 ~ tpts_2))) ... confidence intervals, etc.
-> plot(fitted(lm(tpts_5 ~ rolling_avg_scrap_metal))), etc.