JPWO2017222030A1 - Event prediction apparatus and event prediction method - Google Patents

Event prediction apparatus and event prediction method Download PDF

Info

Publication number
JPWO2017222030A1
JPWO2017222030A1 JP2018524169A JP2018524169A JPWO2017222030A1 JP WO2017222030 A1 JPWO2017222030 A1 JP WO2017222030A1 JP 2018524169 A JP2018524169 A JP 2018524169A JP 2018524169 A JP2018524169 A JP 2018524169A JP WO2017222030 A1 JPWO2017222030 A1 JP WO2017222030A1
Authority
JP
Japan
Prior art keywords
function
event
prediction
state
time
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.)
Granted
Application number
JP2018524169A
Other languages
Japanese (ja)
Other versions
JP6935928B2 (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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Publication of JPWO2017222030A1 publication Critical patent/JPWO2017222030A1/en
Application granted granted Critical
Publication of JP6935928B2 publication Critical patent/JP6935928B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Marketing (AREA)
  • Game Theory and Decision Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

【課題】過去に発生した特定事象の履歴データに基づいて特定事象の将来の発生密度を予測する事象発生予測装置において、データ数が少ない場合も予測精度を上げる。
【解決手段】事象発生予測装置は、予測式構築部10と、予測部30と、を備える。予測式構築部10は、特定事象の発生時刻tおよび特定事象の発生領域を特定する領域特定変数xの関数ρ(t,x)で前記特定事象の発生密度が与えられ、かつ、関数ρ(t,x)は外部因子{f}と関数ρ(t,x)の写像F[ρ(t,x)+{f}]で与えられるとして、写像F[ρ(t,x)+{f}]を過去に発生した特定事象の履歴データから求めて、関数ρ(t,x)を発生時刻tおよび領域特定変数xの関数として表現する。予測部30は、関数ρに将来の時刻および領域を特定する値を入力して、特定事象の発生密度を予測する予測部と、を有することを特徴とする。
【選択図】図1
In an event occurrence prediction apparatus for predicting the future occurrence density of a specific event based on historical data of a specific event that occurred in the past, the prediction accuracy is improved even when the number of data is small.
An event occurrence prediction apparatus includes a prediction formula construction unit 10 and a prediction unit 30. The prediction formula construction unit 10 is given the occurrence density of the specific event by the function ρ (t, x) of the region specific variable x that specifies the occurrence time t of the specific event and the region of occurrence of the specific event, and the function ρ ( t, x) is given by the external factor {f} and the map F [ρ (t, x) + {f}] of the function ρ (t, x), and the map F [ρ (t, x) + {f }] Is obtained from the history data of the specific event that occurred in the past, and the function ρ (t, x) is expressed as a function of the occurrence time t and the region specific variable x. The prediction unit 30 includes a prediction unit that inputs a value specifying a future time and region to the function ρ and predicts the occurrence density of a specific event.
[Selection] Figure 1

Description

本発明は、事象予測装置および事象予測方法に関する。   The present invention relates to an event prediction apparatus and an event prediction method.

米国カリフォルニア州で犯罪発生の可能性が高い要注意エリアを警察官が重点的に見回ることにより、その後の犯罪発生が抑制されたという事例がある。また、地震の余震予測のために提案されたアルゴリズムを用いて犯罪の発生密度の予測精度を向上させた例がある。   In some cases, police officers focused on areas requiring attention in California, where crimes are likely to occur, and subsequent crimes were suppressed. In addition, there is an example in which the prediction accuracy of crime occurrence density is improved by using an algorithm proposed for predicting aftershocks of an earthquake.

犯罪発生の予測は、たとえばSelf−Exciting Point Process(SEPP)モデルを用いて行われる。SEPPモデルでは、特定の位置、将来の特定の時刻における犯罪の発生密度を、過去の犯罪発生事象による影響の総和で示す。より具体的には、ある過去の犯罪発生事象の犯罪の特定の位置および特定の時刻における発生密度に対する影響は、予測する位置と過去の犯罪の発生位置との距離と、予測する時刻と過去の犯罪の発生時刻との間の時間との関数で表されるとする。さらに、過去に発生した犯罪のそれぞれについて求めた発生密度に対する影響を足し合わせることによって求める。   Crime occurrence is predicted using, for example, a Self-Exciting Point Process (SEPP) model. In the SEPP model, the crime occurrence density at a specific position and a specific time in the future is indicated by the sum of the effects of past crime occurrence events. More specifically, the influence of a certain crime occurrence event on the occurrence density of a crime at a specific position and a specific time is determined by the distance between the predicted position and the past crime occurrence position, the predicted time and the past It is assumed that it is expressed as a function of the time between the crime occurrence time. Furthermore, it calculates | requires by adding the influence with respect to the occurrence density calculated | required about each crime which occurred in the past.

過去の犯罪発生事象の犯罪発生密度に対する影響を、現在着目する位置・時刻と過去の事象の位置・時刻との間の距離Δxおよび時間Δt(それぞれ適切な長さ、時間を単位にスケールする)を用いてg(Δt,Δx)と書くことにすると、g(Δt,Δx)=1/((1+Δt)(1+Δx))と表現する方法がある(Prospective Hotspot Method)。期待値最大化法(Expectation Maximization Algorithm)を用いて、過去の犯罪発生事象の履歴データから過去の犯罪発生事象の犯罪の位置xおよび時刻tにおける発生密度に対する影響g(Δt,Δx)を構築する方法もある。   The distance Δx and the time Δt between the position / time of the current focus and the position / time of the past event (the scale is set to an appropriate length and time, respectively). When g (Δt, Δx) is written using, there is a method of expressing as g (Δt, Δx) = 1 / ((1 + Δt) (1 + Δx)) (Prospective Hotspot Method). By using the expectation maximization algorithm, the past crime occurrence event history data is used to construct an influence g (Δt, Δx) of the past crime occurrence event on the crime location x and the occurrence density at time t. There is also a method.

米国特許第8949164号U.S. Pat. No. 8,949,164 米国特許第9129219号U.S. Patent No. 9129219

過去の発生事象である犯罪との距離差Δxおよび時間差Δtにおける発生密度に対する影響g(Δt,Δx)を特定の関数で表現した場合、現実とは異なり、精度が上がらない場合がある。また、期待値最大化法を用いて影響g(Δt,Δx)を構築する場合、機械学習による予測においては、データ数が少ない場合に精度が高く出ない場合がある。   When the influence g (Δt, Δx) on the occurrence density in the distance difference Δx and the time difference Δt from the crime that is a past occurrence event is expressed by a specific function, the accuracy may not be increased unlike in reality. In addition, when the influence g (Δt, Δx) is constructed using the expected value maximization method, in the prediction by machine learning, the accuracy may not be high when the number of data is small.

そこで、本発明は、過去に発生した特定事象の履歴データに基づいて特定事象の将来の発生密度を予測する事象発生予測装置において、データ数が少ない場合も予測精度を上げることを目的とする。   Therefore, an object of the present invention is to increase the prediction accuracy even when the number of data is small in an event occurrence prediction device that predicts the future occurrence density of a specific event based on historical data of a specific event that occurred in the past.

上述の目的を達成するため、本発明は、nを自然数としたときに過去に発生したn種の特定事象の状態を表す値を予測する事象発生装置において、i番目の特定事象の時刻tの状態が関数ρ_i(t)で与えられ、かつ、前記関数ρ_i(t)の関数ベクトルΡ(t)={ρ_1(t),…,ρ_i(t),…,ρ_n(t)}は外部因子ベクトル{f}と前記関数ベクトルΡ(t)の写像F[Ρ(t)+{f}]で与えられるとして、前記写像F[Ρ(t)+{f}]を前記特定事象の状態を表す値と外部因子の履歴データから求める予測式構築部と、前記関数ベクトルΡ(t)に将来の時刻を特定する値を入力して、前記特定事象の発生密度を予測する予測部と、を有することを特徴とする。   In order to achieve the above-described object, the present invention provides an event generator that predicts a value representing the state of n types of specific events that occurred in the past when n is a natural number. The state is given by a function ρ_i (t), and the function vector Ρ (t) = {ρ_1 (t),..., Ρ_i (t), ..., ρ_n (t)} of the function ρ_i (t) is an external factor. Assuming that the vector {f} and the map F [t (t) + {f}] of the function vector Ρ (t) are given, the map F [Ρ (t) + {f}] is expressed as the state of the specific event. A prediction formula construction unit that is obtained from a value to be represented and history data of external factors, and a prediction unit that inputs a value that specifies a future time in the function vector Ρ (t) and predicts the occurrence density of the specific event. It is characterized by having.

また、本発明は、1つの事象が次に引き続く新たな事象に影響を与えるカスケード現象一般に対し、過去の特定事象の状態を表す値の履歴データに基づいて前記特定事象の将来の状態を表す値を予測する事象予測方法において、時刻の関数Ρで前記特定事象の状態を表す値が与えられ、かつ、前記関数Ρは前記外部因子{f}と前記関数Ρの写像F[Ρ+{f}]で与えられるとして、前記写像F[Ρ+{f}]を過去の前記特定事象の状態を表す値の履歴データから求めて、前記関数Ρを前記時刻の関数として表現する予測式構築ステップと、前記関数Ρに将来の時刻および領域を特定する値を入力して、前記特定事象の状態を表す値を予測する予測ステップと、を有することを特徴とする。   Further, the present invention relates to a cascade phenomenon in which one event affects a subsequent new event, and a value representing the future state of the specific event based on historical data of a value representing the state of the specific event in the past. In the event prediction method for predicting, a value representing the state of the specific event is given by a function Ρ of time, and the function Ρ is a map F [Ρ + {f}] of the external factor {f} and the function Ρ. The mapping F [Ρ + {f}] is obtained from the history data of the value representing the state of the specific event in the past, and the prediction formula construction step of expressing the function と し て as a function of the time, And a prediction step of predicting a value representing the state of the specific event by inputting a value specifying a future time and region into the function Ρ.

また、本発明は、1つの事象が次に引き続く新たな事象に影響を与えるカスケード現象一般に対し、過去の特定事象の状態を表す値の履歴データに基づいて前記特定事象の将来の状態を表す値を予測する事象予測プログラムにおいて、計算機に時刻の関数Ρで前記特定事象の状態を表す値が与えられ、かつ、前記関数Ρは前記外部因子{f}と前記関数Ρの写像F[Ρ+{f}]で与えられるとして、前記写像F[Ρ+{f}]を過去の前記特定事象の状態を表す値の履歴データから求めて、前記関数Ρを前記時刻の関数として表現させる予測式構築ステップと、前記関数Ρに将来の時刻および領域を特定する値を入力して、前記特定事象の状態を表す値を予測させる予測ステップと、を実行させることを特徴とする。   Further, the present invention relates to a cascade phenomenon in which one event affects a subsequent new event, and a value representing the future state of the specific event based on historical data of a value representing the state of the specific event in the past. In the event prediction program for predicting, the computer is given a value representing the state of the specific event by a function の of time, and the function Ρ is a map F [Ρ + {f of the external factor {f} and the function Ρ. }], The map F [Ρ + {f}] is obtained from the history data of the value representing the state of the specific event in the past, and the prediction formula construction step for expressing the function と し て as a function of the time; A function of inputting a value specifying a future time and region to the function Ρ and predicting a value representing the state of the specific event is executed.

また、本発明は、1つの事象が次に引き続く新たな事象に影響を与えるカスケード現象一般に対し、過去の特定事象の状態を表す値の履歴データに基づいて前記特定事象の将来の状態を表す値を予測する事象予測システムにおいて、時刻の関数Ρで前記特定事象の状態を表す値が与えられ、かつ、前記関数Ρは前記外部因子{f}と前記関数Ρの写像F[Ρ+{f}]で与えられるとして、前記写像F[Ρ+{f}]を過去の前記特定事象の状態を表す値の履歴データから求めて、前記関数Ρを前記時刻の関数として表現させる予測式構築部と、前記関数Ρに将来の時刻および領域を特定する値を入力して、前記特定事象の状態を表す値を予測する予測部と、前記特定事象の状態に関する値を前記履歴データとして送信する携帯端末及びそれに準じる通信機能を持った端末と、前記携帯端末から新たな前記履歴データを受信してして新たな前記履歴データが得られたときに前記予測式構築部に前記写像F[Ρ+{f}]を再度求めさせるサーバーと、を有することを特徴とする。   Further, the present invention relates to a cascade phenomenon in which one event affects a subsequent new event, and a value representing the future state of the specific event based on historical data of a value representing the state of the specific event in the past. In the event prediction system for predicting the function, a value representing the state of the specific event is given by a function Ρ of time, and the function Ρ is a map F [Ρ + {f}] of the external factor {f} and the function Ρ. A prediction formula construction unit that obtains the map F [Ρ + {f}] from historical data of values representing the state of the specific event in the past, and expresses the function と し て as a function of the time; A value that specifies a future time and region is input to the function 、, a prediction unit that predicts a value representing the state of the specific event, a mobile terminal that transmits a value related to the state of the specific event as the history data, and Associate When the new history data is obtained by receiving new history data from the mobile terminal and the terminal having the communication function, the map F [Ρ + {f}] is sent to the prediction formula construction unit. And a server for requesting again.

また、本発明は、nを自然数としたときに過去に発生したn種の特定事象の状態を表す値を予測する事象発生装置であって、i番目の特定事象の時刻tの状態が関数ρ_i(t)で与えられ、かつ、前記関数ρ_i(t)の関数ベクトルΡ(t)={ρ_1(t),…,ρ_i(t),…,ρ_n(t)}は外部因子ベクトル{f}と前記関数ベクトルΡ(t)の写像F[Ρ(t)+{f}]で与えられるとして、前記写像F[Ρ(t)+{f}]を前記特定事象の状態を表す値と外部因子の履歴データから求める予測式構築部と、前記関数ベクトルΡ(t)に将来の時刻を特定する値を入力して、前記特定事象の発生密度を予測する予測部と、を備えた事象予測装置とともに用いる携帯端末にインストールされるプログラムにおいて、携帯端末に、前記事象予測装置に前記特定事象の状態を表す値を前記履歴データとして送信させる、ことを特徴とする。   The present invention is also an event generating device that predicts a value representing the state of n types of specific events that occurred in the past when n is a natural number, and the state of the i-th specific event at time t is a function ρ_i. And the function vector Ρ (t) = {ρ_1 (t),..., Ρ_i (t),..., Ρ_n (t)} of the function ρ_i (t) is an external factor vector {f} And the map F [Ρ (t) + {f}] of the function vector Ρ (t) and the value representing the state of the specific event and the external An event prediction comprising: a prediction formula construction unit obtained from factor history data; and a prediction unit that inputs a value specifying a future time into the function vector Ρ (t) and predicts the occurrence density of the specific event. In a program installed on a mobile terminal used with the device, To transmit the state values of the specific event to the event prediction device as the history data, characterized in that.

本発明によれば、過去に発生した特定事象の履歴データに基づいて特定事象の将来の発生密度を予測する事象発生予測装置において、データ数が少ない場合も予測精度を上げることができる。   According to the present invention, in an event occurrence prediction device that predicts the future occurrence density of a specific event based on historical data of a specific event that occurred in the past, the prediction accuracy can be increased even when the number of data is small.

本発明に係る事象発生予測装置の一実施の形態におけるブロック図である。It is a block diagram in one embodiment of an event occurrence prediction device according to the present invention. 本発明に係る事象発生予測装置の一実施の形態を用いた事象発生予測方法のフローチャートである。It is a flowchart of the event occurrence prediction method using one embodiment of the event occurrence prediction device according to the present invention. 本発明に係る事象発生予測装置の一実施の形態を用いた犯罪発生率予測におけるカーネル関数gの時刻t依存性の評価例を示す図である。It is a figure which shows the example of evaluation of the time t dependence of the kernel function g in the crime occurrence rate prediction using one embodiment of the event occurrence prediction device according to the present invention. 本発明に係る事象発生予測装置の一実施の形態を用いた犯罪発生率予測におけるカーネル関数gの領域x依存性の評価例を示す図である。It is a figure which shows the evaluation example of the area | region x dependence of the kernel function g in the crime occurrence rate prediction using one embodiment of the event occurrence prediction device according to the present invention. 本発明に係る事象発生予測装置の一実施の形態を用いた犯罪発生率予測結果を示す等高線である。It is a contour line which shows the crime occurrence rate prediction result using one embodiment of the event occurrence prediction device according to the present invention.

本発明に係る事象発生予測装置の一実施の形態を、図面を参照して説明する。なお、この実施の形態は単なる例示であり、本発明はこれに限定されない。同一または類似の構成には同一の符号を付し、重複する説明は省略する。   An embodiment of an event occurrence prediction apparatus according to the present invention will be described with reference to the drawings. This embodiment is merely an example, and the present invention is not limited to this. The same or similar components are denoted by the same reference numerals, and redundant description is omitted.

図1は、本発明に係る事象発生予測装置の一実施の形態におけるブロック図である。   FIG. 1 is a block diagram of an embodiment of an event occurrence prediction apparatus according to the present invention.

本実施の形態では、特定の時刻および領域における住居侵入などの犯罪(特定事象)の発生密度を予測する。ここで、領域とは、犯罪が発生する地図上の位置を示すものである。発生密度を予測する特定事象が、たとえばインターネット上での詐欺などの犯罪や迷惑電話の場合、領域はURLなどのインターネット上の位置や電話番号を示すものとなる。   In the present embodiment, the occurrence density of crimes (specific events) such as residence intrusion at a specific time and area is predicted. Here, the area indicates a position on the map where the crime occurs. When the specific event for predicting the occurrence density is, for example, a crime such as fraud on the Internet or a nuisance call, the area indicates a location on the Internet such as a URL or a telephone number.

事象発生予測装置は、予測式構築部10と、予測部30とを有している。   The event occurrence prediction apparatus includes a prediction formula construction unit 10 and a prediction unit 30.

予測式構築部10は、特定事象の発生時刻tおよび特定事象の発生領域を特定する領域特定変数xの関数ρ(t,x)で特定事象の発生密度が与えられ、かつ、関数ρ(t,x)は外部因子{f}と関数ρ(t,x)の写像F[ρ(t,x)+{f}]で与えられるとして、写像F[ρ(t,x)+{f}]を過去に発生した特定事象の履歴データから求めて、関数ρ(t,x)を発生時刻tおよび領域特定変数xの関数として表現する。外部因子{f}とは、天候条件、地理構造などの環境要因や、パトロール状況などの特定事象の履歴データ以外の特定事象の発生密度に影響を与える因子である。予測式構築部10は、たとえば履歴データ群記憶部12とカーネル関数構築部14とカーネル関数記憶部16と予測式構築部18と予測式記憶部20と履歴データ受信部22を有している。   The prediction formula construction unit 10 is given the occurrence density of the specific event by the function ρ (t, x) of the region specific variable x that specifies the generation time t of the specific event and the generation region of the specific event, and the function ρ (t , X) is given by the external factor {f} and the map F [ρ (t, x) + {f}] of the function ρ (t, x), and the map F [ρ (t, x) + {f} ] Is obtained from the history data of the specific event that occurred in the past, and the function ρ (t, x) is expressed as a function of the occurrence time t and the region specific variable x. The external factor {f} is a factor that affects the occurrence density of specific events other than environmental factors such as weather conditions and geographical structure, and historical data of specific events such as patrol status. The prediction formula construction unit 10 includes, for example, a history data group storage unit 12, a kernel function construction unit 14, a kernel function storage unit 16, a prediction formula construction unit 18, a prediction formula storage unit 20, and a history data reception unit 22.

予測部30は、予測式構築部10が構築した予測式、すなわち、関数g(t,x)とデータを使って、特定事象の発生密度を日時t、場所xの関数として算出する。   The prediction unit 30 uses the prediction formula constructed by the prediction formula construction unit 10, that is, the function g (t, x) and data, to calculate the occurrence density of the specific event as a function of the date and time t and the location x.

事象発生予測装置は、表示部40を備えていてもよい。表示部40は、予測部30が算出した特定事象の発生密度を認識可能な状態に表示する。表示部40は、たとえば地図上に、特定事象の発生密度の等高線を表示する。   The event occurrence prediction device may include the display unit 40. The display unit 40 displays the occurrence density of the specific event calculated by the prediction unit 30 in a recognizable state. The display unit 40 displays, for example, contour lines of the occurrence density of specific events on a map.

次に、本実施の形態における、特定事象の発生密度の算出方法について説明する。時刻t、領域xにおける特定事象の発生密度ρ(t,x)を考える。この発生密度ρ(t,x)は、それ自身の写像Fで表されると仮定する。つまり、
ρ(t,x)=F[ρ(t,x)+{f}] …(1)
であるとする。ここで{f}は季節、環境要因、パトロール状況などの外部因子が考えられる。
Next, a method for calculating the occurrence density of a specific event in the present embodiment will be described. Consider the occurrence density ρ (t, x) of a specific event at time t and region x. This generation density ρ (t, x) is assumed to be represented by its own map F. That means
ρ (t, x) = F [ρ (t, x) + {f}] (1)
Suppose that Here, {f} may be external factors such as season, environmental factors, and patrol status.

予測式構築部10は、過去の特定事象の発生履歴データ群を用いて、写像Fを構築するという問題を解く。ここで、特定事象の発生履歴データ群とは、複数の履歴データの集合である。履歴データとは、特定事象が発生した時刻tおよび発生した領域を特定する領域特定変数xの組である。履歴データに特定事象の種類を示すインデックスを含めてもよい。   The prediction formula construction unit 10 solves the problem of constructing the mapping F using the past history data group of specific events. Here, the specific event occurrence history data group is a set of a plurality of history data. The history data is a set of a time t when a specific event has occurred and a region specifying variable x that specifies the region in which the specific event has occurred. An index indicating the type of specific event may be included in the history data.

写像F[ρ(t,x)+{f}]は、たとえば、偏微分方程式の解と仮定することができる。   The map F [ρ (t, x) + {f}] can be assumed to be a solution of a partial differential equation, for example.

また、偏微分方程式の解は、たとえばカーネル関数(グリーン関数、応答関数、積分核と呼ばれることもある)を使って書くことができる。   The solution of the partial differential equation can be written using, for example, a kernel function (sometimes called a Green function, a response function, or an integral kernel).

発生密度ρ(t,x)は、たとえば、以下の方程式を満たすと仮定することができる。   The generation density ρ (t, x) can be assumed to satisfy the following equation, for example.

ρ(t,x)=∫dx’∫dt’g(t,t’,x,x’)(ρ(t’,x’)+{f})+∫dx’g(t,0,x,x’)ρ(t=0,x) …(2)     ρ (t, x) = ∫dx′∫dt′g (t, t ′, x, x ′) (ρ (t ′, x ′) + {f}) + ∫dx′g (t, 0, x , X ′) ρ (t = 0, x) (2)

なお、右辺第一項のカーネル関数にかかるρ(t,x)に関する関数形や外部因子依存性と、第二項の形式には、幾つかのバリエーションが可能である。ここでは、最も単純な一例を挙げて議論を進めることにする。   It should be noted that several variations are possible for the function form and external factor dependency related to ρ (t, x) applied to the kernel function of the first term on the right side and the form of the second term. Here, the discussion will be made with the simplest example.

また、カーネル関数g(t,t’,x,x’)は過去の事象と現在の時空間との時間差、距離差の関数で書けると仮定できる場合がある。たとえば、発展方程式の形は、
ρ(t,x)=∫dx’∫dt’g(t−t’,x−x’)ρ(t’,x’)+∫dx’g(t,x−x’)ρ(t=0,x) …(2)
などが考えられる。
Further, it may be assumed that the kernel function g (t, t ′, x, x ′) can be written as a function of a time difference and a distance difference between a past event and the current time space. For example, the form of the evolution equation is
ρ (t, x) = ∫dx′∫dt′g (t−t ′, xx ′) ρ (t ′, x ′) + ∫dx′g (t, xx ′) ρ (t = 0, x) (2)
And so on.

この発展方程式を領域xについてフーリエ変換し、時刻tについてラプラス変換すると、以下の式が得られる。   When this evolution equation is Fourier transformed for region x and Laplace transformed for time t, the following equation is obtained.

g(z,k)=φ(z,k)/(φ(z,k)+1) …(3)
φ(t,k)=<ρ(t+t,k)/ρ(t,k)>t0 …(4)
g (z, k) = φ (z, k) / (φ (z, k) +1) (3)
φ (t, k) = <ρ (t + t 0 , k) / ρ (t 0 , k)> t0 (4)

過去の犯罪履歴データ、すなわち、ある時刻t、領域kにおける発生密度ρ(t,k)を用いて、φ(t,k)が求められる。ここで、過去の時刻t、領域xにおける発生密度ρ(t,x)は、正の整数値をとる。i番目に発生した特定事象について、発生時刻をt、発生領域xとすると、同時刻、同領域に他の事象が発生しなかった場合、ρ(t,x)=1となる。φ(t,k)が得られたら、そのφを用いてg(z,k)が得られる。g(t,x)が求められたら、予測犯罪密度λ(t,x)は、Self−Exciting Point Process(SEPP)モデルの時空間カーネル項に代入することで与えられる。Φ (t, k) is obtained using past crime history data, that is, the occurrence density ρ (t, k) at a certain time t and region k. Here, the generation density ρ (t, x) in the past time t and region x takes a positive integer value. Assuming that the occurrence time is t i and the occurrence region x i for the i-th specific event, ρ (t i , x i ) = 1 when no other event occurs at the same time and the same region. . When φ (t, k) is obtained, g (z, k) is obtained using the φ. Once g (t, x) is determined, the predicted crime density λ (t, x) is given by substituting it into the space-time kernel term of the Self-Exciting Point Process (SEPP) model.

λ(t,x)の計算方法の一例を以下に挙げる。
λ(t,x)=Σti<tg(t−t,x−x)+ρ(x) …(5)
ρ(x)=∫ dtρ(t,x)/T …(6)
ここで、(5)の和は時刻tより前に発生した特定事象についてのすべての和である。また、ρはバックグラウンドの寄与項であり、Tは元データの時間領域である。なお、ここでは、バックグラウンドの寄与項ρは上式で求めたが、これに限定されるものではない。
An example of a method for calculating λ (t, x) is given below.
λ (t, x) = Σti <t g (t−t i , x−x i ) + ρ b (x) (5)
ρ b (x) = ∫ 0 T dtρ (t, x) / T (6)
Here, the sum of (5) is the sum of all the specific events that occurred before time t. Ρ b is a background contribution term, and T is the time domain of the original data. Here, the background contribution term ρ b is obtained by the above equation, but is not limited thereto.

図2は、本実施の形態の事象発生予測装置を用いた事象発生予測方法のフローチャートである。   FIG. 2 is a flowchart of an event occurrence prediction method using the event occurrence prediction apparatus of the present embodiment.

この事象発生予測方法は、予測式構築フェーズと、予測フェーズに分けられる。予測式構築フェーズでは、特定事象の発生時刻tおよび特定事象の発生領域を特定する領域特定変数xの関数ρ(t,x)で特定事象の発生密度が与えられ、かつ、関数ρ(t,x)は前記外部因子{f}と前記関数ρ(t,x)の写像F[ρ(t,x)+{f}]で与えられるとして、前記写像F[ρ(t,x)+{f}]を過去に発生した前記特定事象の履歴データから求めて、関数ρ(t,x)を前記発生時刻tおよび前記領域特定変数xの関数として表現する。予測フェーズでは、写像F[ρ(t,x)+{f}]と過去データを使って、前記特定事象の発生密度を予測する。より、具体的には次のとおりである。   This event occurrence prediction method is divided into a prediction formula construction phase and a prediction phase. In the prediction formula construction phase, the occurrence density of the specific event is given by the function ρ (t, x) of the region specific variable x that specifies the generation time t of the specific event and the generation region of the specific event, and the function ρ (t, x) is given by the external factor {f} and the map F [ρ (t, x) + {f}] of the function ρ (t, x), and the map F [ρ (t, x) + { f}] is obtained from the history data of the specific event that occurred in the past, and the function ρ (t, x) is expressed as a function of the generation time t and the region specific variable x. In the prediction phase, the occurrence density of the specific event is predicted using the map F [ρ (t, x) + {f}] and past data. More specifically, it is as follows.

予測式構築フェーズでは、まず、過去の履歴データを履歴データ群記憶部12に蓄積する(S11)。過去の履歴データは、複数であってもよい。   In the prediction formula construction phase, first, past history data is accumulated in the history data group storage unit 12 (S11). There may be a plurality of past history data.

次に、カーネル関数構築部14は、履歴データ群記憶部12に蓄積された過去の履歴データ群を用いて、式(3)および式(4)からカーネル関数g(z,k)を求める(S12)。カーネル関数構築部14が導出したカーネル関数g(z,k)は、カーネル関数記憶部16に記憶される。このカーネル関数は、離散化した時刻z、および、離散化した領域を特定する領域特定変数kにおける値の表として記憶される。   Next, the kernel function construction unit 14 uses the past history data group stored in the history data group storage unit 12 to obtain the kernel function g (z, k) from the equations (3) and (4) ( S12). The kernel function g (z, k) derived by the kernel function construction unit 14 is stored in the kernel function storage unit 16. This kernel function is stored as a table of values at the discretized time z and the region specifying variable k that specifies the discretized region.

カーネル関数g(z,k)が得られたら、ラプラス逆変換、フーリエ逆変換を適用すると、予測式構築部18は、式(5)および式(6)から、特定事象の発生密度を予測する予測犯罪密度λ(t,x)を与える関数を構築する(S13)。予測式構築部18が構築した予測犯罪密度λ(t,x)は、予測式記憶部20に記憶される。この予測犯罪密度λは、離散化した時刻t、および、離散化した領域を特定する領域特定変数における値の表として記憶される。   When the kernel function g (z, k) is obtained, the Laplace inverse transform and the Fourier inverse transform are applied, and the prediction formula construction unit 18 predicts the occurrence density of the specific event from Formula (5) and Formula (6). A function that gives the predicted crime density λ (t, x) is constructed (S13). The predicted crime density λ (t, x) constructed by the prediction formula construction unit 18 is stored in the prediction formula storage unit 20. The predicted crime density λ is stored as a table of values in the discretized time t and the region specifying variable that specifies the discretized region.

予測犯罪密度λの関数が構築されたら、履歴データ受信部22は、新たな履歴データが入力されるのを監視する(S14)。たとえば一般ユーザがスマートフォンなどの携帯端末24を通じて犯罪が発生したという事実を入力すると、その履歴データが履歴データ受信部22に受信される。携帯端末24には、一般ユーザが、犯罪を目撃した際などに、新たな履歴データを入力し、その履歴データを履歴データ受信部22に対して送信するアプリケーションがインストールされているものとする。あるいは、警察などの情報提供機関が有する警察システム26から履歴データ受信部22に履歴データが伝達されてもよい。   When the function of the predicted crime density λ is constructed, the history data receiving unit 22 monitors the input of new history data (S14). For example, when a general user inputs the fact that a crime has occurred through a portable terminal 24 such as a smartphone, the history data is received by the history data receiving unit 22. It is assumed that an application for inputting new history data and transmitting the history data to the history data receiving unit 22 is installed in the mobile terminal 24 when a general user witnesses a crime. Alternatively, history data may be transmitted to the history data receiving unit 22 from a police system 26 possessed by an information provider such as the police.

履歴データ受信部22は、新たな履歴データが入力されるたびに繰り返し、監視を続ける。履歴データ受信部22に新たな履歴データが入力されたら、その履歴データは履歴データ群記憶部12に記憶され、工程S11に戻る。その結果、再び、工程S11〜工程S13が行われ、新たな予測犯罪密度λの関数が構築されて、予測式記憶部20に記憶される。   The history data receiving unit 22 repeats monitoring every time new history data is input. When new history data is input to the history data receiving unit 22, the history data is stored in the history data group storage unit 12, and the process returns to step S11. As a result, steps S11 to S13 are performed again, and a new function of the predicted crime density λ is constructed and stored in the prediction formula storage unit 20.

予測フェーズでは、まず、予測犯罪密度を予測する時刻および領域を設定する(S21)。予測する時刻は、特定の時刻でもよいし、幅をもっていてもよい。予測時間として、たとえば現在から所定の期間とする。予測する領域としては、特定の位置でもよいし、広がりをもっていてもよい。予測する領域として、たとえば予測式の構築のために履歴データを収集している領域全体とする。予測する時刻および領域の設定は、予測部30が行う。   In the prediction phase, first, the time and area for predicting the predicted crime density are set (S21). The predicted time may be a specific time or may have a width. The predicted time is, for example, a predetermined period from the present time. The area to be predicted may be a specific position or may have a spread. The region to be predicted is, for example, the entire region in which history data is collected for the construction of a prediction formula. The prediction unit 30 sets the time and area to be predicted.

次に、工程S21で設定した時刻および領域において、予測部30は予測式記憶部20から予測式を受信して、予測犯罪密度λを計算する(S22)。予測する時刻および領域が幅あるいは広がりを持っている場合には、離散化された時刻および領域について、それぞれ予測犯罪密度λの計算を繰り返す。   Next, at the time and region set in step S21, the prediction unit 30 receives the prediction formula from the prediction formula storage unit 20 and calculates the predicted crime density λ (S22). If the predicted time and area have a width or spread, the calculation of the predicted crime density λ is repeated for each discretized time and area.

計算された予測犯罪密度λは、表示部40において表示される(S23)。人間が認識可能なように表示された予測犯罪密度λは、警察官などによって参照され、パトロール活動の参考とされる。あるいは、携帯端末24に表示された予測犯罪密度λは、携帯端末24のユーザが犯罪を避ける行動に用いられる。   The calculated predicted crime density λ is displayed on the display unit 40 (S23). The predicted crime density λ displayed so as to be recognizable by human beings is referred to by a police officer or the like and used as a reference for patrol activities. Alternatively, the predicted crime density λ displayed on the mobile terminal 24 is used for the action of the user of the mobile terminal 24 avoiding the crime.

過去に発生した犯罪は、時間的・空間的に離れた位置において、次の犯罪を誘発する。時間差および距離が大きくなるほど、誘発する影響度は小さくなる傾向にある。たとえば、カリフォルニア州ロサンゼルスにおける不法侵入事件では、犯罪が発生すると、その1〜2日後と7日後に、1m圏内での再発生密度が高くなる。過去の犯罪の履歴データそれぞれが、将来のある時刻・ある位置における犯罪を誘発すると考えると、将来の犯罪の発生密度ρは、過去の犯罪の発生密度ρの写像Fとして与えられることになる。そこで、本実施の形態では、ρ(t,x)=F[ρ(t,x)+{f}]と仮定し、犯罪の発生密度の密度場に対するカーネル関数を履歴データから構築している。   Crimes that have occurred in the past induce the following crimes at positions that are separated in time and space. The greater the time difference and distance, the smaller the degree of influence induced. For example, in a trespassing incident in Los Angeles, California, when a crime occurs, the recurrence density within a 1 meter range increases 1-2 days and 7 days later. Assuming that each past crime history data induces a crime at a certain future time and a certain position, the future crime occurrence density ρ is given as a map F of the past crime occurrence density ρ. Therefore, in this embodiment, it is assumed that ρ (t, x) = F [ρ (t, x) + {f}], and a kernel function for the density field of crime occurrence density is constructed from history data. .

図3は、本実施の形態の事象発生予測装置を用いた犯罪発生率予測におけるカーネル関数gの時刻t依存性の評価例を示す図である。図4は、本実施の形態の事象発生予測装置を用いた犯罪発生密度予測におけるカーネル関数gの領域x依存性の評価例を示す図である。図3の横軸は時刻tの間の時間、縦軸はカーネル関数gの値である。図4の横軸は領域xの間の距離、縦軸はカーネル関数gの値である。   FIG. 3 is a diagram illustrating an evaluation example of the dependency of the kernel function g on the time t in the crime occurrence rate prediction using the event occurrence prediction apparatus of the present embodiment. FIG. 4 is a diagram illustrating an evaluation example of the dependency of the kernel function g on the region x in the crime occurrence density prediction using the event occurrence prediction apparatus of the present embodiment. The horizontal axis in FIG. 3 is the time between time t, and the vertical axis is the value of the kernel function g. The horizontal axis in FIG. 4 is the distance between the regions x, and the vertical axis is the value of the kernel function g.

図3に示すように、カーネル関数gは時間が大きくなるにしたがって、減少する傾向にある。しかし、図3に示すように、カーネル関数gは時間に対して単調減少ではない。また、図4に示すように、カーネル関数gは領域xの間の距離が大きくなるにしたがって、急激に小さくなっている。   As shown in FIG. 3, the kernel function g tends to decrease as the time increases. However, as shown in FIG. 3, the kernel function g is not monotonically decreasing with respect to time. Also, as shown in FIG. 4, the kernel function g decreases rapidly as the distance between the regions x increases.

図5は、本実施の形態の事象発生予測装置を用いた犯罪発生予測結果を示す等高線である。図5は、ある都市の738件の犯罪の履歴データを用いて計算した予測犯罪密度の等高線である。予測犯罪密度の計算には予測したい日の前日までのデータを使った。予測したい日のデータは含まれていない。また、図5には、予測した期間(1日)において、実際に発生した犯罪イベントを×記号で併せて示した。   FIG. 5 is a contour line showing a crime occurrence prediction result using the event occurrence prediction apparatus of the present embodiment. FIG. 5 is a contour line of predicted crime density calculated using history data of 738 crimes in a certain city. The data up to the day before the day of the prediction was used to calculate the predicted crime density. It does not include data for the day you want to predict. Further, in FIG. 5, crime events that actually occurred during the predicted period (one day) are also shown by crosses.

図5に示す通り、予測犯罪密度が高い位置において、実際に犯罪が発生している場合が多いことがわかる。   As shown in FIG. 5, it can be seen that there are many cases where crime actually occurs at a position where the predicted crime density is high.

本実施の形態の事象発生予測装置を用いた犯罪発生予測結果の精度は、過去の犯罪発生事象の犯罪の距離Δx(例えば空間分解能の半分の長さでスケールする)および時間Δt(例えば7日間でスケールする)における発生密度に対する影響g(Δt,Δx)をg(t,Δx)=1/((1+Δt)(1+Δx))と表現する方法(Prospective Hotspot Method)による予測結果の精度と比較して、概ね30%の精度向上となる。精度は、実際の犯罪のうち予測できた犯罪の件数を実際の犯罪件数で除したもので比較し、予測においては、予測対象領域を1km四方のセルに分割し、そのうちの10%を犯罪危険地域に指定した。   The accuracy of the crime occurrence prediction result using the event occurrence prediction device of the present embodiment is the crime distance Δx (for example, scaled by half the spatial resolution) and time Δt (for example, 7 days) of the past crime occurrence event. Compared with the accuracy of the prediction result by the method (Proactive Hotspot Method) of expressing g (Δt, Δx) as g (t, Δx) = 1 / ((1 + Δt) (1 + Δx)). Thus, the accuracy is improved by approximately 30%. The accuracy is compared by dividing the number of crimes that could be predicted among actual crimes by dividing the actual number of crimes, and in the prediction, the prediction target area was divided into 1 km square cells, and 10% of them were criminal risk Designated to the area.

このように本実施の形態によれば、予測精度は、Prospective Hotspot Methodよりも向上していることがわかる。   Thus, according to the present embodiment, it can be seen that the prediction accuracy is improved as compared to the prospective hotspot method.

なお、上述の実施の形態は、特定事象として犯罪を例として説明したが、ある1件の事象が新たに続く次の事象を誘発してカスケードする現象一般に適用可能である。特に、この誘発する割合が時間tおよび空間xの関数として見たときに時空間相関があり、かつその法則性にある程度の定常性があると予測される場合に、より高い精度を示すと期待される。このような事象として、たとえば地震の余震現象や空爆の被害状況などについては、過去にカスケード現象との関連が議論されている。さらには、自殺、疫病のパンデミック、ソーシャルネットワーキングサービス(SNS)上での勧誘・宗教・ねずみ講などの拡散、迷惑電話、服装などの流行、株価・金融商品の動向、選挙における投票、免疫系の異常事象、消費者の購買行動、ウェブサイト上の広告のクリック、ネットワーク上でのショッピングなどの需要予測、お見合いなどのマッチングにおける適合、サイバーテロを含むテロの予測、インフラの劣化予測、ソフトウェア・ハードウェアの故障や異常検知などについても、このようなカーネル関数を考えられると期待される。   In addition, although the above-mentioned embodiment demonstrated the crime as an example as a specific event, it is applicable to the phenomenon in which a certain event induces and cascades the next event which continues newly. In particular, this triggering rate is expected to show higher accuracy when there is a spatiotemporal correlation when viewed as a function of time t and space x, and the lawness is expected to have some degree of stationarity. Is done. As such events, for example, the aftershock phenomenon of an earthquake and the damage situation of an aerial bomb have been discussed in relation to the cascade phenomenon in the past. Furthermore, suicide, plague pandemic, social networking service (SNS) solicitation / religion / mouse lectures, nuisance calls, fashion trends, stock price / financial product trends, voting in the election, immune system Abnormal events, consumer purchasing behavior, clicks on advertisements on websites, demand forecasts such as shopping on the network, matching in matching, matching terrorism including cyber terrorism, infrastructure deterioration prediction, software / It is expected that such a kernel function can be considered for hardware failure and abnormality detection.

上述の実施の形態は、一つの特定事象を対象としてその将来の発生密度を予測するものであるが、多変数が影響しあうカスケード現象の予測にも応用できる。たとえば金融データに着目すると、それぞれの国における不動産の売買、国債・社債などの債権や株式の取引、現実の物の生産・販売、あるいはサービスの提供などによる現金の収受などが互いに影響しあい、外国為替取引におけるレート(為替レート)が変化する。この変化は、取引の瞬間における各国の不動産の状態に影響を受けるというよりも、過去の取引状態などにカスケード的に影響を受けている。   The above-described embodiment is for predicting the future occurrence density of one specific event, but it can also be applied to the prediction of a cascade phenomenon in which multiple variables affect. Focusing on financial data, for example, the purchase and sale of real estate in each country, the transaction of bonds and stocks such as government bonds and corporate bonds, the production and sale of real things, or the receipt of cash from the provision of services, etc. affect each other. The rate in exchange transactions (exchange rate) changes. This change is not influenced by the state of real estate in each country at the moment of transaction, but rather by the past transaction state.

n種の特定事象についてそれぞれの状態を予測する方法について考える。nは自然数である。ここで、特定事象の状態とは、発生密度でもよいし、為替レートなどの値であってもよい。i番目の特定事象の状態変数が関数ρ_i(t)で与えられるとする。   Consider a method for predicting each state for n types of specific events. n is a natural number. Here, the state of the specific event may be an occurrence density or a value such as an exchange rate. Assume that the state variable of the i-th specific event is given by the function ρ_i (t).

n種の要素を持つ状態変数を与えるの関数ベクトルΡ(t)={ρ_1(t),…,ρ_i(t),…,ρ_n(t)}を考える。カスケード現象であれば、この状態変数ベクトルΡ(t)は、この関数ベクトルΡ(t)および外部因子ベクトル{f}の写像F[Ρ(t)+{f}]で与えられる。   Consider a function vector Ρ (t) = {ρ_1 (t),..., ρ_i (t),..., ρ_n (t)} that gives a state variable having n types of elements. In the case of a cascade phenomenon, this state variable vector Ρ (t) is given by this function vector Ρ (t) and a map F [Ρ (t) + {f}] of the external factor vector {f}.

上述の実施の形態と同様にして、予測式構築部は、写像F[Ρ(t)+{f}]を過去に発生した特定事象の状態値と外部因子の履歴データから求めることができる。関数ベクトルΡ(t)が得られたら、予測部は、この関数ベクトルΡ(t)に将来の時刻を特定する値を入力して、特定事象の状態値を予測する。   Similar to the above-described embodiment, the prediction formula construction unit can obtain the map F [Ρ (t) + {f}] from the state value of the specific event that occurred in the past and the history data of the external factor. When the function vector Ρ (t) is obtained, the prediction unit inputs a value specifying a future time to the function vector Ρ (t), and predicts the state value of the specific event.

上述の実施の形態は、関数ベクトルΡ(t)が時刻tおよび特定事象の発生領域を特定する領域特定変数xの関数ρ(t,x)であると限定した場合である。xを離散化した場合には、xが示す各領域の発生密度ρ(t,x)の全体を関数ベクトルΡ(t)とみなすことができる。   The above-described embodiment is a case where the function vector t (t) is limited to the function ρ (t, x) of the region specifying variable x that specifies the time t and the occurrence region of the specific event. When x is discretized, the entire generation density ρ (t, x) of each region indicated by x can be regarded as a function vector Ρ (t).

たとえば2つの特定事象の状態値を予測する場合は、以下の通りである。2つの特定事象の状態値を与える状態関数をそれぞれρ(t)、ρ(t)とする。また、それぞれの過去の状態がそれ自身および他の特定事象に与える影響を表現したグリーン関数としてgaa(t)、gab(t)、gba(t)、gbb(t)を導入する。このとき、発展方程式は、以下のように書ける。For example, when predicting the state values of two specific events, it is as follows. State functions that give state values of two specific events are denoted by ρ a (t) and ρ b (t), respectively. In addition, g aa (t), g ab (t), g ba (t), and g bb (t) are introduced as green functions expressing the influence of each past state on itself and other specific events. . At this time, the evolution equation can be written as follows.

ただし、外部因子の影響はないとした。 However, external factors were not affected.

このように、本実施の形態では、複数の事象の過去の状態が影響しあうカスケード現象の予測にも用いることができる。   Thus, in this embodiment, it can also be used to predict a cascade phenomenon in which past states of a plurality of events affect each other.

10…予測式構築部、12…履歴データ群記憶部、14…カーネル関数構築部、16…カーネル関数記憶部、18…予測式構築部、20…予測式記憶部、22…履歴データ受信部、24…携帯端末、26…警察システム、30…予測部、40…表示部
DESCRIPTION OF SYMBOLS 10 ... Prediction formula construction part, 12 ... History data group storage part, 14 ... Kernel function construction part, 16 ... Kernel function storage part, 18 ... Prediction formula construction part, 20 ... Prediction formula storage part, 22 ... History data reception part, 24 ... mobile terminal, 26 ... police system, 30 ... prediction unit, 40 ... display unit

Claims (11)

nを自然数としたときに過去に発生したn種の特定事象の状態を表す値を予測する事象発生装置において、
i番目の特定事象の時刻tの状態が関数ρ_i(t)で与えられ、かつ、前記関数ρ_i(t)の関数ベクトルΡ(t)={ρ_1(t),…,ρ_i(t),…,ρ_n(t)}は外部因子ベクトル{f}と前記関数ベクトルΡ(t)の写像F[Ρ(t)+{f}]で与えられるとして、前記写像F[Ρ(t)+{f}]を前記特定事象の状態を表す値と外部因子の履歴データから求める予測式構築部と、
前記関数ベクトルΡ(t)に将来の時刻を特定する値を入力して、前記特定事象の発生密度を予測する予測部と、
を有することを特徴とする事象予測装置。
In an event generator for predicting a value representing the state of n types of specific events that occurred in the past when n is a natural number,
The state of the i-th specific event at time t is given by a function ρ_i (t), and the function vector Ρ (t) = {ρ_1 (t),..., ρ_i (t),. , Ρ_n (t)} is given by the external factor vector {f} and the map F [Ρ (t) + {f}] of the function vector Ρ (t), and the map F [Ρ (t) + {f }] From the value representing the state of the specific event and the history data of the external factor,
A prediction unit that inputs a value specifying a future time into the function vector Ρ (t) and predicts the occurrence density of the specific event;
An event prediction apparatus characterized by comprising:
前記関数ベクトルΡ(t)は、前記時刻tおよび特定事象の発生領域を特定する領域特定変数xの関数ρ(t,x)であることを特徴とする請求項1に記載の事象予測装置。   The event prediction apparatus according to claim 1, wherein the function vector t (t) is a function ρ (t, x) of a region specifying variable x that specifies the time t and a generation region of a specific event. 前記予測式構築部は、前記関数ρ(t,x)の偏微分方程式の解として写像F[ρ(t,x)+{f}]を前記履歴データを用いて再構築する問題を解く導出部を備えることを特徴とする請求項2に記載の事象発生予測装置。   The prediction formula construction unit derives a solution for reconstructing a map F [ρ (t, x) + {f}] using the history data as a solution of the partial differential equation of the function ρ (t, x). The event occurrence prediction device according to claim 2, further comprising a unit. 前記予測式構築部は、前記履歴データを用いて前記関数ρ(t,x)のカーネル関数を導出する導出部を備え、
前記関数ρ(t,x)は、ρ(t,x)=∫dx’∫dt’g(t,t’,x,x’)M(ρ(t’,x’),{f},ρb)+N(g(t,t’,x,x’),ρ(t,x),{f},ρb)で表現される発展方程式を満足する、
但し、M(ρ(t’,x’),{f},ρb)、及びN(g(t,t’,x,x’),ρ(t,x),{f},ρb)はρ(t,x)に関する線形項であり、ρbはρ(t,x)の時間平均によって決められる量である、
ことを特徴とする請求項3に記載の事象発生予測装置。
The prediction formula construction unit includes a derivation unit that derives a kernel function of the function ρ (t, x) using the history data,
The function ρ (t, x) is expressed as ρ (t, x) = ∫dx′∫dt′g (t, t ′, x, x ′) M (ρ (t ′, x ′), {f}, ρb) + N (g (t, t ′, x, x ′), ρ (t, x), {f}, ρb) is satisfied,
Where M (ρ (t ′, x ′), {f}, ρb) and N (g (t, t ′, x, x ′), ρ (t, x), {f}, ρb) are is a linear term for ρ (t, x), and ρb is a quantity determined by the time average of ρ (t, x).
The event occurrence prediction device according to claim 3.
過去の事象と現在の時空間との時間差t−t’、位置の差x−x’の関数であるカーネル関数を用いて前記関数ρ(t,x)はρ(t,x)=∫dx’∫dt’g(t−t’,x−x’)M(ρ(t’,x’),{f},ρb)+N(g(t,x),ρ(t,x),{f},ρb)で表現される発展方程式を満足するとする、
但し、M(ρ(t’,x’),{f},ρb)、及びN(g(t,x),ρ(t,x),{f},ρb)はρ(t,x)に関する線形項であり、ρbはρ(t,x)の時間平均によって決められる量である、
ことを特徴とする請求項4に記載の事象発生予測装置。
The function ρ (t, x) is expressed as ρ (t, x) = ∫dx using a kernel function that is a function of a time difference tt ′ between a past event and the current space-time and a position difference xx ′. '∫dt'g (t−t ′, xx ′) M (ρ (t ′, x ′), {f}, ρb) + N (g (t, x), ρ (t, x), { f}, ρb) satisfy the evolution equation expressed by
However, M (ρ (t ′, x ′), {f}, ρb) and N (g (t, x), ρ (t, x), {f}, ρb) are ρ (t, x). Ρb is a quantity determined by the time average of ρ (t, x),
The event occurrence prediction device according to claim 4, wherein:
前記カーネル関数g(t,x)は、特定事象が発生した前記時刻tおよび発生した領域を特定する前記領域特定変数xの組である履歴データを前記発展方程式の解に複数入力してサンプル平均を測定することで定められ、
前記予測部は、前記時刻tおよび前記領域特定変数xにおける事象発生率λ(x,t)を、Self−Exciting Point Process(SEPP)モデルの時空間カーネル項に代入することで求める、
ことを特徴とする請求項5に記載の事象発生予測装置。
The kernel function g (t, x) is a sample average obtained by inputting a plurality of history data, which is a set of the region specifying variable x specifying the time t at which a specific event has occurred and the region in which the specific event has occurred, to the solution of the evolution equation. Is determined by measuring
The prediction unit obtains the event occurrence rate λ (x, t) at the time t and the region specifying variable x by substituting it into a space-time kernel term of a Self-Exciting Point Process (SEPP) model.
The event occurrence prediction device according to claim 5, wherein:
新たな前記履歴データを収集して新たな前記履歴データが得られたときに前記予測式構築部に前記写像F[Ρ(t)+{f}]を再度求めさせるサーバーをさらに有することを特徴とする請求項1ないし請求項6のいずれか1項に記載の事象発生予測装置。   The server further includes a server that collects the new history data and obtains the map F [Ρ (t) + {f}] again when the new history data is obtained. The event occurrence prediction device according to any one of claims 1 to 6. 1つの事象が次に引き続く新たな事象に影響を与えるカスケード現象一般に対し、過去の特定事象の状態を表す値の履歴データに基づいて前記特定事象の将来の状態を表す値を予測する事象予測方法において、
時刻の関数Ρで前記特定事象の状態を表す値が与えられ、かつ、前記関数Ρは前記外部因子{f}と前記関数Ρの写像F[Ρ+{f}]で与えられるとして、前記写像F[Ρ+{f}]を過去の前記特定事象の状態を表す値の履歴データから求めて、前記関数Ρを前記時刻の関数として表現する予測式構築ステップと、
前記関数Ρに将来の時刻および領域を特定する値を入力して、前記特定事象の状態を表す値を予測する予測ステップと、
を有することを特徴とする事象予測方法。
An event prediction method for predicting a value representing a future state of a specific event based on historical data of a value representing a state of a specific event in the past, in general for a cascade phenomenon in which one event affects a subsequent new event In
A value representing the state of the specific event is given by a function Ρ of time, and the function F is given by the external factor {f} and a map F [Ρ + {f}] of the function 、. Obtaining [Ρ + {f}] from historical data of values representing the state of the specific event in the past, and expressing the function と し て as a function of the time;
A prediction step of inputting a value specifying a future time and region to the function Ρ and predicting a value representing the state of the specific event;
An event prediction method characterized by comprising:
1つの事象が次に引き続く新たな事象に影響を与えるカスケード現象一般に対し、過去の特定事象の状態を表す値の履歴データに基づいて前記特定事象の将来の状態を表す値を予測する事象予測プログラムにおいて、計算機に
時刻の関数Ρで前記特定事象の状態を表す値が与えられ、かつ、前記関数Ρは前記外部因子{f}と前記関数Ρの写像F[Ρ+{f}]で与えられるとして、前記写像F[Ρ+{f}]を過去の前記特定事象の状態を表す値の履歴データから求めて、前記関数Ρを前記時刻の関数として表現させる予測式構築ステップと、
前記関数Ρに将来の時刻および領域を特定する値を入力して、前記特定事象の状態を表す値を予測させる予測ステップと、
を実行させることを特徴とする事象予測プログラム。
An event prediction program for predicting a value representing a future state of a specific event based on historical data of a value representing a state of a specific event in the past, in general for a cascade phenomenon in which one event affects a subsequent new event , The computer is given a value representing the state of the specific event by a function Ρ of time, and the function Ρ is given by the external factor {f} and the map F [Ρ + {f}] of the function Ρ Obtaining the map F [Ρ + {f}] from historical data of values representing the state of the specific event in the past and expressing the function と し て as a function of the time;
A prediction step of inputting a value specifying a future time and region to the function 、 and predicting a value representing the state of the specific event;
An event prediction program characterized by causing
1つの事象が次に引き続く新たな事象に影響を与えるカスケード現象一般に対し、過去の特定事象の状態を表す値の履歴データに基づいて前記特定事象の将来の状態を表す値を予測する事象予測システムにおいて、
時刻の関数Ρで前記特定事象の状態を表す値が与えられ、かつ、前記関数Ρは前記外部因子{f}と前記関数Ρの写像F[Ρ+{f}]で与えられるとして、前記写像F[Ρ+{f}]を過去の前記特定事象の状態を表す値の履歴データから求めて、前記関数Ρを前記時刻の関数として表現させる予測式構築部と、
前記関数Ρに将来の時刻および領域を特定する値を入力して、前記特定事象の状態を表す値を予測する予測部と、
前記特定事象の状態に関する値を前記履歴データとして送信する携帯端末及びそれに準じる通信機能を持った端末と、
前記携帯端末から新たな前記履歴データを受信してして新たな前記履歴データが得られたときに前記予測式構築部に前記写像F[Ρ+{f}]を再度求めさせるサーバーと、
を有することを特徴とする事象予測システム。
An event prediction system for predicting a value representing a future state of a specific event based on historical data of a value representing a state of a specific event in the past, in general for a cascade phenomenon in which one event affects a subsequent new event In
A value representing the state of the specific event is given by a function Ρ of time, and the function F is given by the external factor {f} and a map F [Ρ + {f}] of the function 、. [Ρ + {f}] is obtained from history data of a value representing the state of the specific event in the past, and a prediction formula construction unit that expresses the function と し て as a function of the time;
A value that specifies a future time and region in the function 領域, and a prediction unit that predicts a value representing the state of the specific event;
A mobile terminal that transmits a value related to the state of the specific event as the history data, and a terminal having a communication function equivalent thereto;
A server that receives the new history data from the mobile terminal and obtains the map F [Ρ + {f}] again when the prediction formula construction unit obtains the new history data;
An event prediction system comprising:
nを自然数としたときに過去に発生したn種の特定事象の状態を表す値を予測する事象発生装置であって、i番目の特定事象の時刻tの状態が関数ρ_i(t)で与えられ、かつ、前記関数ρ_i(t)の関数ベクトルΡ(t)={ρ_1(t),…,ρ_i(t),…,ρ_n(t)}は外部因子ベクトル{f}と前記関数ベクトルΡ(t)の写像F[Ρ(t)+{f}]で与えられるとして、前記写像F[Ρ(t)+{f}]を前記特定事象の状態を表す値と外部因子の履歴データから求める予測式構築部と、前記関数ベクトルΡ(t)に将来の時刻を特定する値を入力して、前記特定事象の発生密度を予測する予測部と、を備えた事象予測装置とともに用いる携帯端末にインストールされるプログラムにおいて、携帯端末に、
前記事象予測装置に前記特定事象の状態を表す値を前記履歴データとして送信させる、
ことを特徴とする携帯端末にインストールされるプログラム。

An event generator that predicts a value representing the state of n types of specific events that occurred in the past when n is a natural number, and the state at time t of the i-th specific event is given by a function ρ_i (t) And the function vector の (t) = {ρ_1 (t),..., Ρ_i (t),..., Ρ_n (t)} of the function ρ_i (t) is an external factor vector {f} and the function vector Ρ ( t) Given the map F [Ρ (t) + {f}] of t), the map F [) (t) + {f}] is obtained from the value representing the state of the specific event and the history data of the external factor. A mobile terminal used together with an event prediction apparatus comprising: a prediction formula construction unit; and a prediction unit that inputs a value specifying a future time into the function vector Ρ (t) and predicts the occurrence density of the specific event In the program to be installed,
Causing the event prediction device to transmit a value representing the state of the specific event as the history data;
The program installed in the portable terminal characterized by this.

JP2018524169A 2016-06-22 2017-06-22 Event prediction device and event prediction method Active JP6935928B2 (en)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2016123259 2016-06-22
JP2016123259 2016-06-22
PCT/JP2017/023107 WO2017222030A1 (en) 2016-06-22 2017-06-22 Event prediction device and event prediction method

Publications (2)

Publication Number Publication Date
JPWO2017222030A1 true JPWO2017222030A1 (en) 2019-04-11
JP6935928B2 JP6935928B2 (en) 2021-09-15

Family

ID=60784128

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2018524169A Active JP6935928B2 (en) 2016-06-22 2017-06-22 Event prediction device and event prediction method

Country Status (2)

Country Link
JP (1) JP6935928B2 (en)
WO (1) WO2017222030A1 (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109102113A (en) * 2018-07-27 2018-12-28 阿里巴巴集团控股有限公司 Event prediction method and device, electronic equipment
JP7151382B2 (en) * 2018-11-01 2022-10-12 日本電信電話株式会社 Event prediction device, event prediction method and event prediction program
CN112005254A (en) * 2019-03-27 2020-11-27 株式会社奇异摄动 Event prediction device and event prediction method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04372046A (en) * 1991-06-20 1992-12-25 Hitachi Ltd Method and device for predicting demand amount
JPH07151369A (en) * 1993-11-30 1995-06-13 Toshiba Corp Heat load predicting apparatus and plant heat load predicting apparatus
JP2015069309A (en) * 2013-09-27 2015-04-13 Kddi株式会社 Information terminal device and scheduling method and program therefore

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04372046A (en) * 1991-06-20 1992-12-25 Hitachi Ltd Method and device for predicting demand amount
JPH07151369A (en) * 1993-11-30 1995-06-13 Toshiba Corp Heat load predicting apparatus and plant heat load predicting apparatus
JP2015069309A (en) * 2013-09-27 2015-04-13 Kddi株式会社 Information terminal device and scheduling method and program therefore

Also Published As

Publication number Publication date
JP6935928B2 (en) 2021-09-15
WO2017222030A1 (en) 2017-12-28

Similar Documents

Publication Publication Date Title
US10679260B2 (en) Cross-device message touchpoint attribution
US10755202B1 (en) Integrated risk analysis management
US8880097B1 (en) Systems and methods for statistically associating mobile devices to households
US20100114954A1 (en) Realtime popularity prediction for events and queries
US9208512B1 (en) Generating content for promotional messages based on distance
US9414199B2 (en) Predictive delivery of information based on device history
WO2017222030A1 (en) Event prediction device and event prediction method
CN111008335B (en) Information processing method, device, equipment and storage medium
US20230143484A1 (en) System, Method, and Computer Program Product for Analyzing Multivariate Time Series Using a Convolutional Fourier Network
US20170330231A1 (en) Method and system to display targeted ads based on ranking output of transactions
US20150309962A1 (en) Method and apparatus for modeling a population to predict individual behavior using location data from social network messages
US20160117737A1 (en) Preference Mapping for Automated Attribute-Selection in Campaign Design
US20140279053A1 (en) System and method for applying spatially indexed data to digital advertising bids
Mokhsin et al. Probability model: Malaysian consumer online shopping behavior towards online shopping scam
WO2020194642A1 (en) Event prediction device and event prediction method
US9306958B2 (en) Methods, systems and media for detecting non-intended traffic using co-visitation information
US20130203443A1 (en) Providing information about a location to a mobile device based on the location of the mobile device
Desai Mobile cloud computing in business
WO2015184199A9 (en) Systems and methods for generating contextual insight on electronic devices
US20140156388A1 (en) Tracking deal provider performance
Davydenko et al. Adaptation of cluster analysis methods in respect to vector space of social network analysis indicators for revealing suspicious government contracts
KR20200127531A (en) Method and appratus for providing online advertisement service
US20220051273A1 (en) Telecommunications Data Used For Lookalike Analysis
US20220358347A1 (en) Computerized system and method for distilled deep prediction for personalized stream ranking
Shapiro et al. The Employment Effects of Advances in Internet and Wireless Technology: Evaluating the Transitions from 2G to 3G and from 3G to 4G

Legal Events

Date Code Title Description
A711 Notification of change in applicant

Free format text: JAPANESE INTERMEDIATE CODE: A711

Effective date: 20200408

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A821

Effective date: 20200408

A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20200620

A871 Explanation of circumstances concerning accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A871

Effective date: 20200625

A975 Report on accelerated examination

Free format text: JAPANESE INTERMEDIATE CODE: A971005

Effective date: 20200722

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20200929

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20201130

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20210302

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20210506

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: 20210720

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20210819

R150 Certificate of patent or registration of utility model

Ref document number: 6935928

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250