基本统计 - 基于 RDD 的 API

概要统计
相关性
分层抽样
假设检验
- 流式显著性检验
随机数据生成
核密度估计

\[ \newcommand{\R}{\mathbb{R}} \newcommand{\E}{\mathbb{E}} \newcommand{\x}{\mathbf{x}} \newcommand{\y}{\mathbf{y}} \newcommand{\wv}{\mathbf{w}} \newcommand{\av}{\mathbf{\alpha}} \newcommand{\bv}{\mathbf{b}} \newcommand{\N}{\mathbb{N}} \newcommand{\id}{\mathbf{I}} \newcommand{\ind}{\mathbf{1}} \newcommand{\0}{\mathbf{0}} \newcommand{\unit}{\mathbf{e}} \newcommand{\one}{\mathbf{1}} \newcommand{\zero}{\mathbf{0}} \]

概要统计

我们通过 Statistics 中提供的函数 colStats 为 RDD[Vector] 提供列概要统计信息。

colStats() 返回 MultivariateStatisticalSummary 的一个实例，其中包含按列排列的最大值、最小值、平均值、方差和非零值的数量，以及总计数。

有关 API 的更多详细信息，请参阅 MultivariateStatisticalSummary Python 文档。

import numpy as np

from pyspark.mllib.stat import Statistics

mat = sc.parallelize(
    [np.array([1.0, 10.0, 100.0]), np.array([2.0, 20.0, 200.0]), np.array([3.0, 30.0, 300.0])]
)  # an RDD of Vectors

# Compute column summary statistics.
summary = Statistics.colStats(mat)
print(summary.mean())  # a dense vector containing the mean value for each column
print(summary.variance())  # column-wise variance
print(summary.numNonzeros())  # number of nonzeros in each column

在 Spark 仓库中的 "examples/src/main/python/mllib/summary_statistics_example.py" 中查找完整的示例代码。

colStats() 返回 MultivariateStatisticalSummary 的一个实例，其中包含按列排列的最大值、最小值、平均值、方差和非零值的数量，以及总计数。

有关 API 的详细信息，请参阅 MultivariateStatisticalSummary Scala 文档。

import org.apache.spark.mllib.linalg.Vectors
import org.apache.spark.mllib.stat.{MultivariateStatisticalSummary, Statistics}

val observations = sc.parallelize(
  Seq(
    Vectors.dense(1.0, 10.0, 100.0),
    Vectors.dense(2.0, 20.0, 200.0),
    Vectors.dense(3.0, 30.0, 300.0)
  )
)

// Compute column summary statistics.
val summary: MultivariateStatisticalSummary = Statistics.colStats(observations)
println(summary.mean)  // a dense vector containing the mean value for each column
println(summary.variance)  // column-wise variance
println(summary.numNonzeros)  // number of nonzeros in each column

在 Spark 仓库中的 "examples/src/main/scala/org/apache/spark/examples/mllib/SummaryStatisticsExample.scala" 中查找完整的示例代码。

colStats() 返回 MultivariateStatisticalSummary 的一个实例，其中包含按列排列的最大值、最小值、平均值、方差和非零值的数量，以及总计数。

有关 API 的详细信息，请参阅 MultivariateStatisticalSummary Java 文档。

import java.util.Arrays;

import org.apache.spark.api.java.JavaRDD;
import org.apache.spark.mllib.linalg.Vector;
import org.apache.spark.mllib.linalg.Vectors;
import org.apache.spark.mllib.stat.MultivariateStatisticalSummary;
import org.apache.spark.mllib.stat.Statistics;

JavaRDD<Vector> mat = jsc.parallelize(
  Arrays.asList(
    Vectors.dense(1.0, 10.0, 100.0),
    Vectors.dense(2.0, 20.0, 200.0),
    Vectors.dense(3.0, 30.0, 300.0)
  )
); // an RDD of Vectors

// Compute column summary statistics.
MultivariateStatisticalSummary summary = Statistics.colStats(mat.rdd());
System.out.println(summary.mean());  // a dense vector containing the mean value for each column
System.out.println(summary.variance());  // column-wise variance
System.out.println(summary.numNonzeros());  // number of nonzeros in each column

在 Spark 仓库中的 "examples/src/main/java/org/apache/spark/examples/mllib/JavaSummaryStatisticsExample.java" 中查找完整的示例代码。

分层抽样

与其他位于 spark.mllib 中的统计函数不同，分层抽样方法 sampleByKey 和 sampleByKeyExact 可以在键值对的 RDD 上执行。对于分层抽样，可以将键视为标签，将值视为特定属性。例如，键可以是男人或女人，或文档 ID，而相应的值可以是人口中人的年龄列表或文档中的单词列表。 sampleByKey 方法将掷硬币以决定是否对观察值进行抽样，因此需要一次遍历数据，并提供预期样本大小。 sampleByKeyExact 比 sampleByKey 中使用的每层简单随机抽样需要更多的资源，但将以 99.99% 的置信度提供精确的抽样大小。 Python 目前不支持 sampleByKeyExact。

sampleByKey() 允许用户对大约 $\lceil f_k \cdot n_k \rceil \, \forall k \in K$ 个项目进行抽样，其中 $f_k$ 是键 $k$ 的所需比例，$n_k$ 是键 $k$ 的键值对数，$K$ 是键的集合。

注意： Python 目前不支持 sampleByKeyExact()。

# an RDD of any key value pairs
data = sc.parallelize([(1, 'a'), (1, 'b'), (2, 'c'), (2, 'd'), (2, 'e'), (3, 'f')])

# specify the exact fraction desired from each key as a dictionary
fractions = {1: 0.1, 2: 0.6, 3: 0.3}

approxSample = data.sampleByKey(False, fractions)

在 Spark 仓库中的 "examples/src/main/python/mllib/stratified_sampling_example.py" 中查找完整的示例代码。

sampleByKeyExact() 允许用户精确地对 $\lceil f_k \cdot n_k \rceil \, \forall k \in K$ 个项目进行抽样，其中 $f_k$ 是键 $k$ 的所需比例，$n_k$ 是键 $k$ 的键值对数，$K$ 是键的集合。无放回抽样需要额外一次遍历 RDD 以保证样本大小，而有放回抽样需要额外两次遍历。

// an RDD[(K, V)] of any key value pairs
val data = sc.parallelize(
  Seq((1, 'a'), (1, 'b'), (2, 'c'), (2, 'd'), (2, 'e'), (3, 'f')))

// specify the exact fraction desired from each key
val fractions = Map(1 -> 0.1, 2 -> 0.6, 3 -> 0.3)

// Get an approximate sample from each stratum
val approxSample = data.sampleByKey(withReplacement = false, fractions = fractions)
// Get an exact sample from each stratum
val exactSample = data.sampleByKeyExact(withReplacement = false, fractions = fractions)

在 Spark 仓库中的 "examples/src/main/scala/org/apache/spark/examples/mllib/StratifiedSamplingExample.scala" 中查找完整的示例代码。

import java.util.*;

import scala.Tuple2;

import org.apache.spark.api.java.JavaPairRDD;

List<Tuple2<Integer, Character>> list = Arrays.asList(
    new Tuple2<>(1, 'a'),
    new Tuple2<>(1, 'b'),
    new Tuple2<>(2, 'c'),
    new Tuple2<>(2, 'd'),
    new Tuple2<>(2, 'e'),
    new Tuple2<>(3, 'f')
);

JavaPairRDD<Integer, Character> data = jsc.parallelizePairs(list);

// specify the exact fraction desired from each key Map<K, Double>
ImmutableMap<Integer, Double> fractions = ImmutableMap.of(1, 0.1, 2, 0.6, 3, 0.3);

// Get an approximate sample from each stratum
JavaPairRDD<Integer, Character> approxSample = data.sampleByKey(false, fractions);
// Get an exact sample from each stratum
JavaPairRDD<Integer, Character> exactSample = data.sampleByKeyExact(false, fractions);

在 Spark 仓库中的 "examples/src/main/java/org/apache/spark/examples/mllib/JavaStratifiedSamplingExample.java" 中查找完整的示例代码。

假设检验

假设检验是统计学中一个强大的工具，用于确定结果是否具有统计意义，以及该结果是否偶然发生。 spark.mllib 目前支持 Pearson 的卡方 ($\chi^2$) 检验，用于拟合优度和独立性。输入数据类型决定是进行拟合优度检验还是独立性检验。拟合优度检验需要 Vector 的输入类型，而独立性检验需要 Matrix 作为输入。

spark.mllib 还支持输入类型 RDD[LabeledPoint]，以支持通过卡方独立性检验进行特征选择。

Statistics 提供了运行 Pearson 卡方检验的方法。以下示例演示了如何运行和解释假设检验。

有关 API 的更多详细信息，请参阅 Statistics Python 文档。

from pyspark.mllib.linalg import Matrices, Vectors
from pyspark.mllib.regression import LabeledPoint
from pyspark.mllib.stat import Statistics

vec = Vectors.dense(0.1, 0.15, 0.2, 0.3, 0.25)  # a vector composed of the frequencies of events

# compute the goodness of fit. If a second vector to test against
# is not supplied as a parameter, the test runs against a uniform distribution.
goodnessOfFitTestResult = Statistics.chiSqTest(vec)

# summary of the test including the p-value, degrees of freedom,
# test statistic, the method used, and the null hypothesis.
print("%s\n" % goodnessOfFitTestResult)

mat = Matrices.dense(3, 2, [1.0, 3.0, 5.0, 2.0, 4.0, 6.0])  # a contingency matrix

# conduct Pearson's independence test on the input contingency matrix
independenceTestResult = Statistics.chiSqTest(mat)

# summary of the test including the p-value, degrees of freedom,
# test statistic, the method used, and the null hypothesis.
print("%s\n" % independenceTestResult)

obs = sc.parallelize(
    [LabeledPoint(1.0, [1.0, 0.0, 3.0]),
     LabeledPoint(1.0, [1.0, 2.0, 0.0]),
     LabeledPoint(1.0, [-1.0, 0.0, -0.5])]
)  # LabeledPoint(label, feature)

# The contingency table is constructed from an RDD of LabeledPoint and used to conduct
# the independence test. Returns an array containing the ChiSquaredTestResult for every feature
# against the label.
featureTestResults = Statistics.chiSqTest(obs)

for i, result in enumerate(featureTestResults):
    print("Column %d:\n%s" % (i + 1, result))

在 Spark 仓库中的 "examples/src/main/python/mllib/hypothesis_testing_example.py" 中查找完整的示例代码。

Statistics 提供了运行 Pearson 卡方检验的方法。以下示例演示了如何运行和解释假设检验。

import org.apache.spark.mllib.linalg._
import org.apache.spark.mllib.regression.LabeledPoint
import org.apache.spark.mllib.stat.Statistics
import org.apache.spark.mllib.stat.test.ChiSqTestResult
import org.apache.spark.rdd.RDD

// a vector composed of the frequencies of events
val vec: Vector = Vectors.dense(0.1, 0.15, 0.2, 0.3, 0.25)

// compute the goodness of fit. If a second vector to test against is not supplied
// as a parameter, the test runs against a uniform distribution.
val goodnessOfFitTestResult = Statistics.chiSqTest(vec)
// summary of the test including the p-value, degrees of freedom, test statistic, the method
// used, and the null hypothesis.
println(s"$goodnessOfFitTestResult\n")

// a contingency matrix. Create a dense matrix ((1.0, 2.0), (3.0, 4.0), (5.0, 6.0))
val mat: Matrix = Matrices.dense(3, 2, Array(1.0, 3.0, 5.0, 2.0, 4.0, 6.0))

// conduct Pearson's independence test on the input contingency matrix
val independenceTestResult = Statistics.chiSqTest(mat)
// summary of the test including the p-value, degrees of freedom
println(s"$independenceTestResult\n")

val obs: RDD[LabeledPoint] =
  sc.parallelize(
    Seq(
      LabeledPoint(1.0, Vectors.dense(1.0, 0.0, 3.0)),
      LabeledPoint(1.0, Vectors.dense(1.0, 2.0, 0.0)),
      LabeledPoint(-1.0, Vectors.dense(-1.0, 0.0, -0.5)
      )
    )
  ) // (label, feature) pairs.

// The contingency table is constructed from the raw (label, feature) pairs and used to conduct
// the independence test. Returns an array containing the ChiSquaredTestResult for every feature
// against the label.
val featureTestResults: Array[ChiSqTestResult] = Statistics.chiSqTest(obs)
featureTestResults.zipWithIndex.foreach { case (k, v) =>
  println(s"Column ${(v + 1)} :")
  println(k)
}  // summary of the test

在 Spark 仓库中的 "examples/src/main/scala/org/apache/spark/examples/mllib/HypothesisTestingExample.scala" 中查找完整的示例代码。

Statistics 提供了运行 Pearson 卡方检验的方法。以下示例演示了如何运行和解释假设检验。

有关 API 的详细信息，请参阅 ChiSqTestResult Java 文档。

import java.util.Arrays;

import org.apache.spark.api.java.JavaRDD;
import org.apache.spark.mllib.linalg.Matrices;
import org.apache.spark.mllib.linalg.Matrix;
import org.apache.spark.mllib.linalg.Vector;
import org.apache.spark.mllib.linalg.Vectors;
import org.apache.spark.mllib.regression.LabeledPoint;
import org.apache.spark.mllib.stat.Statistics;
import org.apache.spark.mllib.stat.test.ChiSqTestResult;

// a vector composed of the frequencies of events
Vector vec = Vectors.dense(0.1, 0.15, 0.2, 0.3, 0.25);

// compute the goodness of fit. If a second vector to test against is not supplied
// as a parameter, the test runs against a uniform distribution.
ChiSqTestResult goodnessOfFitTestResult = Statistics.chiSqTest(vec);
// summary of the test including the p-value, degrees of freedom, test statistic,
// the method used, and the null hypothesis.
System.out.println(goodnessOfFitTestResult + "\n");

// Create a contingency matrix ((1.0, 2.0), (3.0, 4.0), (5.0, 6.0))
Matrix mat = Matrices.dense(3, 2, new double[]{1.0, 3.0, 5.0, 2.0, 4.0, 6.0});

// conduct Pearson's independence test on the input contingency matrix
ChiSqTestResult independenceTestResult = Statistics.chiSqTest(mat);
// summary of the test including the p-value, degrees of freedom...
System.out.println(independenceTestResult + "\n");

// an RDD of labeled points
JavaRDD<LabeledPoint> obs = jsc.parallelize(
  Arrays.asList(
    new LabeledPoint(1.0, Vectors.dense(1.0, 0.0, 3.0)),
    new LabeledPoint(1.0, Vectors.dense(1.0, 2.0, 0.0)),
    new LabeledPoint(-1.0, Vectors.dense(-1.0, 0.0, -0.5))
  )
);

// The contingency table is constructed from the raw (label, feature) pairs and used to conduct
// the independence test. Returns an array containing the ChiSquaredTestResult for every feature
// against the label.
ChiSqTestResult[] featureTestResults = Statistics.chiSqTest(obs.rdd());
int i = 1;
for (ChiSqTestResult result : featureTestResults) {
  System.out.println("Column " + i + ":");
  System.out.println(result + "\n");  // summary of the test
  i++;
}

在 Spark 仓库中的 "examples/src/main/java/org/apache/spark/examples/mllib/JavaHypothesisTestingExample.java" 中查找完整的示例代码。

此外，spark.mllib 提供了 Kolmogorov-Smirnov (KS) 检验的 1 样本、双侧实现，用于检验概率分布的相等性。通过提供理论分布的名称（目前仅支持正态分布）及其参数，或者提供根据给定的理论分布计算累积分布的函数，用户可以检验其样本是否从该分布中抽取的零假设。如果用户针对正态分布进行检验 (distName="norm")，但不提供分布参数，则检验将初始化为标准正态分布并记录适当的消息。

Statistics 提供了运行单样本、双侧 Kolmogorov-Smirnov 检验的方法。以下示例演示了如何运行和解释假设检验。

有关 API 的更多详细信息，请参阅 Statistics Python 文档。

from pyspark.mllib.stat import Statistics

parallelData = sc.parallelize([0.1, 0.15, 0.2, 0.3, 0.25])

# run a KS test for the sample versus a standard normal distribution
testResult = Statistics.kolmogorovSmirnovTest(parallelData, "norm", 0, 1)
# summary of the test including the p-value, test statistic, and null hypothesis
# if our p-value indicates significance, we can reject the null hypothesis
# Note that the Scala functionality of calling Statistics.kolmogorovSmirnovTest with
# a lambda to calculate the CDF is not made available in the Python API
print(testResult)

在 Spark 仓库的 "examples/src/main/python/mllib/hypothesis_testing_kolmogorov_smirnov_test_example.py" 中查找完整的示例代码。

Statistics 提供了运行单样本、双侧 Kolmogorov-Smirnov 检验的方法。以下示例演示了如何运行和解释假设检验。

有关 API 的详细信息，请参阅 Statistics Scala 文档。

import org.apache.spark.mllib.stat.Statistics
import org.apache.spark.rdd.RDD

val data: RDD[Double] = sc.parallelize(Seq(0.1, 0.15, 0.2, 0.3, 0.25))  // an RDD of sample data

// run a KS test for the sample versus a standard normal distribution
val testResult = Statistics.kolmogorovSmirnovTest(data, "norm", 0, 1)
// summary of the test including the p-value, test statistic, and null hypothesis if our p-value
// indicates significance, we can reject the null hypothesis.
println(testResult)
println()

// perform a KS test using a cumulative distribution function of our making
val myCDF = Map(0.1 -> 0.2, 0.15 -> 0.6, 0.2 -> 0.05, 0.3 -> 0.05, 0.25 -> 0.1)
val testResult2 = Statistics.kolmogorovSmirnovTest(data, myCDF)
println(testResult2)

在 Spark 仓库的 "examples/src/main/scala/org/apache/spark/examples/mllib/HypothesisTestingKolmogorovSmirnovTestExample.scala" 中查找完整的示例代码。

Statistics 提供了运行单样本、双侧 Kolmogorov-Smirnov 检验的方法。以下示例演示了如何运行和解释假设检验。

有关 API 的详细信息，请参阅 Statistics Java 文档。

import java.util.Arrays;

import org.apache.spark.api.java.JavaDoubleRDD;
import org.apache.spark.mllib.stat.Statistics;
import org.apache.spark.mllib.stat.test.KolmogorovSmirnovTestResult;

JavaDoubleRDD data = jsc.parallelizeDoubles(Arrays.asList(0.1, 0.15, 0.2, 0.3, 0.25));
KolmogorovSmirnovTestResult testResult =
  Statistics.kolmogorovSmirnovTest(data, "norm", 0.0, 1.0);
// summary of the test including the p-value, test statistic, and null hypothesis
// if our p-value indicates significance, we can reject the null hypothesis
System.out.println(testResult);

在 Spark 仓库的 "examples/src/main/java/org/apache/spark/examples/mllib/JavaHypothesisTestingKolmogorovSmirnovTestExample.java" 中查找完整的示例代码。

流式显著性检验

spark.mllib 提供了某些检验的在线实现，以支持 A/B 测试等用例。这些检验可以在 Spark Streaming DStream[(Boolean, Double)] 上执行，其中每个元组的第一个元素指示对照组 (false) 或处理组 (true)，第二个元素是观察值。

流式显著性检验支持以下参数

peacePeriod - 要忽略的流中的初始数据点数，用于缓解新奇效应。
windowSize - 执行假设检验的过去批次的数量。设置为 0 将使用所有先前批次执行累积处理。

StreamingTest 提供流式假设检验。

val data = ssc.textFileStream(dataDir).map(line => line.split(",") match {
  case Array(label, value) => BinarySample(label.toBoolean, value.toDouble)
})

val streamingTest = new StreamingTest()
  .setPeacePeriod(0)
  .setWindowSize(0)
  .setTestMethod("welch")

val out = streamingTest.registerStream(data)
out.print()

在 Spark 仓库的 "examples/src/main/scala/org/apache/spark/examples/mllib/StreamingTestExample.scala" 中查找完整的示例代码。

StreamingTest 提供流式假设检验。

import org.apache.spark.mllib.stat.test.BinarySample;
import org.apache.spark.mllib.stat.test.StreamingTest;
import org.apache.spark.mllib.stat.test.StreamingTestResult;

JavaDStream<BinarySample> data = ssc.textFileStream(dataDir).map(line -> {
  String[] ts = line.split(",");
  boolean label = Boolean.parseBoolean(ts[0]);
  double value = Double.parseDouble(ts[1]);
  return new BinarySample(label, value);
});

StreamingTest streamingTest = new StreamingTest()
  .setPeacePeriod(0)
  .setWindowSize(0)
  .setTestMethod("welch");

JavaDStream<StreamingTestResult> out = streamingTest.registerStream(data);
out.print();

在 Spark 仓库的 "examples/src/main/java/org/apache/spark/examples/mllib/JavaStreamingTestExample.java" 中查找完整的示例代码。

随机数据生成

随机数据生成对于随机算法、原型设计和性能测试非常有用。spark.mllib 支持生成具有从给定分布（均匀分布、标准正态分布或泊松分布）中提取的 i.i.d. 值的随机 RDD。

RandomRDDs 提供了生成随机 double RDD 或 vector RDD 的工厂方法。以下示例生成一个随机 double RDD，其值遵循标准正态分布 N(0, 1)，然后将其映射到 N(1, 4)。

有关 API 的更多详细信息，请参阅 RandomRDDs Python 文档。

from pyspark.mllib.random import RandomRDDs

sc = ... # SparkContext

# Generate a random double RDD that contains 1 million i.i.d. values drawn from the
# standard normal distribution `N(0, 1)`, evenly distributed in 10 partitions.
u = RandomRDDs.normalRDD(sc, 1000000L, 10)
# Apply a transform to get a random double RDD following `N(1, 4)`.
v = u.map(lambda x: 1.0 + 2.0 * x)

有关 API 的详细信息，请参阅 RandomRDDs Scala 文档。

import org.apache.spark.SparkContext
import org.apache.spark.mllib.random.RandomRDDs._

val sc: SparkContext = ...

// Generate a random double RDD that contains 1 million i.i.d. values drawn from the
// standard normal distribution `N(0, 1)`, evenly distributed in 10 partitions.
val u = normalRDD(sc, 1000000L, 10)
// Apply a transform to get a random double RDD following `N(1, 4)`.
val v = u.map(x => 1.0 + 2.0 * x)

有关 API 的详细信息，请参阅 RandomRDDs Java 文档。

import org.apache.spark.SparkContext;
import org.apache.spark.api.JavaDoubleRDD;
import static org.apache.spark.mllib.random.RandomRDDs.*;

JavaSparkContext jsc = ...

// Generate a random double RDD that contains 1 million i.i.d. values drawn from the
// standard normal distribution `N(0, 1)`, evenly distributed in 10 partitions.
JavaDoubleRDD u = normalJavaRDD(jsc, 1000000L, 10);
// Apply a transform to get a random double RDD following `N(1, 4)`.
JavaDoubleRDD v = u.mapToDouble(x -> 1.0 + 2.0 * x);

核密度估计

核密度估计是一种有用的技术，用于可视化经验概率分布，而无需对观察到的样本所来自的特定分布做出假设。它计算随机变量的概率密度函数的估计值，并在给定的点集处进行评估。它通过将特定点处的经验分布的 PDF 表示为以每个样本为中心的正态分布的 PDF 的平均值来实现此估计。

KernelDensity 提供了从样本的 RDD 计算核密度估计的方法。以下示例演示了如何执行此操作。

有关 API 的更多详细信息，请参阅 KernelDensity Python 文档。

from pyspark.mllib.stat import KernelDensity

# an RDD of sample data
data = sc.parallelize([1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 5.0, 5.0, 6.0, 7.0, 8.0, 9.0, 9.0])

# Construct the density estimator with the sample data and a standard deviation for the Gaussian
# kernels
kd = KernelDensity()
kd.setSample(data)
kd.setBandwidth(3.0)

# Find density estimates for the given values
densities = kd.estimate([-1.0, 2.0, 5.0])

在 Spark 仓库的 "examples/src/main/python/mllib/kernel_density_estimation_example.py" 中查找完整的示例代码。

KernelDensity 提供了从样本的 RDD 计算核密度估计的方法。以下示例演示了如何执行此操作。

有关 API 的详细信息，请参阅 KernelDensity Scala 文档。

import org.apache.spark.mllib.stat.KernelDensity
import org.apache.spark.rdd.RDD

// an RDD of sample data
val data: RDD[Double] = sc.parallelize(Seq(1, 1, 1, 2, 3, 4, 5, 5, 6, 7, 8, 9, 9))

// Construct the density estimator with the sample data and a standard deviation
// for the Gaussian kernels
val kd = new KernelDensity()
  .setSample(data)
  .setBandwidth(3.0)

// Find density estimates for the given values
val densities = kd.estimate(Array(-1.0, 2.0, 5.0))

在 Spark 仓库的 "examples/src/main/scala/org/apache/spark/examples/mllib/KernelDensityEstimationExample.scala" 中查找完整的示例代码。

KernelDensity 提供了从样本的 RDD 计算核密度估计的方法。以下示例演示了如何执行此操作。

有关 API 的详细信息，请参阅 KernelDensity Java 文档。

import java.util.Arrays;

import org.apache.spark.api.java.JavaRDD;
import org.apache.spark.mllib.stat.KernelDensity;

// an RDD of sample data
JavaRDD<Double> data = jsc.parallelize(
  Arrays.asList(1.0, 1.0, 1.0, 2.0, 3.0, 4.0, 5.0, 5.0, 6.0, 7.0, 8.0, 9.0, 9.0));

// Construct the density estimator with the sample data
// and a standard deviation for the Gaussian kernels
KernelDensity kd = new KernelDensity().setSample(data).setBandwidth(3.0);

// Find density estimates for the given values
double[] densities = kd.estimate(new double[]{-1.0, 2.0, 5.0});

System.out.println(Arrays.toString(densities));

在 Spark 仓库的 "examples/src/main/java/org/apache/spark/examples/mllib/JavaKernelDensityEstimationExample.java" 中查找完整的示例代码。

MLlib: 主要指南

MLlib：基于 RDD 的 API 指南

基本统计 - 基于 RDD 的 API

概要统计

相关性

分层抽样

假设检验

流式显著性检验

随机数据生成

核密度估计