 
     %Suite{suite | scenarios: new_scenarios}
   end
-
-  defp iterations_per_second(0.0, _), do: 0.0
-
-  defp iterations_per_second(_, :memory), do: 0.0
-
-  defp iterations_per_second(average_microseconds, :run_time) do
-    Duration.microseconds({1, :second}) / average_microseconds
-  end
-
-  defp standard_deviation(_samples, _average, 1), do: 0
-
-  defp standard_deviation(samples, average, sample_size) do
-    total_variance =
-      Enum.reduce(samples, 0, fn sample, total ->
-        total + :math.pow(sample - average, 2)
-      end)
-
-    variance = total_variance / (sample_size - 1)
-    :math.sqrt(variance)
-  end
 end

  def statistics(suite = %Suite{scenarios: scenarios}) do
    scenarios_with_statistics =
      Parallel.map(scenarios, fn scenario ->
        run_time_stats = job_statistics(scenario.run_times, :run_time)
        memory_stats = job_statistics(scenario.memory_usages, :memory)

        %Scenario{
          scenario


    %__MODULE__{
      average: average,
      std_dev: deviation,
      std_dev_ratio: standard_dev_ratio,
      median: median,
      percentiles: percentiles,
      mode: mode,
