gas_model <- UniGASSpec(Dist = "std", ScalingType = "Identity",
                GASPar = list(location = TRUE, scale = TRUE, shape = FALSE))

gas_fit <- UniGASFit(gas_model, time_series)

gas_forecast <- UniGASFor(gas_fit, H = 1, ReturnDraws =  T)

prediction_point <- mean(gas_forecast@Draws)