aggregate_number: Aggregate replicate simulations
In rscherrer/egstools:

Aggregate a statistic over replicates by counting how many have values above a threshold, within values of one or two parameters.

aggregate_number(
  d,
  theta = 0.9,
  xname = "hsymmetry",
  yname = "ecosel",
  zname = "summary",
  colname = NULL,
  keep = NULL
)

`d`	Data frame with coordinates per timepoint per simulation
`theta`	The threshold to pass
`xname`	Column name of the first grouping parameter
`yname`	Column name of the second grouping parameter (set to NULL if only one grouping parameter)
`zname`	Column name of the statistic to aggregate
`colname`	Optional name of the aggregated statistic column
`keep`	Optional names of columns to keep in the aggregated data frame