Deprecated: The each() function is deprecated. This message will be suppressed on further calls in /home/zhenxiangba/zhenxiangba.com/public_html/phproxy-improved-master/index.php on line 456
JP4839358B2 - Purchase prediction method, purchase prediction apparatus, purchase prediction program, and recording medium - Google Patents
[go: Go Back, main page]

JP4839358B2 - Purchase prediction method, purchase prediction apparatus, purchase prediction program, and recording medium - Google Patents

Purchase prediction method, purchase prediction apparatus, purchase prediction program, and recording medium Download PDF

Info

Publication number
JP4839358B2
JP4839358B2 JP2008227235A JP2008227235A JP4839358B2 JP 4839358 B2 JP4839358 B2 JP 4839358B2 JP 2008227235 A JP2008227235 A JP 2008227235A JP 2008227235 A JP2008227235 A JP 2008227235A JP 4839358 B2 JP4839358 B2 JP 4839358B2
Authority
JP
Japan
Prior art keywords
probability
purchase
class
user
degree
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
JP2008227235A
Other languages
Japanese (ja)
Other versions
JP2010061456A (en
Inventor
具治 岩田
晋治 渡部
武士 山田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NTT Inc
NTT Inc USA
Original Assignee
Nippon Telegraph and Telephone Corp
NTT Inc USA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nippon Telegraph and Telephone Corp, NTT Inc USA filed Critical Nippon Telegraph and Telephone Corp
Priority to JP2008227235A priority Critical patent/JP4839358B2/en
Publication of JP2010061456A publication Critical patent/JP2010061456A/en
Application granted granted Critical
Publication of JP4839358B2 publication Critical patent/JP4839358B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Description

本発明は、複数のユーザの複数の商品に対する時系列の購買データに基づいて、その後の購買状況を予測する技術に関する。   The present invention relates to a technique for predicting a subsequent purchase situation based on time-series purchase data for a plurality of products of a plurality of users.

近年、オンライン技術やデータベース技術の発展により、膨大な購買データを日々収集できるようになった。購買データに基づいてその後の購買状況を予測することで、ユーザにお勧めの商品を教示したり、購買傾向を解析したりすることが可能となる。   In recent years, with the development of online technology and database technology, a huge amount of purchase data can be collected every day. By predicting the subsequent purchase status based on the purchase data, it is possible to teach the recommended product to the user and analyze the purchase tendency.

従来の購買状況の予測法として、購買データから頻出する購買順序パターンを抽出し、そのパターンを基に予測する手法がある(非特許文献1参照)。この手法では、購買の時間情報(何月何日に購入したか等)が考慮されない。また、購買の時間情報を考慮して購買状況を予測する手法もある(非特許文献2)。
岩田具治、山田武士、上田修功、“購買順序を効率的に用いた協調フィルタリング”、情報処理学会論文誌:数理モデル化と応用、Vol.49、No.SIG4(TOM20)、p.125−134、2008 岩田具治、田中利幸、山田武士、上田修功、“分布が変化するデータにおけるモデル学習法”、情報論的学習理論ワークショップ(IBIS2007)、p.170−175、2007
As a conventional method for predicting a purchase situation, there is a method of extracting a purchase order pattern frequently appearing from purchase data and predicting based on the pattern (see Non-Patent Document 1). This method does not take into account purchase time information (month, date, etc.). There is also a method for predicting a purchase situation in consideration of purchase time information (Non-patent Document 2).
Tomoharu Iwata, Takeshi Yamada, Nobuyoshi Ueda, “Collaborative Filtering Using Purchase Order Efficiently”, IPSJ Transactions: Mathematical Modeling and Applications, Vol. 49, no. SIG4 (TOM20), p. 125-134, 2008 Tomoharu Iwata, Toshiyuki Tanaka, Takeshi Yamada, Nobuyoshi Ueda, “Model learning method for data with varying distributions”, Information theory learning theory workshop (IBIS2007), p. 170-175, 2007

しかし、非特許文献1の技術では、購買の時間情報が考慮されないため、流行の変化や興味の変化などに対応できないという問題がある。また、非特許文献2の技術では、データが追加されるたびに、全データを読み込み、計算を行う必要があるため、日々膨大なデータが蓄積される場合、逐次対応することが困難という問題がある。   However, the technique of Non-Patent Document 1 has a problem in that it cannot cope with a change in fashion or a change in interest because purchase time information is not taken into consideration. Further, in the technique of Non-Patent Document 2, since it is necessary to read all data and perform calculation every time data is added, there is a problem that it is difficult to sequentially cope with a huge amount of data accumulated every day. is there.

そこで、本発明は、前記問題に鑑みてなされたものであり、購買の時間情報を考慮して購買データの経時的な変化に対応しながらも少ない計算量で購買状況を予測することを課題とする。   Therefore, the present invention has been made in view of the above problems, and it is an object to predict a purchase situation with a small amount of calculation while dealing with changes in purchase data over time in consideration of purchase time information. To do.

前記課題を解決するために、本発明は、複数のユーザの複数の商品に対する購買のデータが時系列に沿って予め集められた購買データに基づいて、任意のユーザの任意の商品に対するその後の購入の可能性を予測する場合に、複数のユーザを複数のクラスに分け、同じクラスに属するユーザは同じ確率で任意の商品を購入すると仮定する確率モデルを用いて購入の可能性を予測する購入予測装置であって、購買データ、ユーザごとのクラスそれぞれに属する確率の集合であるクラス確率集合、ユーザのクラスに属する確率が変わらない度合いを示すクラス確率不変化度をパラメータとして持つ、ユーザのクラス確率の事前確率、クラスごとの商品それぞれがそのクラスのユーザによって購入される確率の集合である購入確率集合、および、クラスに属するユーザによって商品を購入される確率が変わらない度合いを示す購入確率不変化度をパラメータとして持つ、クラスの購入確率の事前確率、を記憶する記憶部と、クラス確率不変化度の尤度と購入確率不変化度の尤度とが最大になるように、クラス確率不変化度と購入確率不変化度とを繰り返し計算し、その計算したクラス確率不変化度と購入確率不変化度とを用いて、クラス確率集合および購入確率集合を更新するモデル推定部と、更新したクラス確率集合および購入確率集合を用いて、ユーザが商品を購入する確率を算出する購入予測部と、を備えることを特徴とする。また、この購入予測装置を用いて購入予測方法を実行することもできる。 In order to solve the above-described problem, the present invention is based on purchase data in which purchase data for a plurality of products of a plurality of users is pre-collected in time series, and subsequent purchases for any product of any user. When predicting the probability of purchase, purchase prediction is performed using a probability model that divides multiple users into multiple classes and assumes that users belonging to the same class will purchase any product with the same probability. The class probability of a user, which is a device, having as parameters a purchase probability, a class probability set that is a set of probabilities belonging to each class for each user, and a class probability invariance indicating the degree to which the probability that the user belongs to a class does not change prior probability, each item of each class purchase probability set is a set of probability to be purchased by the user of the class, and, click With purchase probability not change degree indicating a degree that does not change probability of purchase items by a user belonging to the scan as a parameter, a storage unit for storing the prior probability, the purchase probability of class, the class probability immutabilization degree likelihood And the probability of purchase probability unchanged so that the likelihood of purchase probability unchanged is maximized, and the class probability unchanged and purchase probability unchanged are repeatedly calculated, and the calculated class probability unchanged and purchase probability unchanged are calculated. used, and the model estimation unit you update the class probability set and purchase probability set, using the updated class probability set and purchase probability set, and the purchase prediction unit that to calculate the probability that a user purchases the goods It is characterized by providing. Moreover, a purchase prediction method can also be executed using this purchase prediction apparatus.

かかる発明によれば、複数のユーザが複数のクラスに分けられるという前提でモデル推定部がクラス確率集合および購入確率集合を推定(更新)し、購入予測部がそのクラス確率集合および購入確率集合に基づいて購買状況(購入の可能性)を予測することで、時系列の購買データの経時的な変化に対応しながらも少ない計算量で購買状況を予測することができる。 According to the invention, the model estimator assuming that a plurality of users are divided into a plurality of classes class probabilities set and purchase probability set estimated (updated), purchased prediction unit within that class probabilities set and purchase probability set By predicting the purchase situation (probability of purchase) based on this, it is possible to predict the purchase situation with a small amount of calculation while responding to changes over time in purchase data in time series.

また、本発明は、記憶部が、購買データの最新時刻以前の所定時間分のクラス確率集合、ユーザのクラス確率の事前確率、購入確率集合、および、クラスの購入確率の事前確率、を記憶しており、モデル推定部が、所定時間分のクラス確率不変化度の尤度と、所定時間分の購入確率不変化度の尤度と、が最大になるように、所定時間分のクラス確率不変化度と、所定時間分の購入確率不変化度と、を繰り返し計算し、その計算した所定時間分のクラス確率不変化度と、所定時間分の購入確率不変化度と、を用いて、所定時間分のクラス確率集合、および、所定時間分の購入確率集合を更新し、購入予測部が、更新した所定時間分のクラス確率集合、および、所定時間分の購入確率集合を用いて、ユーザが商品を購入する確率を算出することが好ましい。 Further, in the present invention, the storage unit stores a class probability set , a user class probability prior probability, a purchase probability set , and a class purchase probability prior probability for a predetermined time before the latest time of purchase data. The class estimation unit has a class probability for a predetermined time so that the likelihood of the class probability unchanged for a predetermined time and the likelihood of the purchase probability unchanged for a predetermined time are maximized. The degree of change and the purchase probability unchanged degree for a predetermined time are repeatedly calculated, and the calculated class probability unchanged degree for the predetermined time and the purchase probability unchanged degree for the predetermined time are used. The class probability set for a predetermined time and the purchase probability set for a predetermined time are updated, and the purchase predicting unit uses the updated class probability set for the predetermined time and the purchase probability set for the predetermined time. There possible to calculate the probability to purchase goods Preferred.

かかる発明によれば、クラス確率集合および購入確率集合が、購買データの最新時刻分だけでなく、それ以前の所定時間分を考慮したものであることで、購買データの最新時刻分だけのときに比べて予測の精度を高くすることができる。   According to this invention, when the class probability set and the purchase probability set are not only for the latest time of the purchase data but also for a predetermined time before that, only when the latest time of the purchase data is Compared with this, the prediction accuracy can be increased.

また、本発明は、コンピュータを、購入予測装置の各部として機能させるための購入予測プログラムである。これにより、このプログラムをインストールされたコンピュータは、このプログラムに基づいた各機能を実現することができる。   Moreover, this invention is a purchase prediction program for functioning a computer as each part of a purchase prediction apparatus. Thereby, the computer installed with this program can realize each function based on this program.

また、本発明は、前記購入予測プログラムが記録されたことを特徴とするコンピュータに読み取り可能な記録媒体である。これにより、この記録媒体を装着されたコンピュータは、この記録媒体に記録されたプログラムに基づいた各機能を実現することができる。   The present invention is a computer-readable recording medium in which the purchase prediction program is recorded. Thereby, the computer equipped with this recording medium can realize each function based on the program recorded on this recording medium.

本発明によれば、購買の時間情報を考慮して購買データの経時的な変化に対応しながらも少ない計算量で購買状況を予測することができる。   According to the present invention, it is possible to predict a purchase situation with a small amount of calculation while dealing with changes in purchase data over time in consideration of purchase time information.

以下、本発明を実施するための最良の形態(以下、「実施形態」という。)について、詳細に説明する。なお、本実施形態では、複数のユーザの複数の商品に対する購買のデータが時系列に沿って予め集められた購買データに基づいて、任意のユーザの任意の商品に対するその後の購入の可能性を予測する場合に、複数のユーザを複数のクラスに分け、属するクラスが同じユーザは同じ確率で任意の商品を購入すると仮定する確率モデルを用いることを前提とする。   Hereinafter, the best mode for carrying out the present invention (hereinafter referred to as “embodiment”) will be described in detail. Note that in this embodiment, the possibility of subsequent purchases for any product by any user is predicted based on purchase data in which purchase data for a plurality of products by a plurality of users is collected in time series. In this case, it is assumed that a plurality of users are divided into a plurality of classes, and a probability model that assumes that users belonging to the same class purchase arbitrary products with the same probability is used.

また、以下において、「購入の可能性の予測」、「購入する確率の推定(算出)」、「購買状況の予測」、「興味の予測」などの表現は同義である。さらに、以下の各数式において、太字の文字(θなど)は、複数の成分を有していることを示す。また、それらの太字の文字について、文章中では太字で示していないが、各数式と整合をとったものであるものとする。   In the following, expressions such as “prediction of purchase possibility”, “estimation (calculation) of purchase probability”, “prediction of purchase status”, and “prediction of interest” are synonymous. Furthermore, in each of the following mathematical formulas, a bold character (such as θ) indicates that it has a plurality of components. Also, these bold characters are not shown in bold in the text, but are assumed to be consistent with each mathematical expression.

図1は、本実施形態の興味予測装置の構成を示すブロック図である。図1に示すように、興味予測装置1(購入予測装置)は、演算手段2と、入力手段3と、記憶手段4と、出力手段5とを備えており、それらはバスライン11に接続されている。   FIG. 1 is a block diagram showing the configuration of the interest prediction apparatus of this embodiment. As shown in FIG. 1, the interest prediction device 1 (purchase prediction device) includes a calculation means 2, an input means 3, a storage means 4, and an output means 5, which are connected to a bus line 11. ing.

演算手段2は、モデル推定部21と、興味予測部22(購入予測部)と、メモリ23とを含んで構成される。演算手段2は、例えば、CPU(Central Processing Unit)およびRAM(Random Access Memory)から構成され、CPUが記憶手段4からモデル推定プログラム41と興味予測プログラム42とを読み込み、RAMであるメモリ23に格納し、実行することで、モデル推定部21と興味予測部22とがそれぞれ実現される。   The calculation means 2 includes a model estimation unit 21, an interest prediction unit 22 (purchase prediction unit), and a memory 23. The calculation means 2 is composed of, for example, a CPU (Central Processing Unit) and a RAM (Random Access Memory), and the CPU reads the model estimation program 41 and the interest prediction program 42 from the storage means 4 and stores them in the memory 23 which is a RAM. And by executing, the model estimation part 21 and the interest prediction part 22 are each implement | achieved.

入力手段3は、キーボードやディスクドライブ装置などから構成される。入力データは入力手段3を介して興味予測装置1の管理者などによって入力され、記憶手段4に記憶される。   The input means 3 includes a keyboard and a disk drive device. The input data is input by the administrator of the interest prediction device 1 via the input unit 3 and stored in the storage unit 4.

記憶手段4は、ハードディスク装置などから構成され、プログラム格納部40aと、データ格納部40bと、を有する。プログラム格納部40aは、モデル推定プログラム41と、興味予測プログラム42と、を記憶する。また、データ格納部40bは、入力データ43、クラス集合44、クラス確率集合45、購入確率集合46、クラス確率不変化度集合47、購入確率不変化度集合48を記憶する。   The storage unit 4 includes a hard disk device and the like, and includes a program storage unit 40a and a data storage unit 40b. The program storage unit 40a stores a model estimation program 41 and an interest prediction program 42. The data storage unit 40b stores input data 43, a class set 44, a class probability set 45, a purchase probability set 46, a class probability invariance set 47, and a purchase probability invariance set 48.

入力データ43は、複数のユーザの複数の商品に対する購買のデータが時系列に沿って予め集められた購買データである。クラス集合44は、入力データ43における複数のユーザを分ける複数のクラスの集合データである。   The input data 43 is purchase data in which purchase data for a plurality of products of a plurality of users is collected in advance in time series. The class set 44 is set data of a plurality of classes that divide a plurality of users in the input data 43.

クラス確率集合45は、ユーザごとのクラスそれぞれに属する確率の集合である。購入確率集合46は、クラスごとの、商品それぞれがそのクラスのユーザによって購入される確率の集合である。   The class probability set 45 is a set of probabilities belonging to each class for each user. The purchase probability set 46 is a set of probabilities that each product is purchased by a user of the class for each class.

クラス確率不変化度集合47は、ユーザごとのクラスそれぞれに属する確率が変わらない度合いの集合である。購入確率不変化度集合48は、クラスごとの、商品それぞれがそのクラスのユーザによって購入される確率が変わらない度合いの集合である。   The class probability invariance degree set 47 is a set to which the probability of belonging to each class for each user does not change. The purchase probability invariance degree set 48 is a set for each class to the extent that the probability that each product is purchased by a user of the class does not change.

出力手段5は、例えば、グラフィックボード及びそれに接続されたモニタであり、興味予測を行った結果などを表示するものである。   The output means 5 is, for example, a graphic board and a monitor connected to the graphic board, and displays the result of interest prediction.

本実施形態では、ユーザuが時刻tにおいて商品iを購入する確率P(i|u,t)を推定(算出)する。つまり、確率P(i|u,t)は、ユーザuの時刻tにおける商品iに対する興味の度合いを表す。時刻tは離散値とし、単位時間として1日、1週間など任意の期間を用いることができる。   In the present embodiment, the probability P (i | u, t) that the user u purchases the product i at time t is estimated (calculated). That is, the probability P (i | u, t) represents the degree of interest in the product i at the time t of the user u. The time t is a discrete value, and an arbitrary period such as one day or one week can be used as a unit time.

確率P(i|u,t)を推定する場合、ユーザの合計数をU、商品の合計数をIとしたとき、各時刻において、U(I−1)個のパラメータが必要となる。なお、「U(I−1)」において、「I」ではなく「I−1」となっているのは、「I−1」個のパラメータが決まれば、残りの1つのパラメータは、パラメータの全体合計値から「I−1」個のパラメータの合計値を減算することで自動的に決定するからである。通常、U(I−1)はUやIに比べて非常に多く、推定が困難である。そこで、隠れ変数であるクラスzを導入して、式(1)のように確率P(i|u,t)をモデル化する。   When estimating the probability P (i | u, t), if the total number of users is U and the total number of products is I, U (I-1) parameters are required at each time. In “U (I-1)”, “I-1” instead of “I” means that if “I-1” parameters are determined, the remaining one parameter is the parameter This is because it is automatically determined by subtracting the total value of “I−1” parameters from the total total value. Usually, U (I-1) is much larger than U and I and is difficult to estimate. Therefore, a class z that is a hidden variable is introduced to model the probability P (i | u, t) as shown in Equation (1).

Figure 0004839358
Figure 0004839358

Zをクラスの合計数とすると、このとき各時刻のパラメータの合計数はU(Z−1)+Z(I−1)となり、Z≪U、Z≪Iであるため、パラメータの合計数を大幅に削減でき、頑健な推定が可能となる。ここで、P(,t)は時刻tでのユーザuがクラスzに帰属する確率、P(i|z,t)は時刻tでのクラスzで商品iが購入される確率を表す。 If Z is the total number of classes, the total number of parameters at each time is U (Z-1) + Z (I-1). Since Z << U and Z << I, the total number of parameters is greatly increased. Can be reduced to a robust estimate. Here, P ( z | u , t) is a probability that the user u belongs to the class z at the time t, and P (i | z, t) is a probability that the product i is purchased at the class z at the time t. To express.

商品に対する興味の変化、人気の変化に対応するため、P(z|u,t)、および、P(i|z,t)は、以下の式(2)、式(3)のように、1時刻前の状態を引数とするある関数Fで時間発展するものとする。   P (z | u, t) and P (i | z, t) are expressed by the following formulas (2) and (3) in order to cope with changes in interest and popularity of products. It is assumed that the time evolves with a function F that takes the state one time ago as an argument.

Figure 0004839358
Figure 0004839358

時間発展として、例えば、状態空間モデルを用いて、カルマンフィルタで推定することができる。また、φtuz=P(z|t,u)(=P(z|u,t))、θtzi=P(i|z,t)とし、以下のような時間発展も考えられる。 For example, the time evolution can be estimated by a Kalman filter using a state space model. Further, φ tuz = P (z | t, u) (= P (z | u, t)) and θ tzi = P (i | z, t) are assumed, and the following time development is also conceivable.

ここで、時刻tのユーザuのクラス確率φtuの事前確率は、1時刻前のユーザuのクラス確率φt−1u(「t−1u」は「(t−1)u」の意。以下、他の変数についても同様)およびユーザuのクラス確率不変化度αtuをパラメータとして持つ、以下の式(4)のようなディリクレ分布であるとする。 Here, the prior probability of the class u φ tu of the user u at time t is the class probability φ t−1u of the user u one time before (“t−1u” means “(t−1) u”. The same applies to other variables) and the Dirichlet distribution as in the following expression (4) having the class u invariance α tu of the user u as a parameter.

Figure 0004839358
Figure 0004839358

なお、この分布の平均値Eは1時刻前のクラス確率に一致し
(E(φtu|φt−1u,αtu)=φt−1u)、
分散varはαtu+1に反比例する(var(φtuz|φt−1u,αtu)=
(φt−1uz(1−φt−1uz))/αtu+1)。
The average value E of this distribution coincides with the class probability one time ago (E (φ tu | φ t−1u , α tu ) = φ t−1u ),
The variance var is inversely proportional to α tu +1 (var (φ tuz | φ t−1u , α tu ) =
([ Phi] t-1uz (1- [phi] t-1uz )) / [alpha] tu + 1).

また、時刻tのクラスzの購入確率θtzの事前確率も同様に、1時刻前のクラスzの購入確率θt−1zおよびクラスzの購入確率不変化度βtzをパラメータとして持つ、以下の式(5)のようなディリクレ分布であるとする。なお、Γ()はガンマ関数である。 Similarly, the prior probability of the purchase probability θ tz of the class z at the time t similarly has the purchase probability θ t−1z of the class z one time ago and the purchase probability unchanged degree β tz of the class z as parameters as follows: It is assumed that the Dirichlet distribution is as in Expression (5). Note that Γ () is a gamma function.

Figure 0004839358
Figure 0004839358

<モデル推定部>
図2、図3を参照しながら、モデル推定部21について説明する。図2は、本実施形態の興味予測装置のモデル推定部のブロック図を含む説明図である。図2に示すように、モデル推定部21は、入力データ読み込み部211と、パラメータ読み込み部212と、初期化部213と、クラス推定部214と、クラス確率不変化度推定部215と、購入確率不変化度推定部216と、パラメータ書き込み部217と、を備えている。
<Model estimation unit>
The model estimation unit 21 will be described with reference to FIGS. 2 and 3. FIG. 2 is an explanatory diagram including a block diagram of a model estimation unit of the interest prediction apparatus of the present embodiment. As shown in FIG. 2, the model estimation unit 21 includes an input data reading unit 211, a parameter reading unit 212, an initialization unit 213, a class estimation unit 214, a class probability invariance estimation unit 215, and a purchase probability. A non-change degree estimation unit 216 and a parameter writing unit 217 are provided.

図3は、モデル推定部の処理の流れを示すフローチャートである。まず、ステップS1において、入力データ読み込み部211により入力データ43を読み込み、パラメータ読み込み部212により1時刻前のデータを用いて推定した(ステップS1〜S7の処理を行った)クラス確率集合45、購入確率集合46を読み込む。なお、1時刻前のクラス確率集合45および購入確率集合46がない場合、つまりt=0の場合、クラス確率および購入確率は一様分布とする。   FIG. 3 is a flowchart showing a process flow of the model estimation unit. First, in step S1, the input data reading unit 211 reads the input data 43, and the parameter reading unit 212 estimates using the data one hour before (the processing in steps S1 to S7 is performed). Read probability set 46. When there is no class probability set 45 and purchase probability set 46 one time ago, that is, when t = 0, the class probability and purchase probability are uniformly distributed.

次に、ステップS2において、初期化部213により、クラス集合44、クラス確率不変化度集合47、購入確率不変化度集合48を初期化する。例えば、クラス集合44はランダムな値となるように初期化し、クラス確率不変化度集合47および購入確率不変化度集合48は全て「1」に初期化する。   Next, in step S2, the initialization unit 213 initializes the class set 44, the class probability unchanged degree set 47, and the purchase probability unchanged degree set 48. For example, the class set 44 is initialized to be a random value, and the class probability invariance set 47 and the purchase probability invariance set 48 are all initialized to “1”.

次に、クラス推定部214によるクラス集合44の計算(ステップS3)、クラス確率不変化度推定部215によるクラス確率不変化度集合47の計算(ステップS4)、購入確率不変化度推定部216による購入確率不変化度集合48の計算(ステップS5)を繰り返す(ステップS6(終了条件を満たすか否かの判断)でNo→ステップS3)ことにより、クラス集合44、クラス確率不変化度集合47および購入確率不変化度集合48の現在の推定値が、どのくらい入力データ43をもっともらしく説明できているかを表す尤度を最大化する。例えば、商品に対して同じ興味を持つユーザの購買(購買データ。以下同様)が同じクラスに帰属していたり、共通ユーザ(クラスが共通であるユーザ)に購入されやすい商品の購買が同じクラスに帰属していたりする場合、尤度は高くなる。尤度計算は、入力データ43、1時刻前のクラス確率集合45および購入確率集合46を用いることで実現できる。例えば、次の式(6)を尤度として用いることができる。   Next, calculation of the class set 44 by the class estimation unit 214 (step S3), calculation of the class probability invariance set 47 by the class probability invariance estimation unit 215 (step S4), and by the purchase probability invariance estimation unit 216 By repeating the calculation (step S5) of the purchase probability invariance set 48 (step S6 (determining whether or not the end condition is satisfied) → step S3), the class set 44, the class probability invariance set 47 and The likelihood that the current estimated value of the purchase probability unchanged set 48 can explain the input data 43 is maximized. For example, purchases of users who have the same interest in products (purchase data; the same applies hereinafter) belong to the same class, or purchases of products that are easily purchased by common users (users with the same class) belong to the same class. If it belongs, the likelihood is high. The likelihood calculation can be realized by using the input data 43, the class probability set 45 and the purchase probability set 46 one time before. For example, the following equation (6) can be used as the likelihood.

Figure 0004839358
Figure 0004839358

ここで、xtujは時刻tでユーザuがj番目に購入した商品であり、xtu={xtujは時刻tでユーザuが購入した商品集合であるとすると、X={xtuは時刻tでの全ユーザの購入商品集合である。また、φtuz=P(z|t,u)は時刻tでのユーザuがクラスzに帰属する確率(式(4)参照)、φtu={φtuzは時刻tでのユーザuのクラス確率、Φ={φtuは時刻tでのクラス確率集合、θtzi=P(i|z,t)は時刻tでのクラスzで商品iが購入される確率(式(5)参照)、θtz={θtziは時刻tでのクラスzの購入確率、Θ={θtzは時刻tでの購入確率集合とする。また、α={αはクラス確率不変化度集合、β={βは購入確率不変化度集合を表す。 Here, x tuj is a product purchased by user u at time t at j time, and x tu = {x tuj } j is a set of products purchased by user u at time t, X t = {x tu } u is a set of purchased products of all users at time t. Also, φ tuz = P (z | t, u) is the probability that user u at time t belongs to class z (see equation (4)), and φ tu = {φ tuz } z is user u at time t. Φ t = {φ tu } u is a set of class probabilities at time t, and θ tzi = P (i | z, t) is a probability that a product i is purchased with class z at time t (formula ( 5), θ tz = {θ tzi } i is a purchase probability of class z at time t, and Θ t = {θ tz } z is a purchase probability set at time t. Also, α t = {α u } u represents a class probability invariance set, and β t = {β z } z represents a purchase probability invariance set.

ステップS3において、クラス推定部214では、各購買mがどのクラスに帰属するかを計算する。クラス推定部214によって、クラス計算は、入力として、購買mを除いた時刻tでユーザuの購買がクラスkである数ntuk\m、購買mを除いた時刻tでユーザuの購買数ntu\m、購買mを除いた時刻tで商品iの購買がクラスkである数ntik\m、購買mを除いた時刻tで購買がクラスkである数ntk\mをとり、例えば、ギブスサンプリングに基づく次の式(7)により計算される。ここで、Zt\m={ztnn≠mは購買mを除いた時刻tでのクラス集合を表す。 In step S3, the class estimation unit 214 calculates which class each purchase m belongs to. The class calculation by the class estimator 214 is as follows: the number n tuk \ m that the purchase of the user u is class k at the time t excluding the purchase m, and the number of purchases n of the user u at the time t excluding the purchase m. take tu\m, number n tik\m purchasing of goods i at time t, except for the purchase m is a class k, the number n tk\m purchase is a class k at time t, except for the purchase m, for example The following equation (7) based on Gibbs sampling is used. Here, Z t \ m = {z tn } n ≠ m represents a class set at time t excluding purchase m.

Figure 0004839358
Figure 0004839358

ステップS4において、クラス確率不変化度推定部215では、各ユーザuのクラス確率不変化度αtuを、尤度が高くなるように計算する。例えば、次の式(8)の計算を繰り返すことで尤度が高いαtuを計算することができる。ここで、「←」は代入、Ψ(・)はディガンマ関数(ガンマ関数の対数の一階導関数)を表す。 In step S4, the class probability unchanged degree estimation unit 215 calculates the class probability unchanged degree α tu of each user u so as to increase the likelihood. For example, α tu with a high likelihood can be calculated by repeating the calculation of the following equation (8). Here, “←” represents substitution, and Ψ (•) represents a digamma function (the first derivative of the logarithm of the gamma function).

Figure 0004839358
Figure 0004839358

ステップS5において、購入確率不変化度推定部216では、各クラスzの購入確率不変化度βtzを、尤度が高くなるように計算する。例えば、次の式(9)の計算を繰り返すことで尤度が高いβtzを計算することができる。 In step S5, the purchase probability unchanged degree estimation unit 216 calculates the purchase probability unchanged degree βtz of each class z so that the likelihood becomes high. For example, β tz having a high likelihood can be calculated by repeating the calculation of the following equation (9).

Figure 0004839358
Figure 0004839358

繰り返しの終了条件(ステップS6)としては、クラス集合44、クラス確率不変化度集合47、購入確率不変化度集合48が収束する、あるいは、繰り返し数がある一定以上になる、などが考えられる。本実施形態では、例えば、前記収束を終了条件とする。   As the repetition end condition (step S6), the class set 44, the class probability invariance set 47, and the purchase probability invariance set 48 converge, or the number of repetitions is more than a certain value. In the present embodiment, for example, the convergence is set as an end condition.

終了条件を満たした後(ステップS6でYes)、パラメータ書き込み部217によって、次の式(10)で計算される時刻tのユーザuのクラスzに属する確率をクラス確率集合45に格納し、次の式(11)で計算される時刻tのクラスzの商品iの購入確率を購入確率集合46に格納する(ステップS7)。   After the end condition is satisfied (Yes in step S6), the parameter writing unit 217 stores the probability belonging to the class z of the user u at time t calculated by the following equation (10) in the class probability set 45, and The purchase probability of the product i of the class z at the time t calculated by the equation (11) is stored in the purchase probability set 46 (step S7).

Figure 0004839358
Figure 0004839358

<興味予測部>
図4を参照しながら、興味予測部22について説明する。図4は、本実施形態の興味予測装置の興味予測部のブロック図を含む説明図である。図4に示すように、興味予測部22は、パラメータ読み込み部221と、興味計算部222と、興味出力部223とを備えている。
<Interest Prediction Department>
The interest prediction unit 22 will be described with reference to FIG. FIG. 4 is an explanatory diagram including a block diagram of an interest prediction unit of the interest prediction apparatus of the present embodiment. As shown in FIG. 4, the interest prediction unit 22 includes a parameter reading unit 221, an interest calculation unit 222, and an interest output unit 223.

興味予測部22では、まず、パラメータ読み込み部221により、モデル推定部21によって計算された最新のクラス確率集合45および購入確率集合46を読み込む。次に、興味計算部222において、クラス確率集合45および購入確率集合46を用いて、次の式(12)に基づき、時刻tにおけるユーザuの商品iに対する興味、つまり、そのユーザがその商品を購入する確率を計算する。   In the interest prediction unit 22, first, the parameter reading unit 221 reads the latest class probability set 45 and purchase probability set 46 calculated by the model estimation unit 21. Next, the interest calculation unit 222 uses the class probability set 45 and the purchase probability set 46, based on the following equation (12), interest of the user u at the time t, that is, the user selects the product. Calculate the probability of buying.

Figure 0004839358
Figure 0004839358

最後に、興味出力部223により、出力手段5に、予測した興味の結果を出力する。   Finally, the interest output unit 223 outputs the predicted interest result to the output unit 5.

(変形例)
興味予測装置1では、前記したように1時刻前のパラメータ(クラス確率集合45および購入確率集合46)のみ用いるのではなく、L時刻前から1時刻前までのパラメータ(所定時間分のパラメータ)を用いて、興味を予測することも可能である。つまり、φtuがl(エル)時刻前のユーザuのクラス確率φt−1uおよびl(エル)時刻前からのユーザuのクラス確率不変化度αtulをパラメータとして持つ、次の式(13)のようなディリクレ分布に従うとする。
(Modification)
In the interest prediction device 1, as described above, not only the parameters one class before (class probability set 45 and purchase probability set 46) are used, but the parameters from L time to one time before (parameters for a predetermined time) are used. It can also be used to predict interest. In other words, φ tu has the class probability φ t−1u of user u before l (el) time and the class probability invariance α tu of user u from before l (el) time as parameters (13 ).

Figure 0004839358
Figure 0004839358

また、θtzも同様に、l(エル)時刻前のクラスzの購入確率θt−1uおよびl(エル)時刻前からのクラスzの購入確率不変化度βtzlをパラメータとして持つ、次の式(14)のようなディリクレ分布に従うとする。 Similarly, θ tz has the following parameters as parameters: purchase probability θ t−1u of class z before l (el) time and purchase probability unchanged degree β tzl of class z from l (el) time before Suppose that it follows a Dirichlet distribution as in equation (14).

Figure 0004839358
Figure 0004839358

L時刻前から1時刻前までのパラメータを用いることにより、1時刻前のパラメータのみを用いる場合に比べ、モデル推定が安定化し、予測精度が高まることが期待できる。この変形例において、モデル推定部21による処理の際、式(7)、式(8)、式(9)、式(10)および式(11)を、それぞれ次のような式(15)、式(16)、式(17)、式(18)および式(19)に変更する。この変形例において、興味予測部22による処理の内容に変更はない。   By using parameters from time L to one hour before, it can be expected that model estimation is stabilized and prediction accuracy is increased as compared to the case of using only the parameters one time before. In this modified example, in the processing by the model estimation unit 21, Expression (7), Expression (8), Expression (9), Expression (10), and Expression (11) are changed into the following expressions (15), It changes into Formula (16), Formula (17), Formula (18), and Formula (19). In this modification, the content of the processing by the interest prediction unit 22 is not changed.

Figure 0004839358
Figure 0004839358

このように、本実施形態の興味予測装置1によれば、複数のユーザが複数のクラスに分けられるという前提でモデル推定部21が確率モデルを推定し、興味予測部22がその確率モデルに基づいてユーザの興味(その後の購買状況(購入確率))を予測することで、入力データ43(時系列の購買データ)の経時的な変化に対応しながらも少ない計算量でユーザの興味を予測することができる。つまり、この予測により、各ユーザに対して適したお勧めの商品を教示したり、どの商品が同じユーザに購入されやすいか、どの商品がいつ人気があったかなどの購買傾向を解析したりすることが可能となる。   Thus, according to the interest prediction apparatus 1 of the present embodiment, the model estimation unit 21 estimates the probability model on the assumption that a plurality of users are divided into a plurality of classes, and the interest prediction unit 22 is based on the probability model. By predicting the user's interest (subsequent purchase status (purchase probability)), the user's interest can be predicted with a small amount of calculation while responding to changes in the input data 43 (time-series purchase data) over time. be able to. In other words, this prediction can be used to teach recommended products suitable for each user, and to analyze purchasing trends such as which products are likely to be purchased by the same user and which products were popular when. Is possible.

また、クラス確率集合45および購入確率集合46が、入力データ43の最新時刻分だけでなく、それ以前の所定時間分を考慮したものであることで、予測の精度を高くすることができる。   Further, since the class probability set 45 and the purchase probability set 46 consider not only the latest time of the input data 43 but also a predetermined time before that, the prediction accuracy can be increased.

さらに、クラス確率集合45および購入確率集合46を用いて表した入力データ43に対する尤度(式(6)参照)を計算し、その尤度の最大化を図ることで確率モデルを推定するので、確率モデルひいては予測の精度を高くすることができる。   Furthermore, the likelihood model is estimated by calculating the likelihood (see Expression (6)) for the input data 43 expressed using the class probability set 45 and the purchase probability set 46, and maximizing the likelihood. The probability model and thus the accuracy of prediction can be increased.

また、興味予測装置1は、一般的なコンピュータに、前記した各処理のプログラムを実行させることで実現することができる。このプログラムは、通信回線を介して提供することが可能であるし、CD−ROM(Compact Disc−Read Only Memory)等の記録媒体に書き込んで配布することも可能である。   In addition, the interest prediction device 1 can be realized by causing a general computer to execute the above-described processing programs. This program can be provided through a communication line, or can be written and distributed on a recording medium such as a CD-ROM (Compact Disc-Read Only Memory).

以上で本実施形態の説明を終えるが、本発明の態様はこれらに限定されるものではない。例えば、本発明は、任意の確率モデルや尤度計算式に適用が可能である。また、本実施形態では、uをユーザ、iを商品としたが、時間発展する離散変数であれば、uおよびiは任意の情報でよい。例えば、uを文書、iを単語とし、各文書でどの単語が使われやすいかの時間発展のモデル化にも適用できる。その他、ハードウェアやフローチャート等の具体的な構成について、本発明の趣旨を逸脱しない範囲で適宜変更が可能である。   Although description of this embodiment is finished above, the aspect of the present invention is not limited to these. For example, the present invention can be applied to any probability model or likelihood calculation formula. In this embodiment, u is a user and i is a product. However, u and i may be arbitrary information as long as they are discrete variables that develop over time. For example, u is a document, i is a word, and it can be applied to modeling of time evolution of which word is likely to be used in each document. In addition, specific configurations such as hardware and flowcharts can be appropriately changed without departing from the spirit of the present invention.

次に、本実施形態の興味予測装置1の実施例について説明する。本実施形態の興味予測装置1を評価するため、購買データの予測実験を行った。単位時間は1日とした。前日までのデータを学習データとして用い、各日の各ユーザの購買を予測する問題の正答率で評価を行った。   Next, an example of the interest prediction device 1 of the present embodiment will be described. In order to evaluate the interest prediction apparatus 1 of the present embodiment, a purchase data prediction experiment was performed. The unit time was 1 day. Using the data up to the previous day as learning data, we evaluated the correct answer rate of the problem of predicting the purchase of each user on each day.

購買データとして、動画および漫画配信のオンラインストアでのログを用いた。動画データは、2007年5月14日から2007年8月31日までの110日間のログで、ユーザ数は70,122、商品数は7,469、全購買数は11,243,935であった。漫画データは、2006年1月1日から2006年5月31日までの151日間のログで、ユーザ数は143,212、商品数は206、全購買数は12,642,505であった、なお、各データともに、購買数が10未満の商品、10未満のユーザは省いた。また、各データともに、新商品が1週間ごとに発売されている。   As purchase data, we used videos and logs from the online store for manga distribution. The video data is a log for 110 days from May 14, 2007 to August 31, 2007. The number of users is 70,122, the number of products is 7,469, and the total number of purchases is 11,243,935. It was. The comic data is a log for 151 days from January 1, 2006 to May 31, 2006. The number of users was 143,212, the number of products was 206, and the total number of purchases was 12,642,505. For each data, products with a purchase quantity of less than 10 were omitted. Also, for each data, new products are released every week.

比較手法として、比較例1(LDAall)、比較例2(LDAonline)、比較例3(LDAone)を用いた。LDAは、本実施形態の興味予測装置1の処理内容からクラス確率集合45および購入確率集合46の時間発展を除いた手法である。比較例1(LDAall)ではログ開始日から前日までの全データを用いてLDAを学習し、比較例2(LDAonline)では前日のデータおよびそれ以前のデータで推定したクラス集合を用いてLDAを学習し、比較例3(LDAone)では前日のデータのみを用いてLDAを学習する。   As comparative methods, Comparative Example 1 (LDAall), Comparative Example 2 (LDAonline), and Comparative Example 3 (LDAone) were used. LDA is a method in which the time evolution of the class probability set 45 and the purchase probability set 46 is removed from the processing contents of the interest prediction device 1 of the present embodiment. In Comparative Example 1 (LDAall), LDA is learned using all data from the log start date to the previous day, and in Comparative Example 2 (LDAonline), LDA is learned using the class set estimated from the previous day's data and previous data. In Comparative Example 3 (LDAone), LDA is learned using only the previous day's data.

図5は、(a)が動画、(b)が漫画について、日毎の3ベスト正答率を示すグラフであり、縦軸が正答率(%)を、横軸が日を示している。3ベスト正答率とは、予測により得られた正解候補の上位3つのいずれかが実際に正解していた場合の割合(%)のことである。また、実施例1(OurMethod7)は、7日前から前日までのパラメータを用いた場合の興味予測装置1による予測結果を示し、実施例2(OurMethod1)は、1日前のパラメータのみを用いた場合の興味予測装置1による予測結果を示す。   FIG. 5 is a graph showing the three best answer rates per day for (a) a moving image and (b) a comic, where the vertical axis shows the correct answer rate (%) and the horizontal axis shows the day. The 3-best correct answer rate is a ratio (%) when any of the top three correct answer candidates obtained by prediction is actually correct. In addition, Example 1 (OurMethod 7) shows a prediction result by the interest prediction device 1 when parameters from 7 days ago to the previous day are used, and Example 2 (OurMethod 1) shows a case where only the parameters of the previous day are used. The prediction result by the interest prediction apparatus 1 is shown.

図5の(a)において、実施例1(OurMethod7)および実施例2(OurMethod1)の3ベスト正答率は、比較例1〜3の3ベスト正答率と比べて、総じて高いことがわかる。また、図5の(b)において、実施例1(OurMethod7)の3ベスト正答率は、比較例1〜3の3ベスト正答率と比べて、総じて高いことがわかる。さらに、図5の(a)および(b)において、特に、実施例1(OurMethod7)、実施例2(OurMethod1)および比較例3(LDAone)について、1週間ごとの谷の部分は、新商品が発売されたことにより正答率が低くなっていることを示している。   In FIG. 5A, it can be seen that the 3 best correct answer rates of Example 1 (OurMethod 7) and Example 2 (Our Method 1) are generally higher than the 3 best correct answer rates of Comparative Examples 1-3. Moreover, in FIG.5 (b), it turns out that the 3 best correct answer rate of Example 1 (OurMethod7) is generally high compared with the 3 best correct answer rate of Comparative Examples 1-3. Further, in FIGS. 5A and 5B, in particular, in Example 1 (OurMethod 7), Example 2 (OurMethod 1), and Comparative Example 3 (LDAone), a new product is shown in the valley portion for each week. It shows that the correct answer rate is low due to the release.

また、その全時間で平均した正答率(%)を表1に示す。表1において、(a)は動画、(b)は漫画の場合の結果である。なお、正答率に付随する括弧書きの数字(「0.9」など)は標準偏差を表す。(a)および(b)の表において、最左列の「1」〜「5」は、「1ベスト正答率」〜「5ベスト正答率」を表している。(a)および(b)の表において、実施例1(OurMethod7)の正答率が最も高く、正確にユーザの興味を予測できていることが確認できる。実施例2(OurMethod1)の正答率も、動画ではいずれの比較例よりも高く、漫画では比較例1とは同等であるが比較例2,3よりは高い。   In addition, Table 1 shows the correct answer rate (%) averaged over the entire time. In Table 1, (a) is a moving image, and (b) is a result of a comic. A number in parentheses (such as “0.9”) accompanying the correct answer rate represents a standard deviation. In the tables of (a) and (b), “1” to “5” in the leftmost column represent “1 best correct answer rate” to “5 best correct answer rate”. In the tables of (a) and (b), it can be confirmed that the correct answer rate of Example 1 (OurMethod 7) is the highest and the user's interest can be accurately predicted. The correct answer rate of Example 2 (OurMethod 1) is also higher than any of the comparative examples in the moving image, and is equal to that of the comparative example 1 in the comic, but higher than those of the comparative examples 2 and 3.

Figure 0004839358
Figure 0004839358

また、この場合の、Xeon(登録商標) 2.93GHzのCPUのPC(Personal Computer)における各手法の1時刻あたりの平均計算時間(秒)を表2に示す。表2において、(a)は動画、(b)は漫画の場合の結果である。なお、平均計算時間の値に付随する括弧書きの数字(「40.4」など)は標準偏差を表す。比較例1は、漫画データにおける正答率(表1の(b)参照)が実施例1および実施例2とほぼ同等であった。この比較例1は、全データを用いるため、動画および漫画の場合ともに、他手法と比べ、多くの計算量が必要となっていることがわかる。一方、動画および漫画の場合ともに、実施例1および実施例2では、比較例1と比べ、1/10以下の時間で計算が終っており、大規模データでも高速に計算できていることがわかる。   Table 2 shows the average calculation time (seconds) per time of each method in a PC (Personal Computer) of a CPU of Xeon (registered trademark) 2.93 GHz in this case. In Table 2, (a) is the result for a moving image, and (b) is the result for a comic. The numbers in parentheses (such as “40.4”) attached to the average calculation time value represent standard deviation. In Comparative Example 1, the correct answer rate in comic data (see (b) of Table 1) was almost equal to that in Example 1 and Example 2. Since this comparative example 1 uses all data, it can be seen that a large amount of calculation is required for both moving images and comics compared to other methods. On the other hand, in both the case of moving images and comics, the calculation in Example 1 and Example 2 was completed in 1/10 or less time compared to Comparative Example 1, and it can be seen that high-speed calculation was possible even for large-scale data. .

Figure 0004839358
Figure 0004839358

このように、本実施例から、本実施形態の興味予測装置1による計算は、計算精度と計算時間の2つの観点から見た場合、比較例1〜3の計算と比べて優れていることがわかる。   Thus, from this example, the calculation by the interest prediction device 1 of the present embodiment is superior to the calculations of Comparative Examples 1 to 3 when viewed from two viewpoints of calculation accuracy and calculation time. Recognize.

本実施形態の興味予測装置の構成を示すブロック図である。It is a block diagram which shows the structure of the interest prediction apparatus of this embodiment. 本実施形態の興味予測装置のモデル推定部のブロック図を含む説明図である。It is explanatory drawing including the block diagram of the model estimation part of the interest prediction apparatus of this embodiment. モデル推定部の処理の流れを示すフローチャートである。It is a flowchart which shows the flow of a process of a model estimation part. 本実施形態の興味予測装置の興味予測部のブロック図を含む説明図である。It is explanatory drawing containing the block diagram of the interest prediction part of the interest prediction apparatus of this embodiment. 本実施形態の興味予測方法を含む各手法の日毎の3ベスト正答率を示すグラフである。It is a graph which shows the 3-best correct answer rate for every method of each method including the interest prediction method of this embodiment.

符号の説明Explanation of symbols

1 興味予測装置(購入予測装置)
2 演算手段
3 入力手段
4 記憶手段
5 出力手段
11 バスライン
21 モデル推定部
22 興味予測部(購入予測部)
23 メモリ
40a プログラム格納部
40b データ格納部
41 モデル推定プログラム
42 興味予測プログラム
43 入力データ
44 クラス集合
45 クラス確率集合
46 購入確率集合
47 クラス確率不変化度集合
48 購入確率不変化度集合
1 Interest prediction device (purchase prediction device)
2 Calculation means 3 Input means 4 Storage means 5 Output means 11 Bus line 21 Model estimation part 22 Interest prediction part (purchase prediction part)
23 memory 40a program storage unit 40b data storage unit 41 model estimation program 42 interest prediction program 43 input data 44 class set 45 class probability set 46 purchase probability set 47 class probability invariant set 48 purchase probability invariant set

Claims (8)

複数のユーザの複数の商品に対する購買のデータが時系列に沿って予め集められた購買データに基づいて、任意の前記ユーザの任意の前記商品に対するその後の購入の可能性を予測する場合に、前記複数のユーザを複数のクラスに分け、同じ前記クラスに属する前記ユーザは同じ確率で任意の前記商品を購入すると仮定する確率モデルを用いて前記購入の可能性を予測する購入予測装置による購入予測方法であって、
前記購入予測装置は、
前記購買データ、
前記ユーザごとの前記クラスそれぞれに属する確率の集合であるクラス確率集合、
前記ユーザの前記クラスに属する確率が変わらない度合いを示すクラス確率不変化度をパラメータとして持つ、前記ユーザのクラス確率の事前確率、
前記クラスごとの前記商品それぞれがそのクラスのユーザによって購入される確率の集合である購入確率集合、および、
前記クラスに属する前記ユーザによって商品を購入される確率が変わらない度合いを示す購入確率不変化度をパラメータとして持つ、前記クラスの購入確率の事前確率、を記憶する記憶部と、
モデル推定部と、購入予測部と、を備えており、
前記モデル推定部は、前記クラス確率不変化度の尤度と前記購入確率不変化度の尤度とが最大になるように、前記クラス確率不変化度と前記購入確率不変化度とを繰り返し計算し、その計算した前記クラス確率不変化度と前記購入確率不変化度とを用いて、前記クラス確率集合および前記購入確率集合を更新し
前記購入予測部は、前記更新した前記クラス確率集合および前記購入確率集合を用いて、前記ユーザが前記商品を購入する確率を算出す
とを特徴とする購入予測方法。
When predicting the possibility of subsequent purchase for any of the products of any user based on purchase data collected in advance in time series for purchase data for a plurality of products of a plurality of users, A purchase prediction method by a purchase prediction device that predicts the purchase possibility using a probability model that divides a plurality of users into a plurality of classes and assumes that the users belonging to the same class purchase any of the products with the same probability Because
The purchase prediction device
The purchase data,
A class probability set that is a set of probabilities belonging to each of the classes for each user;
A prior probability of the class probability of the user, having as a parameter a class probability invariance indicating the degree to which the probability that the user belongs to the class does not change,
A purchase probability set that is a set of probabilities that each of the products for each class is purchased by a user of that class; and
A storage unit that stores , as a parameter, a purchase probability invariance indicating a degree that the probability of purchasing a product by the user belonging to the class does not change, and a prior probability of the purchase probability of the class ;
A model estimation unit and a purchase prediction unit,
The model estimation unit repeatedly calculates the class probability unchanged degree and the purchase probability unchanged degree so that the likelihood of the class probability unchanged degree and the likelihood of the purchase probability unchanged degree are maximized. and, using said purchase probability not change degree and the class probability invariable degree that the calculation, further new the class probability set and the purchase probability set,
The purchase prediction unit, the updated using the class probability set and the purchase probability set, you calculate the probability that before Symbol user purchases the merchandise
Purchase prediction method wherein a call.
前記ユーザのクラス確率の事前確率は、ディリクレ分布として定義されており、The prior probability of the user class probability is defined as a Dirichlet distribution,
前記クラスの購入確率の事前確率は、ディリクレ分布として定義されているPrior probability of purchase probability of the class is defined as Dirichlet distribution
ことを特徴とする請求項1に記載の購入予測方法。The purchase prediction method according to claim 1, wherein:
前記記憶部は、前記購買データの最新時刻以前の所定時間分の前記クラス確率集合、前記ユーザのクラス確率の事前確率、前記購入確率集合、および、前記クラスの購入確率の事前確率、を記憶しており、
前記モデル推定部は、前記所定時間分の前記クラス確率不変化度の尤度と、前記所定時間分の前記購入確率不変化度の尤度と、が最大になるように、前記所定時間分の前記クラス確率不変化度と、前記所定時間分の前記購入確率不変化度と、を繰り返し計算し、その計算した前記所定時間分の前記クラス確率不変化度と、前記所定時間分の前記購入確率不変化度と、を用いて、前記所定時間分の前記クラス確率集合、および、前記所定時間分の前記購入確率集合を更新し、
前記購入予測部は、前記更新した前記所定時間分の前記クラス確率集合、および、前記所定時間分の前記購入確率集合を用いて、前記ユーザが前記商品を購入する確率を算出する
ことを特徴とする請求項1または請求項2に記載の購入予測方法。
The storage unit, of the purchase latest time before the predetermined period of time data, the class probability set, the prior probability of class probability of the user, the purchase probability set, and store the prior probability, the probability purchase of the class And
The model estimation unit is configured to increase the likelihood of the class probability unchanged for the predetermined time and the likelihood of the purchase probability unchanged for the predetermined time to the maximum. The class probability unchanged degree and the purchase probability unchanged degree for the predetermined time are repeatedly calculated, the class probability unchanged degree for the calculated predetermined time and the purchase probability for the predetermined time are calculated. Update the class probability set for the predetermined time and the purchase probability set for the predetermined time using the degree of invariance,
The purchase prediction unit calculates the probability that the user purchases the product using the updated class probability set for the predetermined time and the purchase probability set for the predetermined time. The purchase prediction method according to claim 1 or 2 .
複数のユーザの複数の商品に対する購買のデータが時系列に沿って予め集められた購買データに基づいて、任意の前記ユーザの任意の前記商品に対するその後の購入の可能性を予測する場合に、前記複数のユーザを複数のクラスに分け、同じ前記クラスに属する前記ユーザは同じ確率で任意の前記商品を購入すると仮定する確率モデルを用いて前記購入の可能性を予測する購入予測装置であって、
前記購買データ、
前記ユーザごとの前記クラスそれぞれに属する確率の集合であるクラス確率集合、
前記ユーザの前記クラスに属する確率が変わらない度合いを示すクラス確率不変化度をパラメータとして持つ、前記ユーザのクラス確率の事前確率、
前記クラスごとの前記商品それぞれがそのクラスのユーザによって購入される確率の集合である購入確率集合、および、
前記クラスに属する前記ユーザによって商品を購入される確率が変わらない度合いを示す購入確率不変化度をパラメータとして持つ、前記クラスの購入確率の事前確率、を記憶する記憶部と、
前記クラス確率不変化度の尤度と前記購入確率不変化度の尤度とが最大になるように、前記クラス確率不変化度と前記購入確率不変化度とを繰り返し計算し、その計算した前記クラス確率不変化度と前記購入確率不変化度とを用いて、前記クラス確率集合および前記購入確率集合を更新するモデル推定部と、
前記更新した前記クラス確率集合および前記購入確率集合を用いて、前記ユーザが前記商品を購入する確率を算出する購入予測部と、
を備えることを特徴とする購入予測装置。
When predicting the possibility of subsequent purchase for any of the products of any user based on purchase data collected in advance in time series for purchase data for a plurality of products of a plurality of users, A purchase prediction apparatus that divides a plurality of users into a plurality of classes and predicts the possibility of purchase using a probability model that assumes that the users belonging to the same class purchase any of the products with the same probability,
The purchase data,
A class probability set that is a set of probabilities belonging to each of the classes for each user;
A prior probability of the class probability of the user, having as a parameter a class probability invariance indicating the degree to which the probability that the user belongs to the class does not change,
A purchase probability set that is a set of probabilities that each of the products for each class is purchased by a user of that class; and
A storage unit that stores , as a parameter, a purchase probability invariance indicating a degree that the probability of purchasing a product by the user belonging to the class does not change, and a prior probability of the purchase probability of the class ;
The class probability unchanged and the purchase probability unchanged are repeatedly calculated so that the likelihood of the class probability unchanged and the purchase probability unchanged are maximized, and the calculated using class probability not change degree and the said purchase probability immutable degree, a model estimation unit you update the class probability set and the purchase probability set,
Using the class probability set and the purchase probability set was the update, and the purchase prediction unit that to calculate the probability that before Symbol user purchases the merchandise,
A purchase prediction apparatus comprising:
前記ユーザのクラス確率の事前確率は、ディリクレ分布として定義されており、The prior probability of the user class probability is defined as a Dirichlet distribution,
前記クラスの購入確率の事前確率は、ディリクレ分布として定義されている  Prior probability of purchase probability of the class is defined as Dirichlet distribution
ことを特徴とする請求項4に記載の購入予測装置。  The purchase prediction apparatus according to claim 4, wherein:
前記記憶部は、前記購買データの最新時刻以前の所定時間分の前記クラス確率集合、前記ユーザのクラス確率の事前確率、前記購入確率集合、および、前記クラスの購入確率の事前確率、を記憶しており、
前記モデル推定部は、前記所定時間分の前記クラス確率不変化度の尤度と、前記所定時間分の前記購入確率不変化度の尤度と、が最大になるように、前記所定時間分の前記クラス確率不変化度と、前記所定時間分の前記購入確率不変化度と、を繰り返し計算し、その計算した前記所定時間分の前記クラス確率不変化度と、前記所定時間分の前記購入確率不変化度と、を用いて、前記所定時間分の前記クラス確率集合、および、前記所定時間分の前記購入確率集合を更新し、
前記購入予測部は、前記更新した前記所定時間分の前記クラス確率集合、および、前記所定時間分の前記購入確率集合を用いて、前記ユーザが前記商品を購入する確率を算出する
ことを特徴とする請求項4または請求項5に記載の購入予測装置。
The storage unit, of the purchase latest time before the predetermined period of time data, the class probability set, the prior probability of class probability of the user, the purchase probability set, and stores the prior probability, the probability purchase of the class And
The model estimation unit is configured to increase the likelihood of the class probability unchanged for the predetermined time and the likelihood of the purchase probability unchanged for the predetermined time to the maximum. The class probability unchanged degree and the purchase probability unchanged degree for the predetermined time are repeatedly calculated, the class probability unchanged degree for the calculated predetermined time and the purchase probability for the predetermined time are calculated. Update the class probability set for the predetermined time and the purchase probability set for the predetermined time using the degree of invariance,
The purchase prediction unit calculates the probability that the user purchases the product using the updated class probability set for the predetermined time and the purchase probability set for the predetermined time. The purchase prediction device according to claim 4 or 5 .
コンピュータを請求項4から請求項6のいずれか一項に記載の購入予測装置の各部として機能させるための購入予測プログラム。 The purchase prediction program for functioning a computer as each part of the purchase prediction apparatus as described in any one of Claims 4-6 . 請求項に記載の購入予測プログラムが記録されたことを特徴とするコンピュータに読み取り可能な記録媒体。 A computer-readable recording medium on which the purchase prediction program according to claim 7 is recorded.
JP2008227235A 2008-09-04 2008-09-04 Purchase prediction method, purchase prediction apparatus, purchase prediction program, and recording medium Expired - Fee Related JP4839358B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2008227235A JP4839358B2 (en) 2008-09-04 2008-09-04 Purchase prediction method, purchase prediction apparatus, purchase prediction program, and recording medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2008227235A JP4839358B2 (en) 2008-09-04 2008-09-04 Purchase prediction method, purchase prediction apparatus, purchase prediction program, and recording medium

Publications (2)

Publication Number Publication Date
JP2010061456A JP2010061456A (en) 2010-03-18
JP4839358B2 true JP4839358B2 (en) 2011-12-21

Family

ID=42188178

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2008227235A Expired - Fee Related JP4839358B2 (en) 2008-09-04 2008-09-04 Purchase prediction method, purchase prediction apparatus, purchase prediction program, and recording medium

Country Status (1)

Country Link
JP (1) JP4839358B2 (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102346894B (en) * 2010-08-03 2017-03-01 阿里巴巴集团控股有限公司 The output intent of recommendation information, system and server
JP5775425B2 (en) * 2011-11-18 2015-09-09 日本電信電話株式会社 Purchase data analysis apparatus, method, and program
JP5860795B2 (en) * 2012-11-06 2016-02-16 日本電信電話株式会社 Survival rate estimation apparatus, method, and program
JP6078014B2 (en) * 2014-02-27 2017-02-08 日本電信電話株式会社 Purchase motivation learning apparatus, purchase prediction apparatus, method, and program
KR101813805B1 (en) * 2016-09-28 2017-12-29 한양대학교 산학협력단 Method and Apparatus for purchase probability prediction of user using machine learning
CN110796280B (en) * 2018-08-02 2024-08-20 北京京东尚科信息技术有限公司 User behavior prediction method, device and computer-readable storage medium
CN111091416A (en) * 2019-12-10 2020-05-01 北京云迹科技有限公司 Method and device for predicting probability of hotel purchase robot

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002279279A (en) * 2001-03-22 2002-09-27 Just Syst Corp Product recommendation device, product recommendation method, and product recommendation program
JP2004054439A (en) * 2002-07-17 2004-02-19 Nippon Telegr & Teleph Corp <Ntt> Marketing strategy formulation support method, device, and program based on customer behavior analysis
JP4751242B2 (en) * 2006-05-29 2011-08-17 日本電信電話株式会社 RECOMMENDATION DEVICE, RECOMMENDATION METHOD, RECOMMENDATION PROGRAM, AND RECORDING MEDIUM CONTAINING THE PROGRAM

Also Published As

Publication number Publication date
JP2010061456A (en) 2010-03-18

Similar Documents

Publication Publication Date Title
JP4839358B2 (en) Purchase prediction method, purchase prediction apparatus, purchase prediction program, and recording medium
US11935049B2 (en) Graph data processing method and apparatus, computer device, and storage medium
Varoquaux et al. Scikit-learn: Machine learning without learning the machinery
WO2019137104A1 (en) Recommendation method and device employing deep learning, electronic apparatus, medium, and program
Rytgaard et al. Estimation of time‐specific intervention effects on continuously distributed time‐to‐event outcomes by targeted maximum likelihood estimation
JP4847916B2 (en) RECOMMENDATION DEVICE, RECOMMENDATION METHOD, RECOMMENDATION PROGRAM, AND RECORDING MEDIUM CONTAINING THE PROGRAM
US20240256871A1 (en) Finite rank deep kernel learning with linear computational complexity
CN112801706B (en) Game user behavior data mining method and system
US20200143252A1 (en) Finite rank deep kernel learning for robust time series forecasting and regression
US20210042820A1 (en) Extending finite rank deep kernel learning to forecasting over long time horizons
JP7655328B2 (en) Information processing system, information processing method, and computer program
CN108073464A (en) A kind of time series data abnormal point detecting method and device based on speed and acceleration
Jimenez et al. An approximate innovation method for the estimation of diffusion processes from discrete data
JP5560154B2 (en) Model parameter estimation apparatus and program thereof
EP4258183A1 (en) Partial action interval estimation model construction device, partial action interval estimation model construction method, and partial action interval estimation model construction program
Duan et al. Sequential Monte Carlo optimization and statistical inference
JP5210224B2 (en) Annotation data analysis device, annotation data analysis program, and recording medium recording the program
US12475132B2 (en) Computing system and method for applying monte carlo estimation to determine the contribution of dependent input variable groups on the output of a data science model
CN111309706A (en) Model training method and device, readable storage medium and electronic equipment
Martin et al. Bayesian Forecasting in the 21st Century: A Modern Review
CN114418060A (en) Identity keeping confrontation training method, device and medium based on graph representation learning
Zammit et al. A simplified and numerically stable approach to the BG/NBD churn prediction model
Lindström et al. Optimal adaptive sequential calibration of option models
CN113590960A (en) User recognition model training method, electronic device, and computer storage medium
JP7709662B2 (en) Specification-compliant data estimation device, machine learning method, specification-compliant data estimation method, and program

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20100721

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110712

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20110811

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20110826

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20110920

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20111003

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20141007

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

LAPS Cancellation because of no payment of annual fees