def METHOD_NAME(
    mapper: Mapping[str, xr.Dataset],
    initial_time: str,
    limit_quantiles: Mapping[str, float],
    limit_only: Optional[Sequence[str]] = None,
) -> DatasetQuantileLimiter:
    sample_tendencies = mapper[initial_time]
    limiter = DatasetQuantileLimiter(
        limit_quantiles["upper"], limit_quantiles["lower"], limit_only=limit_only,
    )
    logger.debug(f"Fitting dataset limiter with limits={limit_quantiles}")
    return limiter.fit(sample_tendencies, feature_dims=["z", "tile"])