5

rules パッケージの「検査」機能の結果として、さらに 2 つの対策を追加したいと思います。それらはクルチンスキーとインバランス比です。情報、inspect 関数のコードの場所、およびそれを変更する方法を教えてください。

ありがとう

4

2 に答える 2

3

不均衡は非常に簡単です。

library(arules)
data("Income")
rules <- apriori(Income)

suppA <- support(lhs(rules), trans = Income)
suppB <- support(rhs(rules), trans = Income)
suppAB <- quality(rules)$supp
quality(rules)$imbalance <- abs(suppA - suppB)/(suppA + suppB - suppAB)

inspect(head(rules))
  lhs                                 rhs                               support confidence     lift  imbalance
1 {}                               => {language in home=english}      0.9128854  0.9128854 1.000000 0.03082862
2 {occupation=clerical/service}    => {language in home=english}      0.1127109  0.9292566 1.017933 0.69021050
3 {ethnic classification=hispanic} => {education=no college graduate} 0.1096568  0.8636884 1.224731 0.61395923
4 {dual incomes=no}                => {marital status=married}        0.1400524  0.9441176 2.447871 0.35210356
5 {dual incomes=no}                => {language in home=english}      0.1364165  0.9196078 1.007364 0.63837280
6 {occupation=student}             => {marital status=single}         0.1449971  0.8838652 2.160490 0.34123127

Kulczynski 測度 1/2(P(A|B)+P(B|A)) はもう少しトリッキーです。P(A|B) は、A->B の信頼度です。ただし、P(B|A) の場合、B->A の信頼度が必要です。そのため、左側と右側を入れ替えた新しいルール セットを作成し、信頼度を計算する必要があります。

 confAB <- quality(rules)$conf
 BArules <- new("rules", lhs = rhs(rules), rhs = lhs(rules))
 confBA <- interestMeasure(BArules, method = "confidence", trans = Income)
 quality(rules)$kulczynski <- .5*(confAB + confBA)

 inspect(head(rules))
    lhs                                 rhs                               support confidence     lift  imbalance kulczynski
  1 {}                               => {language in home=english}      0.9128854  0.9128854 1.000000 0.03082862  0.9564427
  2 {occupation=clerical/service}    => {language in home=english}      0.1127109  0.9292566 1.017933 0.69021050  0.5263616
  3 {ethnic classification=hispanic} => {education=no college graduate} 0.1096568  0.8636884 1.224731 0.61395923  0.5095922
  4 {dual incomes=no}                => {marital status=married}        0.1400524  0.9441176 2.447871 0.35210356  0.6536199
  5 {dual incomes=no}                => {language in home=english}      0.1364165  0.9196078 1.007364 0.63837280  0.5345211
  6 {occupation=student}             => {marital status=single}         0.1449971  0.8838652 2.160490 0.34123127  0.6191456
于 2015-09-13T14:16:42.400 に答える