Please help in below R programming:

This exercise uses the canadian_gas data (monthly Canadian gas production in billions of cubic metres, January 1960 - February 2005). The values are in the column Volume.

1. Please help in Plot the data using autoplot(), gg_subseries() and gg_season().

Example : autoplot(cangas) + xlab("Month") + ylab("Gas Production (billions of cubic metres)") + ggtitle("Canadian Gas Production (Jan 1960 - Feb 2005)")

2. Please help in STL decomposition of the data. Please help in how to plot a new time series fit that will contain STL model (you do not need any additional parameters for STL), then plot it

3. Please help in How does the seasonal shape change over time? [Hint: Try plotting the seasonal component using gg_season().] To find the name for the seasonal component, we can view fit or look at the Environment panel (top left).