fit_AR1_Gaussian function

Fit Gaussian AR(1) model to time series with missing values and/or outliers