assertThrows(IllegalStateException.class, EMPTY_PAIRED_STATS::pearsonsCorrelationCoefficient);
    assertThrows(
        IllegalStateException.class, ONE_VALUE_PAIRED_STATS::pearsonsCorrelationCoefficient);
    PairedStats infiniteXStats = createSingleStats(Double.POSITIVE_INFINITY, 1.23);
    assertThrows(IllegalStateException.class, infiniteXStats::pearsonsCorrelationCoefficient);

    assertThrows(IllegalStateException.class, EMPTY_PAIRED_STATS::pearsonsCorrelationCoefficient);
    assertThrows(
        IllegalStateException.class, ONE_VALUE_PAIRED_STATS::pearsonsCorrelationCoefficient);
    PairedStats infiniteXStats = createSingleStats(Double.POSITIVE_INFINITY, 1.23);
    assertThrows(IllegalStateException.class, infiniteXStats::pearsonsCorrelationCoefficient);

    assertThrows(
        IllegalStateException.class, () -> emptyAccumulator.pearsonsCorrelationCoefficient());
    assertThrows(
        IllegalStateException.class,
        () -> emptyAccumulatorByAddAllEmptyPairedStats.pearsonsCorrelationCoefficient());
    assertThrows(
        IllegalStateException.class, () -> oneValueAccumulator.pearsonsCorrelationCoefficient());
    assertThrows(
        IllegalStateException.class,

   *
   * @throws IllegalStateException if the dataset is empty or contains a single pair of values, or
   *     either the {@code x} and {@code y} dataset has zero population variance
   */
  public final double pearsonsCorrelationCoefficient() {
    checkState(count() > 1);
    if (isNaN(sumOfProductsOfDeltas)) {
      return NaN;
    }
    double xSumOfSquaresOfDeltas = xStats.sumOfSquaresOfDeltas();

   *
   * @throws IllegalStateException if the dataset is empty or contains a single pair of values, or
   *     either the {@code x} and {@code y} dataset has zero population variance
   */
  public double pearsonsCorrelationCoefficient() {
    checkState(count() > 1);
    if (isNaN(sumOfProductsOfDeltas)) {
      return NaN;
    }
    double xSumOfSquaresOfDeltas = xStats().sumOfSquaresOfDeltas();