頻出パターンマイニング
頻出アイテム、アイテムセット、サブシーケンス、またはその他のサブ構造をマイニングすることは、通常、大規模データセットを分析する最初のステップであり、データマイニングにおける長年の活発な研究トピックです。詳細については、Wikipedia のアソシエーションルール学習を参照してください。
目次
FP-Growth
FP-growth アルゴリズムは、Han et al., Mining frequent patterns without candidate generationという論文で説明されています。ここで「FP」は頻出パターンを表します。トランザクションのデータセットが与えられたとき、FP-growth の最初のステップは、アイテムの頻度を計算し、頻出アイテムを特定することです。同じ目的のために設計されたApriori のようなアルゴリズムとは異なり、FP-growth の 2 番目のステップでは、通常生成にコストがかかる候補セットを明示的に生成せずにトランザクションをエンコードするためにサフィックスツリー(FP-tree)構造を使用します。2 番目のステップの後、頻出アイテムセットを FP-tree から抽出できます。spark.mllib
では、Li et al., PFP: Parallel FP-growth for query recommendationで説明されているように、PFP という FP-growth の並列バージョンを実装しました。PFP はトランザクションのサフィックスに基づいて FP-tree を成長させる作業を分散するため、シングルマシンの実装よりもスケーラブルです。詳細については、論文を参照してください。
FP-growth は *アイテムセット* で動作します。アイテムセットは、一意のアイテムの順序付けられていないコレクションです。Spark には *set* 型がないため、アイテムセットは配列として表現されます。
spark.ml
の FP-growth 実装は、以下の(ハイパー)パラメーターを受け取ります
minSupport
:アイテムセットが頻出と識別されるための最小サポート。たとえば、アイテムが 5 つのトランザクションのうち 3 つに現れる場合、そのサポートは 3/5 = 0.6 です。minConfidence
:アソシエーションルールを生成するための最小確信度。確信度は、アソシエーションルールが真であることが判明した頻度を示す指標です。たとえば、トランザクションでアイテムセットX
が 4 回出現し、X
とY
が共起するのが 2 回のみである場合、ルールX => Y
の確信度は 2/4 = 0.5 になります。このパラメーターは、頻出アイテムセットのマイニングには影響しませんが、頻出アイテムセットからアソシエーションルールを生成するための最小確信度を指定します。numPartitions
:作業を分散するために使用されるパーティションの数。デフォルトでは、パラメーターは設定されておらず、入力データセットのパーティション数が使用されます。
FPGrowthModel
は以下を提供します。
freqItemsets
:次の列を持つ DataFrame 形式の頻出アイテムセットitems: array
: 指定されたアイテムセット。freq: long
: このアイテムセットがモデルの構成パラメーターで観測された回数。
associationRules
:minConfidence
を超える確信度で生成されたアソシエーションルール。次の列を持つ DataFrame 形式です。antecedent: array
: アソシエーションルールの仮説であるアイテムセット。consequent: array
: アソシエーションルールの結論を表す単一の要素を常に含むアイテムセット。confidence: double
:confidence
の定義については、上記のminConfidence
を参照してください。lift: double
: 先行詞が結果をどれだけうまく予測するかを測る指標。support(antecedent U consequent) / (support(antecedent) x support(consequent))
として計算されます。support: double
:support
の定義については、上記のminSupport
を参照してください。
transform
:itemsCol
の各トランザクションについて、transform
メソッドは、そのアイテムを各アソシエーションルールの先行詞と比較します。レコードに特定のアソシエーションルールのすべての先行詞が含まれている場合、そのルールが適用可能であると見なされ、その結果が予測結果に追加されます。transform メソッドは、適用可能なすべてのルールからの結果を予測として要約します。予測列はitemsCol
と同じデータ型を持ち、itemsCol
に既存のアイテムは含まれません。
例
詳細については、Python API ドキュメントを参照してください。
from pyspark.ml.fpm import FPGrowth
df = spark.createDataFrame([
(0, [1, 2, 5]),
(1, [1, 2, 3, 5]),
(2, [1, 2])
], ["id", "items"])
fpGrowth = FPGrowth(itemsCol="items", minSupport=0.5, minConfidence=0.6)
model = fpGrowth.fit(df)
# Display frequent itemsets.
model.freqItemsets.show()
# Display generated association rules.
model.associationRules.show()
# transform examines the input items against all the association rules and summarize the
# consequents as prediction
model.transform(df).show()
詳細については、Scala API ドキュメントを参照してください。
import org.apache.spark.ml.fpm.FPGrowth
val dataset = spark.createDataset(Seq(
"1 2 5",
"1 2 3 5",
"1 2")
).map(t => t.split(" ")).toDF("items")
val fpgrowth = new FPGrowth().setItemsCol("items").setMinSupport(0.5).setMinConfidence(0.6)
val model = fpgrowth.fit(dataset)
// Display frequent itemsets.
model.freqItemsets.show()
// Display generated association rules.
model.associationRules.show()
// transform examines the input items against all the association rules and summarize the
// consequents as prediction
model.transform(dataset).show()
詳細については、Java API ドキュメントを参照してください。
import java.util.Arrays;
import java.util.List;
import org.apache.spark.ml.fpm.FPGrowth;
import org.apache.spark.ml.fpm.FPGrowthModel;
import org.apache.spark.sql.Dataset;
import org.apache.spark.sql.Row;
import org.apache.spark.sql.RowFactory;
import org.apache.spark.sql.SparkSession;
import org.apache.spark.sql.types.*;
List<Row> data = Arrays.asList(
RowFactory.create(Arrays.asList("1 2 5".split(" "))),
RowFactory.create(Arrays.asList("1 2 3 5".split(" "))),
RowFactory.create(Arrays.asList("1 2".split(" ")))
);
StructType schema = new StructType(new StructField[]{ new StructField(
"items", new ArrayType(DataTypes.StringType, true), false, Metadata.empty())
});
Dataset<Row> itemsDF = spark.createDataFrame(data, schema);
FPGrowthModel model = new FPGrowth()
.setItemsCol("items")
.setMinSupport(0.5)
.setMinConfidence(0.6)
.fit(itemsDF);
// Display frequent itemsets.
model.freqItemsets().show();
// Display generated association rules.
model.associationRules().show();
// transform examines the input items against all the association rules and summarize the
// consequents as prediction
model.transform(itemsDF).show();
詳細については、R API ドキュメントを参照してください。
# Load training data
df <- selectExpr(createDataFrame(data.frame(rawItems = c(
"1,2,5", "1,2,3,5", "1,2"
))), "split(rawItems, ',') AS items")
fpm <- spark.fpGrowth(df, itemsCol="items", minSupport=0.5, minConfidence=0.6)
# Extracting frequent itemsets
spark.freqItemsets(fpm)
# Extracting association rules
spark.associationRules(fpm)
# Predict uses association rules to and combines possible consequents
predict(fpm, df)
PrefixSpan
PrefixSpan は、Pei et al., Mining Sequential Patterns by Pattern-Growth: The PrefixSpan Approachで説明されているシーケンシャルパターンマイニングアルゴリズムです。シーケンシャルパターンマイニング問題の形式化については、参照論文を参照してください。
spark.ml
の PrefixSpan 実装は、次のパラメーターを受け取ります。
minSupport
:頻出シーケンシャルパターンと見なされるために必要な最小サポート。maxPatternLength
:頻出シーケンシャルパターンの最大長。この長さを超える頻出パターンは、結果に含まれません。maxLocalProjDBSize
:射影データベースのローカル反復処理が開始される前に、プレフィックス射影データベースで許可されるアイテムの最大数。このパラメーターは、エグゼキューターのサイズに関して調整する必要があります。sequenceCol
:データセット内のシーケンス列の名前(デフォルトは「sequence」)。この列に null がある行は無視されます。
例
詳細については、Python API ドキュメントを参照してください。
from pyspark.ml.fpm import PrefixSpan
df = sc.parallelize([Row(sequence=[[1, 2], [3]]),
Row(sequence=[[1], [3, 2], [1, 2]]),
Row(sequence=[[1, 2], [5]]),
Row(sequence=[[6]])]).toDF()
prefixSpan = PrefixSpan(minSupport=0.5, maxPatternLength=5,
maxLocalProjDBSize=32000000)
# Find frequent sequential patterns.
prefixSpan.findFrequentSequentialPatterns(df).show()
詳細については、Scala API ドキュメントを参照してください。
import org.apache.spark.ml.fpm.PrefixSpan
val smallTestData = Seq(
Seq(Seq(1, 2), Seq(3)),
Seq(Seq(1), Seq(3, 2), Seq(1, 2)),
Seq(Seq(1, 2), Seq(5)),
Seq(Seq(6)))
val df = smallTestData.toDF("sequence")
val result = new PrefixSpan()
.setMinSupport(0.5)
.setMaxPatternLength(5)
.setMaxLocalProjDBSize(32000000)
.findFrequentSequentialPatterns(df)
.show()
詳細については、Java API ドキュメントを参照してください。
import java.util.Arrays;
import java.util.List;
import org.apache.spark.ml.fpm.PrefixSpan;
import org.apache.spark.sql.Dataset;
import org.apache.spark.sql.Row;
import org.apache.spark.sql.RowFactory;
import org.apache.spark.sql.SparkSession;
import org.apache.spark.sql.types.*;
List<Row> data = Arrays.asList(
RowFactory.create(Arrays.asList(Arrays.asList(1, 2), Arrays.asList(3))),
RowFactory.create(Arrays.asList(Arrays.asList(1), Arrays.asList(3, 2), Arrays.asList(1,2))),
RowFactory.create(Arrays.asList(Arrays.asList(1, 2), Arrays.asList(5))),
RowFactory.create(Arrays.asList(Arrays.asList(6)))
);
StructType schema = new StructType(new StructField[]{ new StructField(
"sequence", new ArrayType(new ArrayType(DataTypes.IntegerType, true), true),
false, Metadata.empty())
});
Dataset<Row> sequenceDF = spark.createDataFrame(data, schema);
PrefixSpan prefixSpan = new PrefixSpan().setMinSupport(0.5).setMaxPatternLength(5);
// Finding frequent sequential patterns
prefixSpan.findFrequentSequentialPatterns(sequenceDF).show();
詳細については、R API ドキュメントを参照してください。
# Load training data
df <- createDataFrame(list(list(list(list(1L, 2L), list(3L))),
list(list(list(1L), list(3L, 2L), list(1L, 2L))),
list(list(list(1L, 2L), list(5L))),
list(list(list(6L)))),
schema = c("sequence"))
# Finding frequent sequential patterns
frequency <- spark.findFrequentSequentialPatterns(df, minSupport = 0.5, maxPatternLength = 5L,
maxLocalProjDBSize = 32000000L)
showDF(frequency)