Below is code to load in patients data, which is monthly data for 12 years of the number of patients seen at a dentist office. Suggest a MA model for this data. To fit a moving average model of order m the code is arima(patients,order=c(0,0,m)). In addition, provide and explain a forecast of the model using the plot() function to the fitted model.

patients <- as.ts(read.csv("https://sites.google.com/site/bsherwood/patients.csv",header=T)$patients)