diff --git a/sources.c b/sources.c
index f86f198..094e21c 100644
--- a/sources.c
+++ b/sources.c
@@ -891,8 +891,6 @@ SRC_ReloadSources(void)
       LOG(LOGS_WARN, LOGF_Sources, "Could not open dump file %s", filename);
     } else {
       if (SST_LoadFromFile(sources[i]->stats, in)) {
-        /* We might want to use SST_DoUpdateRegression here, but we
-           need to check it has the same functionality */
         SST_DoNewRegression(sources[i]->stats);
       } else {
         LOG(LOGS_WARN, LOGF_Sources, "Problem loading from file %s", filename);
diff --git a/sourcestats.c b/sourcestats.c
index b338d37..4268fdf 100644
--- a/sourcestats.c
+++ b/sourcestats.c
@@ -468,43 +468,6 @@ SST_DoNewRegression(SST_Stats inst)
 
 }
 
-/* ================================================== */
-/* This function does a simple regression on what is in the register,
-   without trying to optimise the error bounds on the frequency by
-   deleting old samples */
-
-void
-SST_DoUpdateRegression(SST_Stats inst)
-{
-  double times_back[MAX_SAMPLES];
-  double freq_error_bound;
-  double est_intercept, est_slope, est_var_base, est_intercept_sd, est_slope_sd;
-
-  convert_to_intervals(inst, times_back);
-
-  if (inst->n_samples >= 3) { /* Otherwise, we're wasting our time - we
-                                 can't do a useful linear regression
-                                 with less than 3 points */
-    
-    RGR_WeightedRegression(times_back, inst->offsets, inst->weights,
-                           inst->n_samples,
-                           &est_intercept, &est_slope, &est_var_base,
-                           &est_intercept_sd, &est_slope_sd);
-
-    freq_error_bound = est_slope_sd * RGR_GetTCoef(inst->n_samples - 2);
-
-    inst->estimated_frequency = est_slope;
-    inst->skew = freq_error_bound;
-
-  } else {
-    inst->estimated_frequency = 0.0;
-    inst->skew = WORST_CASE_FREQ_BOUND;
-  }
-
-  find_best_sample_index(inst, times_back);
-
-}
-
 /* ================================================== */
 
 void