library(fpp2)
## Loading required package: ggplot2
## Loading required package: forecast
## Warning: package 'forecast' was built under R version 3.5.2
## Loading required package: fma
## Loading required package: expsmooth
gas <- read.csv("montly-av-residential-gas-usage-.csv")
gas.ts <- ts(gas$Montly.av..residential.gas.usage.Iowa..cubic.feet..100..71....79, start = c(1971,1), frequency = 12 )
autoplot(gas.ts)

fit <- nnetar(gas.ts)
autoplot(forecast(fit,h=12))