JP6961312B2 - State change detection auxiliary device, state change detection device, state change detection auxiliary program, and state change detection program - Google Patents
State change detection auxiliary device, state change detection device, state change detection auxiliary program, and state change detection program Download PDFInfo
- Publication number
- JP6961312B2 JP6961312B2 JP2019160751A JP2019160751A JP6961312B2 JP 6961312 B2 JP6961312 B2 JP 6961312B2 JP 2019160751 A JP2019160751 A JP 2019160751A JP 2019160751 A JP2019160751 A JP 2019160751A JP 6961312 B2 JP6961312 B2 JP 6961312B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- state change
- value
- machine learning
- change detection
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims description 88
- 238000010801 machine learning Methods 0.000 claims description 75
- 238000004458 analytical method Methods 0.000 claims description 68
- 238000000034 method Methods 0.000 claims description 35
- 239000011159 matrix material Substances 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 14
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000006870 function Effects 0.000 claims description 8
- 238000010586 diagram Methods 0.000 description 10
- 238000007637 random forest analysis Methods 0.000 description 10
- 230000002159 abnormal effect Effects 0.000 description 9
- 238000004519 manufacturing process Methods 0.000 description 9
- 238000003860 storage Methods 0.000 description 8
- 230000005856 abnormality Effects 0.000 description 7
- 238000005259 measurement Methods 0.000 description 7
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000003066 decision tree Methods 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 238000000611 regression analysis Methods 0.000 description 3
- 238000013075 data extraction Methods 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 230000002547 anomalous effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 239000013598 vector Substances 0.000 description 1
Images
Landscapes
- Testing And Monitoring For Control Systems (AREA)
- Complex Calculations (AREA)
Description
この発明は、状態変動検出補助装置、状態変動検出装置、状態変動検出補助用プログラム、及び状態変動検出用プログラムに関するものである。 The present invention relates to a state change detection auxiliary device, a state change detection device, a state change detection auxiliary program, and a state change detection assisting program.
例えば、機械の電流と電圧と振動とを測定して、異常等の状態変動を検出する場合に、一般的には、人による過去の経験等に基づき異常の判定を行う。即ち、状態変動の閾値は人が設定するものであり、異常データを持っていない機械が状態変動の検出を行うことはほぼ不可能であった。 For example, when a state change such as an abnormality is detected by measuring the current, voltage, and vibration of a machine, the abnormality is generally determined based on the past experience of a person or the like. That is, the threshold value of the state change is set by a person, and it is almost impossible for a machine having no abnormal data to detect the state change.
特許文献1には、電波の異常伝搬が発生した場合の当該電波の到来状況を予測するシステムが開示されている。このシステムは、気象情報を取り込む気象情報取込手段と、この気象情報取込手段で取り込んだ気象情報から、高度毎の気温に関する情報、高度毎の気圧に関する情報、高度毎の湿度に関する情報の少なくとも1つを、気象データとして抽出する気象データ抽出手段を備える。更に、上記気象データ抽出手段で上記気象データが抽出された日にちの電波到来に関し、予め蓄積された測定データに基づき算出された上記電波の異常伝搬を予測するための予測関数を蓄積しており、蓄積されている予測関数に、上記気象データ抽出手段で抽出した気象データを入力して、上記電波の異常伝搬が発生する発生確率である予測値を決定するものである。
特許文献2には、生産物を生産可能な機械設備から、上記生産物の所定の生産単位ごとに上記機械設備の有するセンサの時系列データを取得する取得工程、前記時系列データを複数の分割データに分割した分割データ群を生成する分割工程、上記分割データ群に含まれる所定の正常品質分に対し、故障がないことを示す付加情報である第1付加情報を対応付けると共に、上記分割データ群に含まれる直近生産分に対し、故障があることを示す第2付加情報を無条件に対応付ける設定工程を備える品質異常検知方法が開示されている。
この品質異常検知方法は、更に、上記正常品質分及び上記直近生産分を混在させた2グループを形成するグループ化工程、上記2グループのうちの一方を用いた学習により、入力データに対応する上記付加情報を予測する予測モデルを生成する生成工程、上記2グループのうちの他方を上記入力データとすることによって得られる上記予測モデルの予測結果において、上記付加情報が上記第2付加情報であると予測される確率を算出する算出工程を有しており、上記確率に基づいて上記生産物の品質異常を判定すると言うものである。 This quality abnormality detection method further corresponds to the input data by a grouping step of forming two groups in which the normal quality component and the latest production component are mixed, and learning using one of the two groups. In the generation step of generating the prediction model for predicting the additional information, in the prediction result of the prediction model obtained by using the other of the two groups as the input data, the additional information is the second additional information. It has a calculation process for calculating the predicted probability, and determines the quality abnormality of the product based on the probability.
また、特許文献3には過去の実績に基づいて将来を予測するデータ予測方法が開示されている。このデータ予測方法では、予測モデルは、予測対象日における少なくとも一つの特徴量について予測値を出力する第1予測モデルと、この第1予測モデルから出力される予測値を入力因子に含み、予測対象日の所定時間ごとの予測値を出力する第2予測モデルと、から構成されるものである。
Further,
そして、特許文献3のデータ予測方法では、収集された至近実績データ及び過去実績データを用いて予測モデルを構築する予測モデル構築手段と、構築された予測モデルに予測用入力データを入力して予測を実行し、予測値を得る予測実行手段と、収集された至近実績データと予測値とから予測誤差またはモデル誤差を計算する予測誤差計算手段と、予測誤差またはモデル誤差に基づいて前記予測値を補正する補正係数又は補正量を算出し、軸方向補正予測値を得るものである。
Then, in the data prediction method of
図1(A)、(B)は、ランダムフォレストにより状態変動検出を行う場合の手法を示すための説明図である。一方、ランダムフォレストでは、図1(A)に示されるような領域の境界線によって、正常な場合を「1」により示し、異常な場合を「2」で示すようにした状態変動検出装置が知られている。即ち、この状態変動検出装置では縦軸方向の値が160以上で横軸方向の値が45より小さい領域は全て正常の「1」であり、縦軸の値が160より小さい領域では、横軸方向の値が30より小さければ正常の「1」であり、横軸方向の値が30以上であれば異常の「2」となる。また、縦軸方向の値が160以上で横軸方向の値が45より大きい領域は、縦軸方向の値が260より小さければ異常の「2」であり、縦軸方向の値が260以上であれば正常の「1」となる。図1(A)において、白無地の領域が正常の「1」であり、梨地模様の領域が異常の「2」である。 FIGS. 1 (A) and 1 (B) are explanatory views for showing a method for detecting state fluctuations in a random forest. On the other hand, in a random forest, a state change detection device is known in which a normal case is indicated by "1" and an abnormal case is indicated by "2" by the boundary line of the region as shown in FIG. 1 (A). Has been done. That is, in this state fluctuation detection device, the region where the value in the vertical axis direction is 160 or more and the value in the horizontal axis direction is smaller than 45 is all normal "1", and in the region where the value in the vertical axis is smaller than 160, the horizontal axis If the value in the direction is smaller than 30, it is a normal "1", and if the value in the horizontal axis direction is 30 or more, it is an abnormal "2". Further, the region where the value in the vertical axis direction is 160 or more and the value in the horizontal axis direction is larger than 45 is abnormal "2" when the value in the vertical axis direction is smaller than 260, and the value in the vertical axis direction is 260 or more. If there is, it becomes a normal "1". In FIG. 1 (A), the plain white region is a normal “1”, and the satin-patterned region is an abnormal “2”.
以上のような領域を分けて異常と正常の判定を行う場合には、分類器により(例えばツリー構造で)分岐を行って予測を行う機械学習による状態変動検出装置を実現できる。このような状態変動検出装置において、例えば、図1(B)に示すような解析対象データDが生じた場合には、図1(A)の領域に関する状態変動検出装置によると、縦軸方向の値が160以上で横軸方向の値が45より小さい領域に該当し、図示のように正常の「1」に分類されてしまう。しかしながら、縦軸方向の値が概ね10000であるから、通常の縦軸の値400などから矢印で示すように大きく離れており、明らかに異常の「2」に分類しなければならない解析対象データであり、誤分類が生じている。
When the above-mentioned areas are divided into abnormal and normal judgments, it is possible to realize a state change detection device by machine learning that makes predictions by branching (for example, in a tree structure) by a classifier. In such a state change detection device, for example, when the analysis target data D as shown in FIG. 1 (B) is generated, according to the state change detection device for the region of FIG. 1 (A), it is in the vertical direction. It corresponds to a region where the value is 160 or more and the value in the horizontal axis direction is smaller than 45, and is classified as a normal "1" as shown in the figure. However, since the value in the vertical axis direction is approximately 10,000, it is far from the normal
また、図2(A)、(B)は、ランダムフォレストにより状態変動検出を行う場合の手法の問題点を示すための説明図である。図2(A)と図2(B)に示すような領域が境界線によって分けられている場合において、目的変数が図2(A)に示す如くのダミー変数により規定されているとする。即ち、白無地の領域ではダミー変数が「1」であり、梨地模様の領域ではダミー変数が「99」であり、灰色無地の領域ではダミー変数が「100」であるときに、図2(B)に示すような解析対象データが得られたとする。 Further, FIGS. 2A and 2B are explanatory diagrams for showing the problems of the method when the state change detection is performed by the random forest. When the regions shown in FIGS. 2 (A) and 2 (B) are separated by a boundary line, it is assumed that the objective variable is defined by a dummy variable as shown in FIG. 2 (A). That is, when the dummy variable is "1" in the plain white area, the dummy variable is "99" in the satin pattern area, and the dummy variable is "100" in the plain gray area, FIG. 2 (B). ), It is assumed that the data to be analyzed is obtained.
図2(B)の場合の各領域の数値の平均と誤差は、図2(B)に記載の通りである。即ち、白無地の領域では数値平均が「43.4」で誤差が「42.4」であり、梨地模様の領域では数値平均が「99.4」で誤差が「0.4」であり、灰色無地の領域では数値平均が「100.0」で誤差が「0」である。このようにダミー変数に大小差があると、誤差が大きく異なってしまう。従って、誤差の比較によって状態変動を捕えようとしても適切な検出を行うことができないという問題がある。 The average and error of the numerical values in each region in the case of FIG. 2 (B) are as described in FIG. 2 (B). That is, in the plain white area, the numerical average is "43.4" and the error is "42.4", and in the pear-skin pattern area, the numerical average is "99.4" and the error is "0.4". In the plain gray area, the numerical average is "100.0" and the error is "0". If there is a difference in the size of the dummy variables in this way, the error will be significantly different. Therefore, there is a problem that appropriate detection cannot be performed even if an attempt is made to capture a state change by comparing errors.
特許文献4では、信号を取得する空間として二次元平面を考え、機械学習と識別に利用するデータは、平面内の分布情報とスペクトル情報との合計で三次元の構造をした多次元的情報と呼ばれるもので、この多次元的情報を用いた場合の、機械学習と識別の手順は、スペクトルデータを用いた場合と本質的には同様でとしている。更に、この場合には、データのパターンを記述するのに適切な特徴量を複数取得した上で、それを特徴ベクトルとし、機械学習及び識別処理に用いることもできること、特徴量の代表的例としては、体積や曲率、空間勾配、HLAC(高次局所自己相関)等があることが述べられている。
In
そして、機械学習に用いる特徴量を事前に選別することも可能であり、この場合、例えば、マハラノビス距離を算出して、識別に用いる特徴量を選別すれば良いことが述べられている。マハラノビス距離が大きければ、識別がより容易になるため、注目する各群間のマハラノビス距離が大きくなる様な特徴量を優先的に選別することもできるとしている。 Then, it is also possible to select the feature amount used for machine learning in advance. In this case, for example, it is stated that the Mahalanobis distance may be calculated and the feature amount used for identification may be selected. The larger the Mahalanobis distance, the easier it is to identify, so it is possible to preferentially select features that increase the Mahalanobis distance between each group of interest.
つまり、上記特許文献4に記載されているのは、機械学習及び識別処理をより容易にする特徴量の選定を行う場合にマハラノビス距離が大きくなる様な特徴量を選択することを述べている。
That is, the description in
更に、特許文献5の第0117欄には、「上記所定のアルゴリズムとして決定木を利用した場合には、ランダムフォレスト法等により、各特徴量(説明変数)の関連度(重要度又は寄与度とも称しても良い)を算出することができる。具体的には、決定木を作成する際には、ジニ係数、交差エントロピー等で示される不純度が小さくなるように、各特徴量(説明変数)からノード(分岐)の選択が行われる。そのため、各特徴量(説明変数)をノードに選択した際の不純度の減少量をその特徴量(説明変数)の重要度として利用することができる。上記ステップS102では、この重要度を利用して、複数種類の特徴量の中から異常発生の予測に有効な特徴量を選択するようにしてもよい。」と記載がなされており、基本的には図2を用いて説明したような重要度による選択が示されている。
Further, in column 0117 of
特許文献6には、計算されたマハラノビスの距離D2iの計算値が"1+3*σ"以内であれば、危険率0.3% で製造装置の状態に変化無し、と判断され、一方、"1+3*σ"以上であれば、危険率0.3% で製造装置状態に変化有り或いはその動作が異常の兆候有りと判断することが開示されている。このシステムは、D2INIの値を基準空間の標準偏差σより3σから3.5σの範囲とするものである。
According to
本実施形態は、解析対象データを取得する取得元の数が少ない場合にも、取得元の数が多いときに匹敵するような精度で状態変動を検出することを可能とする状態変動検出補助装置、状態変動検出装置、状態変動検出補助用プログラム、状態変動検出用プログラムを提供することを目的とする。 In this embodiment, even when the number of acquisition sources for acquiring the analysis target data is small, the state change detection assisting device can detect the state change with an accuracy comparable to that when the number of acquisition sources is large. , A state change detection device, a state change detection auxiliary program, and a state change detection program.
本実施形態に係る状態変動検出補助装置は、第1目的変数データと第1説明変数データを含む第1機械学習用データについて、前記第1機械学習用データ中の前記第1目的変数データとは異なるデータを第2目的変数データとし、前記第1機械学習用データ中の前記第1説明変数データとは異なるデータを第2説明変数データとした第2機械学習用データを作成する第2機械学習用データ作成手段と、説明変数から目的変数を予測する、教師データを用いて生成された第1予測モデルと、説明変数から目的変数を予測する、前記教師データから前記第2機械学習用データ作成手段によって作成された第2教師データを用いて生成された第2予測モデルと、複数の取得元から取得された第1目的変数データと第1説明変数データを含む第1解析対象データを前記第1予測モデルに適用し、第1予測値データを得る第1予測値取得手段と、前記第1解析対象データを前記第2機械学習用データ作成手段に適用して得られる第2解析対象データを前記第2予測モデルに適用し、第2予測値データを得る第2予測値取得手段と、前記第1解析対象データ中の前記第1目的変数データと前記第1予測値データとの差である第1予測誤差値と、前記第2解析対象データ中の前記第2目的変数データと前記第2予測値データとの差である第2予測誤差値と、を求める予測誤差値算出手段と、前記予測誤差値算出手段により算出された結果を出力手段へ出力する出力制御手段とを具備することを特徴とする。 The state change detection assisting device according to the present embodiment has the first objective variable data including the first objective variable data and the first explanatory variable data, and the first objective variable data in the first machine learning data. A second machine learning that creates a second machine learning data in which different data is used as the second objective variable data and data different from the first explanatory variable data in the first machine learning data is used as the second explanatory variable data. Data creation means, the objective variable is predicted from the explanatory variables, the first prediction model generated using the teacher data, the objective variable is predicted from the explanatory variables, and the second machine learning data creation from the teacher data. The first analysis target data including the second prediction model generated by using the second teacher data created by the means, the first objective variable data and the first explanatory variable data acquired from a plurality of acquisition sources is described above. The first predicted value acquisition means that applies to one prediction model and obtains the first predicted value data, and the second analysis target data obtained by applying the first analysis target data to the second machine learning data creation means. It is the difference between the second predicted value acquisition means that applies to the second predicted model and obtains the second predicted value data, and the first objective variable data and the first predicted value data in the first analysis target data. A prediction error value calculating means for obtaining a first prediction error value, a second prediction error value which is a difference between the second objective variable data and the second prediction value data in the second analysis target data, and the above. It is characterized by including an output control means for outputting the result calculated by the prediction error value calculation means to the output means.
以下添付図面を参照して本発明に係る状態変動検出補助装置1A、状態変動検出装置1B、状態変動検出補助用プログラム、及び状態変動検出用プログラムの実施形態を説明する。各図において同一の構成要素には、同一の符号を付して重複する説明を省略する。図3は、本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bを実現するコンピュータシステムの構成図である。本発明の実施形態に係る状態変動検出補助装置1A及び状態変動検出装置1Bは、例えば図3に示されるようなパーソナルコンピュータやワークステーション、その他のコンピュータシステムにより構成することができる。このコンピュータシステムは、CPU10が主メモリ11に記憶されている或いは主メモリ11に読み込んだプログラムやデータに基づき各部を制御し、必要な処理を実行することにより状態変動検出補助装置1Aや状態変動検出装置として動作を行うものである。
Hereinafter, embodiments of the state change detection
CPU10には、バス12を介して外部記憶インタフェース13、入力インタフェース14、表示インタフェース15、データ入力インタフェース16が接続されている。外部記憶インタフェース13には、状態変動検出用プログラム等のプログラムと必要なデータ等が記憶されている外部記憶装置23が接続されている。入力インタフェース14には、コマンドやデータを入力するための入力装置としてのキーボードなどの入力装置24とポインティングデバイスとしてのマウス22が接続されている。
An
表示インタフェース15には、LEDやLCDなどの表示画面を有する表示装置25が接続されている。データ入力インタフェース16には、測定データを得るためのセンサ26−1、26−2、・・・、26−mが接続されている。このコンピュータシステムには、他の構成が備えられていても良く、また、図3の構成は一例に過ぎない。
A
図4は、本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bの機能ブロック図である。上記において、CPU10では、外部記憶装置23内の状態変動検出補助用プログラムによって図4に記載の各手段等が実現される。即ち、予測モデル作成手段31、第2機械学習用データ作成手段32、第1予測値取得手段33、第2予測値取得手段34、予測誤差値算出手段35、出力制御手段36、誤差値伸長手段37、ダミー数値変更手段38が実現される。また、外部記憶装置23内には、教師データTが記憶されている。予測モデル作成手段31は、教師データTを用いて第1予測モデル30Aを作成するものである。また、予測モデル作成手段31は、本実施形態では、第1予測モデル30Aを作成するために予測モデル作成手段31を備えるが、他の装置やプログラムによって作成された第1予測モデル30Aを用いるものであっても良く、この場合には、予測モデル作成手段31を備えていなくともよい。
FIG. 4 is a functional block diagram of the state change detection
第1予測モデル30Aは、機械学習により説明変数から目的変数を予測するものである。ここに、機械学習のアルゴリズムとしては、パターンマイニングのランダムフォレストを挙げることができるが、これ以外に、分類木や回帰木などのように分類器により(例えばツリー構造で)分岐を行って予測を行う機械学習によるアリゴリズムを採用することができる。また、第1予測モデル30Aは、重回帰分析による機械学習を行うものであっても良い。
The
第2機械学習用データ作成手段32は、第1目的変数データと第1説明変数データを含む第1機械学習用データについて、上記第1機械学習用データ中の前記第1目的変数データとは異なるデータを第2目的変数データとし、前記第1機械学習用データ中の上記第1説明変数データとは異なるデータを第2説明変数データとした第2機械学習用データを作成するものである。 The second machine learning data creating means 32 is different from the first objective variable data in the first machine learning data with respect to the first machine learning data including the first objective variable data and the first explanatory variable data. The data is used as the second objective variable data, and the data different from the first explanatory variable data in the first machine learning data is used as the second explanatory variable data to create the second machine learning data.
また、第2機械学習用データ作成手段32は、与えられた第1機械学習用データが、第1目的変数データと第1説明変数データを含む教師データTである場合に、この教師データT中の上記第1目的変数データとは異なるデータを第2目的変数データとし、上記教師データT中の上記第1説明変数データとは異なるデータを第2説明変数データとする。そして第2機械学習用データ作成手段32は、上記第2説明変数データから上記第2目的変数データに対応する予測値が前記第2予測モデルにより導き出せる関係の上記第2目的変数データと上記第2説明変数データとを含む第2機械学習用データを作成するものである。この第2機械学習用データは新たな教師データである。従って、第2機械学習用データ作成手段32は、上記教師データTを第1教師データとして新たな第2教師データを作成するものと言える。 Further, the second machine learning data creating means 32 is in the teacher data T when the given first machine learning data is the teacher data T including the first objective variable data and the first explanatory variable data. The data different from the first objective variable data is referred to as the second objective variable data, and the data different from the first explanatory variable data in the teacher data T is referred to as the second explanatory variable data. Then, the second machine learning data creating means 32 has the relationship between the second objective variable data and the second objective variable data in which the predicted value corresponding to the second objective variable data can be derived from the second explanatory variable data by the second predicted model. The second machine learning data including the explanatory variable data is created. This second machine learning data is new teacher data. Therefore, it can be said that the second machine learning data creating means 32 creates new second teacher data using the teacher data T as the first teacher data.
更に、第2機械学習用データ作成手段32は、第1機械学習用データが複数の取得元から取得された第1目的変数データと第1説明変数データを含む第1解析対象データである場合に、上記第1解析対象データ中の上記第1目的変数データとは異なるデータを第2目的変数データとし、上記第1解析対象データ中の上記第1説明変数データとは異なるデータを第2説明変数データとして、新たな第2解析対象データを次のようにして作成する。即ち、第2機械学習用データ作成手段32は、上記第2説明変数データから上記第2目的変数データに対応する予測値が上記第2予測モデルにより導き出せる関係の上記第2目的変数データと上記第2説明変数データとを含む新たな第2解析対象データを作成するものである。 Further, the second machine learning data creating means 32 is when the first machine learning data is the first analysis target data including the first objective variable data and the first explanatory variable data acquired from a plurality of acquisition sources. The second objective variable data is data different from the first objective variable data in the first analysis target data, and the second explanatory variable is data different from the first explanatory variable data in the first analysis target data. As the data, new second analysis target data is created as follows. That is, the second machine learning data creating means 32 has the relationship between the second objective variable data and the second objective variable data in which the predicted value corresponding to the second objective variable data can be derived from the second explanatory variable data by the second predicted model. 2 A new second analysis target data including the explanatory variable data is created.
ここで、取得元とは、説明変数を幾つかのセンサを用いて取得・収集している場合に、その各センサを取得元と言う。また、説明変数を幾つかの装置を用いて取得・収集している場合に、その各装置を取得元と言う。つまり、目的変数や説明変数を取得する元となるところと言う意味で、取得元を用いる。 Here, the acquisition source is referred to as an acquisition source when the explanatory variables are acquired / collected using several sensors. Further, when the explanatory variables are acquired / collected using several devices, each device is referred to as an acquisition source. In other words, the acquisition source is used in the sense that it is the source for acquiring the objective variable and the explanatory variable.
上記第2機械学習用データ作成手段32によって、第1解析対象データから第2解析対象データを作成する具体例を説明する。例えば、「ある物」を製造する装置があり、時系列的にStep1〜Step4による4Stepで「ある物」を製造し完成するものとする。この製造装置にセンサ1、センサ2、センサ3、センサ4が設けられ、各Stepのときに検出したデータ値を取得する。センサ1〜4が取得するデータは、例えば、温度、湿度、電圧、電流、圧力など、「ある物」の生産時に製造装置から取得できる通常は、物理量である。勿論、センサを人の五感などとし、熱い、ぬるい、冷たいなどの感覚的な情報であっても良い。
A specific example of creating the second analysis target data from the first analysis target data by the second machine learning data creating means 32 will be described. For example, it is assumed that there is an apparatus for manufacturing "something", and "something" is manufactured and completed in 4 steps by
以上のようなセンサ1〜4による測定によって、図5(A)に示す如くの第1解析対象データであるデータが得られたものとする。図5(A)は、本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる第2機械学習用データ作成手段による処理を説明するための図である。この図5(A)では、取得された第1解析対象データを、時間経過を示す表題情報である「Step1〜Step4」を縦方向に複数並べ、取得元を示す取得元情報である「センサ1〜4」を横方向に複数並べて、各表題情報と各取得元情報の交差する位置に取得された値情報を配置した第1の行列としたものである。
It is assumed that the data which is the first analysis target data as shown in FIG. 5A is obtained by the measurement by the
予測モデル作成手段31は、上記第2機械学習用データ作成手段32によって教師データTから作成された第2教師データを用いて、第2予測モデル30Bを作成するものである。また、本実施形態では、第2予測モデル30Bを作成するために予測モデル作成手段31を備えるが、他の装置やプログラムによって作成された第2予測モデル30Bを用いるものであっても良く、この場合には、予測モデル作成手段31は不要である。
The prediction model creating means 31 creates the second prediction model 30B using the second teacher data created from the teacher data T by the second machine learning
第2予測モデル30Bは、機械学習により説明変数から目的変数を予測するものである。ここに、機械学習のアルゴリズムとしては、パターンマイニングのランダムフォレストを挙げることができるが、これ以外に、分類木や回帰木などのように分類器により(例えばツリー構造で)分岐を行って予測を行う機械学習によるアリゴリズムを採用することができる。また、第2予測モデル30Bは、重回帰分析による機械学習を行うものであっても良い。 The second prediction model 30B predicts the objective variable from the explanatory variables by machine learning. Here, as a machine learning algorithm, a random forest of pattern mining can be mentioned, but in addition to this, prediction is made by branching with a classifier (for example, in a tree structure) such as a classification tree or a regression tree. It is possible to adopt the algorithm by machine learning to be performed. Further, the second prediction model 30B may perform machine learning by multiple regression analysis.
上記第2機械学習用データ作成手段32は、上記第1の行列を転置して転置行列を作成し第2解析対象データとする。図5(B)は、本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる第2機械学習用データ作成手段による処理を説明するための図であって、上記転置行列を示す。上記第1の行列においては、「Step1〜Step4」が目的変数であるのに対し、転置行列では「センサ1〜4」が目的変数となっている。
The second machine learning data creating means 32 transposes the first matrix to create a transposed matrix and uses it as the second analysis target data. FIG. 5B is a diagram for explaining the processing by the second machine learning data creating means by the state change detection
図6(A)〜図6(C)は、本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる第2機械学習用データ作成手段による別の処理手法を説明するための図である。上記第2機械学習用データ作成手段32において、上記のように転置行列を作成する以外に、図6(A)の如き第1解析対象データが得られたとする。この場合に、第2機械学習用データ作成手段32は、「センサ1〜4」の各縦方向の列に係るデータを一つの塊として、例えば、最も左の縦方向1列のデータ「センサ1」のデータを目的変数とし、「センサ2〜4」の縦方向の各1列のデータを説明変数として、第2解析対象データを作成する。つまり、3列存在する説明変数の中から1列を取り出し、目的変数である1列のデータと位置を交換する。これにより、例えば、図6(B)の如き第2解析対象データを得る。
6 (A) to 6 (C) are for explaining another processing method by the second machine learning data creating means by the state change detection
更に、説明変数の位置を変えて図6(B)の如き第2解析対象データを得るようにしても良い。即ち、図6(A)は、センサ1の縦一列とセンサ2の縦一列を交代させたものであるのに対し、図6(C)は、図6(B)の状態においてセンサ3の一列とセンサ4の一列の位置を交代したものとなっている。このようにして得た第2解析対象データは、重回帰分析による機械学習を行う場合に好適である。
Further, the position of the explanatory variable may be changed to obtain the second analysis target data as shown in FIG. 6 (B). That is, while FIG. 6 (A) shows a vertical row of the
第1予測値取得手段33は、複数の取得元から取得された第1目的変数データと第1説明変数データを含む第1解析対象データを上記第1予測モデル30Aに適用し、第1予測値データを得るものである。第2予測値取得手段34は、上記第1解析対象データを上記第2機械学習用データ作成手段32に適用して得られる第2解析対象データを上記第2予測モデル30Bに適用し、第2予測値データを得るものである。予測誤差値算出手段35は、上記第1解析対象データ中の上記第1目的変数データと上記第1予測値データとの差である第1予測誤差値と、上記第2解析対象データ中の上記第2目的変数データと上記第2予測値データとの差である第2予測誤差値と、を求めるものである。
The first predicted value acquisition means 33 applies the first analysis target data including the first objective variable data and the first explanatory variable data acquired from a plurality of acquisition sources to the
出力制御手段36は、上記予測誤差値算出手段35により算出された結果を出力手段へ出力するものである。ここで出力手段は、既に説明した表示装置25とすることができる。また、表示装置25以外にプリンタなどであっても良い。以上のように実際に取得元から取得した第1解析対象データに基づき得られる第1予測誤差値だけでなく、第2機械学習用データ作成手段32が作成した第2解析対象データに基づき得られる第2予測誤差値を得ているので、取得元の数が多いときに匹敵するデータによって状態変動検出補助を行うことになり状態変動検出の確率が上がることが期待される。しかも、第2解析対象データは、第2説明変数データから前記第2目的変数データに対応する予測値が前記予測モデルにより導き出せる関係の前記第2目的変数データと前記第2説明変数データとを含むものであるから、上記期待は十分大きなものと言える。
The output control means 36 outputs the result calculated by the prediction error value calculation means 35 to the output means. Here, the output means can be the
誤差値伸長手段37は、上記第1予測モデル30Aの第1重要度を求めて、この第1重要度で上記第1予測誤差値に伸長を図る処理を行う共に、上記第2予測モデル30Bの第2重要度を求めて、この第2重要度で上記第2予測誤差値を伸長する処理を行うものである。重要度の求め方は、ランダムフォレストの予測モデルを作成する場合に重要度を求めて適切な決定木を作成するなどに用いる手法など、公知の手法を用いることができる。また、「伸長を図る処理」は、重要度の値を掛ける(乗算する)手法を挙げることができるが、重要度の値によっては除算を行うことも可能である。
The error value extending means 37 obtains the first importance of the
図7は、本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる誤差伸長を行った結果を示す図である。即ち図7(A)には、図5(A)に示された第1解析対象データから得られた上記第1誤差値に伸長を図る処理を行った結果を示し、7(B)には、図5(B)に示された第2解析対象データから得られた上記第2誤差値に伸長を図る処理を行った結果を示す。図7(A)における楕円DAにより示す部分には、他の部分よりも差異が大きく強調された誤差値が見られ、伸長の成果と考えることができる。また、図7(B)における楕円DBにより示す部分には、他の部分よりも差異が大きく強調された誤差値が見られ、伸長の成果と考えることができる。しかも、これは図5(B)に示された第2解析対象データから得られた上記第2誤差値に伸長を図る処理を行った結果であり、図7(A)と比較して差異が大きく、「第2予測誤差値を得ているので、取得元の数が多いときに匹敵するデータによって状態変動検出補助を行うことになり状態変動検出の確率が上がることが期待される。」と前述した効果が立証されている。
FIG. 7 is a diagram showing the results of error expansion by the state change detection
第1解析対象データにおいて、時間経過を示す表題情報が所定時間間隔で繰り返される情報であり、上記出力手段は、上記表題情報を第1の軸とし、上記取得元情報を第2の軸とし、時間を第3の軸として、3次元の出力を行うものとすることができる。しかも、上記3軸を移動させて多方面から見た如くに出力(表示)することができる。図8は、本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる3次元の出力を行った結果を示す図である。図8(A)は、横方向が時間軸、上下方向が先の例で「ステップ」とした表題情報の軸(第1の軸)、図の手前から奥へ延びる軸が取得元情報を割り当てた第2の軸である。この表示で、他の部分よりも差異が大きな部分(楕円で囲って示すP1)が識別されないが、図8(B)に示すように3軸を移動させることにより、他の部分よりも差異が大きな部分(楕円で囲って示すP1)が明瞭になり、オペレータが目視して状態変動検出を行う場合の補助に寄与することができる。
In the first analysis target data, the title information indicating the passage of time is information that is repeated at predetermined time intervals, and the output means has the title information as the first axis and the acquisition source information as the second axis. Three-dimensional output can be performed with time as the third axis. Moreover, the above three axes can be moved to output (display) as seen from various directions. FIG. 8 is a diagram showing the results of three-dimensional output by the state change detection
ダミー数値変更手段38は、上記目的変数の数値及び/または上記説明変数の数値をダミー数値へ変更し、上記目的変数の数値と上記説明変数の数値を近接させるものである。ダミー変数に大小差があると、誤差が大きく異なってしまうという問題については、図2を用いて説明した。ダミー数値変更手段38は、目的変数の数値及び/または上記説明変数の数値を比較して他と所定以上の差がある場合に、数値の近接を行う。例えば、既に説明した「ステップ」の数値が、1、13、21、22、31などと例えば、「1」以上の差がある場合には、10001、10002、10003、10004、10005と値を変更し、近接させる。これによって、ダミー変数の数値自体が大きくても、誤差値を得て比較を行う場合には、誤差が飛び抜けて大きな値を採ることが無くなり、状態変動検出を行う場合に好適となる。 The dummy numerical value changing means 38 changes the numerical value of the objective variable and / or the numerical value of the explanatory variable to a dummy numerical value, and brings the numerical value of the objective variable and the numerical value of the explanatory variable close to each other. The problem that the error is greatly different when there is a difference in the size of the dummy variable has been described with reference to FIG. The dummy numerical value changing means 38 compares the numerical values of the objective variable and / or the numerical values of the above explanatory variables, and when there is a difference of a predetermined value or more from the others, the numerical values are brought close to each other. For example, when the numerical value of the "step" already described has a difference of, for example, "1" or more from 1, 13, 21, 22, 31, etc., the value is changed to 10001, 10002, 1003, 10004, 10005. And bring them closer. As a result, even if the numerical value of the dummy variable itself is large, when the error value is obtained and the comparison is performed, the error does not jump out and a large value is not taken, which is suitable for the state fluctuation detection.
図4において、破線よりも上に記載の手段は、状態変動検出補助用プログラム及び状態変動検出用プログラムの双方で実現されるものであり、破線よりも下に記載の手段は、状態変動検出用プログラムのみで実現されるものである。CPU10では、外部記憶装置23内の状態変動検出用プログラムによって図4に記載の各手段等が実現される。状態変動検出用プログラムのみで実現されるものは、状態変動判定手段41である。
In FIG. 4, the means described above the broken line is realized by both the state change detection assisting program and the state change detection program, and the means described below the broken line is for state change detection. It is realized only by the program. In the
状態変動判定手段41は、上記状態変動検出補助装置1A上記出力制御手段から出力される結果について統計処理を行って、予め設定された判定条件情報に基づき状態変動の有無を判定するものである。「統計処理を行って、予め設定された判定条件情報に基づき状態変動の有無を判定する」ことの具体的一例として、上記状態変動判定手段41は、上記出力制御手段から出力される結果について単位ブレ距離を求め、求めた単位ブレ距離の平均値と偏差値を求め、偏差値に基づき状態変動の有無を判定する基礎情報を求めるという手法を挙げることができる。
The state change determination means 41 performs statistical processing on the result output from the state change detection
更に具体的には、第1解析対象データから得られた偏差値と第2解析対象データから得られた偏差値が共に所定値以上であることを状態変動の有無を判定する基礎情報とすることができる。そして、第1解析対象データの取得が前述の「ステップ」毎に得られたとき、1回分の基礎情報を得るものとし、基礎情報が複数回分得られるまでデータ取得を行う。このように複数の基礎情報を得た場合において、基礎情報が、「連続してN回状態変動有り」となったときに状態変動有りと最終的に判定するようにしても良い。 More specifically, the fact that both the deviation value obtained from the first analysis target data and the deviation value obtained from the second analysis target data are equal to or higher than a predetermined value is used as basic information for determining the presence or absence of state fluctuation. Can be done. Then, when the acquisition of the first analysis target data is obtained for each of the above-mentioned "steps", it is assumed that the basic information for one time is obtained, and the data is acquired until the basic information is obtained for a plurality of times. When a plurality of basic informations are obtained in this way, it may be finally determined that there is a state change when the basic information becomes "continuously N times state change".
以上のような手段等によって構成される状態変動検出補助装置1Aと状態変動検出装置1Bは、図9、図10に示すフローチャートによって処理動作を実行するので、このフローチャートを参照して動作説明を行う。図9のフローチャートは、第1予測モデル30Aと第2予測モデル30Bとの生成に係る処理及び第1重要度と第2重要度を求める処理を示すものであり、CPU10が行う処理である。与えられている教師データTを用いて第2機械学習用データ作成手段32による第2教師データ生成を行う(S11)。第2教師データ生成は、転置行列を求める処理によることができる。次に教師データTを用いて第1予測モデル30Aを生成する(S12)。このステップS12においては、ダミー数値変更手段38による目的変数の数値及び/または上記説明変数の数値をダミー数値へ変更を必要に応じて行うことができる。
Since the state change detection
更に、第2教師データ生成を用いて第2予測モデル30Bを生成する(S13)。このステップS13においても、ダミー数値変更手段38による目的変数の数値及び/または上記説明変数の数値をダミー数値へ変更を必要に応じて行うことができる。次に、教師データと第1予測モデル30Aを用いて第1重要度データを求める(S14)。また、第2教師データと第2予測モデル30Bを用いて第2重要度データを求める(S15)。上記ステップS14、S15の処理において、重要度の求め方は、ランダムフォレストの予測モデルを作成する場合に重要度を求めて適切な決定木を作成するなどに用いる手法など、公知の手法を用いることができる。
Further, the second prediction model 30B is generated by using the second teacher data generation (S13). Also in this step S13, the numerical value of the objective variable and / or the numerical value of the above explanatory variable by the dummy numerical value changing means 38 can be changed to a dummy numerical value as necessary. Next, the first importance data is obtained using the teacher data and the
図10のフローチャートは、状態変動検出に到るまでの処理を示すものであり、CPU10が行う処理である。第1目的変数データと第1説明変数データを含む第1解析対象データを複数の取得元から取得し、第1解析対象データから第2機械学習用データ作成手段32の処理によって第2解析対象データを求める(S21)。第1解析対象データを第1予測モデル30Aに適用して第1予測値を求め(S22)、第2解析対象データを第2予測モデル30Bに適用して第2予測値を求める(S23)。上記ステップS22、S23において、ダミー数値変更手段38による目的変数の数値及び/または上記説明変数の数値をダミー数値へ変更を必要に応じて行うことができる。
The flowchart of FIG. 10 shows the process up to the detection of the state change, and is the process performed by the
第1解析対象データの第1目的変数データと上記第1予測値との差により第1予測誤差値を求め、第2解析対象データの第2目的変数データと上記第2予測値との差により第2予測誤差値を求める(S24)。第1予測誤差値に第1重要度を掛けて第1伸長誤差値を求め、第2予測誤差値に第2重要度を掛けて第2伸長誤差値を求める(S25)。第1伸長誤差値から第1単位ブレ距離を算出し、第2伸長誤差値から第2単位ブレ距離を算出する(S26)。ここに、ブレ距離では、ユークリッド距離、マハラノビス距離などを用いることができる。 The first prediction error value is obtained from the difference between the first objective variable data of the first analysis target data and the first predicted value, and the difference between the second objective variable data of the second analysis target data and the second predicted value is used. The second prediction error value is obtained (S24). The first prediction error value is multiplied by the first importance to obtain the first elongation error value, and the second prediction error value is multiplied by the second importance to obtain the second elongation error value (S25). The first unit blur distance is calculated from the first elongation error value, and the second unit blur distance is calculated from the second elongation error value (S26). Here, the Euclidean distance, the Mahalanobis distance, and the like can be used as the blur distance.
第1単位ブレ距離の平均と偏差から第1偏差値を求め、第2単位ブレ距離の平均と偏差から第2偏差値を求める(S27)。第1単位ブレ距離に関する分布の標準偏差をσAとしたとき第1偏差値が例えば3σA以上であり、且つ、第2単位ブレ距離に関する分布の標準偏差をσBとしたとき第2偏差値が例えば3σB以上である測定対象データをbx=1、それ以外の測定対象データをbx=0とする(S28)。このbx=1、bx=0は、前述の基礎情報とすることができる。時系列で取得する測定対象データ(この測定対象データは、図5において、いずれか1つの「ステップ」といずれか1つの「センサ」によって特定されるデータである。)について基礎情報がbx=1であることが所定回連続したときに状態変動有りと判定する(S29)。 The first deviation value is obtained from the average and deviation of the first unit blur distance, and the second deviation value is obtained from the average and deviation of the second unit blur distance (S27). When the standard deviation of the distribution related to the first unit blur distance is σA, the first deviation value is, for example, 3σA or more, and when the standard deviation of the distribution related to the second unit blur distance is σB, the second deviation value is, for example, 3σB. The above measurement target data is set to bx = 1, and the other measurement target data is set to bx = 0 (S28). The bx = 1 and bx = 0 can be used as the above-mentioned basic information. The basic information about the measurement target data acquired in time series (this measurement target data is the data specified by any one "step" and any one "sensor" in FIG. 5) is bx = 1. It is determined that there is a state change when the above is consecutive a predetermined number of times (S29).
ステップS27以降の処理及び判定は一例に過ぎず、3σB以上は、σB以上、2σB以上、4σB以上、・・・とすることもできる。また、所定回連続しなくとも、発生頻度が所定割合以上に多くなることを条件とすることもできる。 The processing and determination after step S27 are only an example, and 3σB or more can be σB or more, 2σB or more, 4σB or more, and so on. Further, even if it does not continue a predetermined number of times, it can be a condition that the frequency of occurrence increases to a predetermined ratio or more.
以上のように本実施形態は、実際に取得元から取得した第1解析対象データに基づき得られる第1予測誤差値だけでなく、第2機械学習用データ作成手段32が作成した第2解析対象データに基づき得られる第2予測誤差値を得ているので、取得元の数が多いときに匹敵するデータによって状態変動検出補助を行うことになり状態変動検出の確率が上がることが期待されるものである。 As described above, in this embodiment, not only the first prediction error value obtained based on the first analysis target data actually acquired from the acquisition source, but also the second analysis target created by the second machine learning data creation means 32. Since the second prediction error value obtained based on the data is obtained, it is expected that the probability of state change detection will increase because the state change detection assistance will be performed by comparable data when the number of acquisition sources is large. Is.
本発明に係る複数の実施形態を説明したが、これらの実施形態は例として提示するものであり、発明の範囲を限定することは意図していない。これら新規な実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれるとともに、特許請求の範囲に記載された発明とその均等の範囲に含まれる。 Although a plurality of embodiments according to the present invention have been described, these embodiments are presented as examples and are not intended to limit the scope of the invention. These novel embodiments can be implemented in various other embodiments, and various omissions, replacements, and changes can be made without departing from the gist of the invention. These embodiments and modifications thereof are included in the scope and gist of the invention, and are also included in the scope of the invention described in the claims and the equivalent scope thereof.
1A…状態変動検出補助装置、1B…状態変動検出装置、10…CPU、11…主メモリ、12…バス、13…外部記憶インタフェース、14…入力インタフェース、15…表示インタフェース、16…データ入力インタフェース、22…マウス、23…外部記憶装置、24…入力装置、25…表示装置、26…センサ、30A、30B…予測モデル、31…予測モデル作成手段、32…第2機械学習用データ作成手段、33…第1予測値取得手段、34…第2予測値取得手段、35…予測誤差値算出手段、36…出力制御手段、37…誤差値伸長手段、38…ダミー数値変更手段、41…状態変動判定手段。 1A ... state change detection auxiliary device, 1B ... state change detection device, 10 ... CPU, 11 ... main memory, 12 ... bus, 13 ... external storage interface, 14 ... input interface, 15 ... display interface, 16 ... data input interface, 22 ... Mouse, 23 ... External storage device, 24 ... Input device, 25 ... Display device, 26 ... Sensor, 30A, 30B ... Prediction model, 31 ... Prediction model creation means, 32 ... Second machine learning data creation means, 33 ... 1st predicted value acquisition means, 34 ... 2nd predicted value acquisition means, 35 ... prediction error value calculation means, 36 ... output control means, 37 ... error value extension means, 38 ... dummy numerical value changing means, 41 ... state change determination means.
Claims (14)
説明変数から目的変数を予測する、教師データを用いて生成された第1予測モデルと、
説明変数から目的変数を予測する、前記教師データから前記第2機械学習用データ作成手段によって作成された第2教師データを用いて生成された第2予測モデルと、
複数の取得元から取得された第1目的変数データと第1説明変数データを含む第1解析対象データを前記第1予測モデルに適用し、第1予測値データを得る第1予測値取得手段と、
前記第1解析対象データを前記第2機械学習用データ作成手段に適用して得られる第2解析対象データを前記第2予測モデルに適用し、第2予測値データを得る第2予測値取得手段と、
前記第1解析対象データ中の前記第1目的変数データと前記第1予測値データとの差である第1予測誤差値と、前記第2解析対象データ中の前記第2目的変数データと前記第2予測値データとの差である第2予測誤差値と、を求める予測誤差値算出手段と、
前記予測誤差値算出手段により算出された結果を出力手段へ出力する出力制御手段と
を具備することを特徴とする状態変動検出補助装置。 Regarding the first machine learning data including the first objective variable data and the first explanatory variable data, data different from the first objective variable data in the first machine learning data is set as the second objective variable data, and the first 1. A second machine learning data creation means for creating a second machine learning data in which data different from the first explanatory variable data in the machine learning data is used as the second explanatory variable data.
The first prediction model generated using the teacher data, which predicts the objective variable from the explanatory variables,
A second prediction model that predicts the objective variable from the explanatory variables and is generated from the teacher data using the second teacher data created by the second machine learning data creation means.
As a first predicted value acquisition means for obtaining the first predicted value data by applying the first analysis target data including the first objective variable data and the first explanatory variable data acquired from a plurality of acquisition sources to the first predicted model. ,
A second predicted value acquisition means for obtaining second predicted value data by applying the second analysis target data obtained by applying the first analysis target data to the second machine learning data creating means to the second prediction model. When,
The first prediction error value, which is the difference between the first objective variable data and the first predicted value data in the first analysis target data, the second objective variable data in the second analysis target data, and the first. 2 Prediction error value calculation means for obtaining the second prediction error value, which is the difference from the prediction value data, and
A state change detection assisting device including an output control means for outputting the result calculated by the prediction error value calculation means to the output means.
前記出力制御手段は、前記誤差値伸長手段による誤差値伸長の結果を出力することを特徴とする請求項1または2に記載の状態変動検出補助装置。 The first importance of the first prediction model is obtained, and the first importance is used to extend the first prediction error value, and the second importance of the second prediction model is obtained. An error value extending means for performing a process of extending the second prediction error value with a second importance is provided.
The state change detection assisting device according to claim 1 or 2, wherein the output control means outputs the result of the error value extension by the error value extension means.
前記出力手段は、前記表題情報を第1の軸とし、前記取得元を示す取得元情報を第2の軸とし、時間を第3の軸として、3次元の出力を行うことを特徴とする請求項1乃至4のいずれか1項に記載の状態変動検出補助装置。 In the first analysis target data, the title information indicating the passage of time is information that is repeated at predetermined time intervals.
The claim is characterized in that the output means performs three-dimensional output with the title information as the first axis, the acquisition source information indicating the acquisition source as the second axis, and time as the third axis. Item 2. The state change detection assisting device according to any one of Items 1 to 4.
前記状態変動検出補助装置の前記出力制御手段から出力される結果について統計処理を行って、予め設定された判定条件情報に基づき状態変動の有無を判定する状態変動判定手段と、
を具備する特徴とする状態変動検出装置。 The state change detection assisting device according to any one of claims 1 to 4.
A state change determining means for determining the presence or absence of a state change based on preset determination condition information by performing statistical processing on the result output from the output control means of the state change detecting auxiliary device.
A state change detection device characterized by the present invention.
第1目的変数データと第1説明変数データを含む第1機械学習用データについて、前記第1機械学習用データ中の前記第1目的変数データとは異なるデータを第2目的変数データとし、前記第1機械学習用データ中の前記第1説明変数データとは異なるデータを第2説明変数データとした第2機械学習用データを作成する第2機械学習用データ作成手段、
説明変数から目的変数を予測する、教師データを用いて生成された第1予測モデル、
説明変数から目的変数を予測する、前記教師データから前記第2機械学習用データ作成手段によって作成された第2教師データを用いて生成された第2予測モデル、
複数の取得元から取得された第1目的変数データと第1説明変数データを含む第1解析対象データを前記第1予測モデルに適用し、第1予測値データを得る第1予測値取得手段、
前記第1解析対象データを前記第2機械学習用データ作成手段に適用して得られる第2解析対象データを前記第2予測モデルに適用し、第2予測値データを得る第2予測値取得手段、
前記第1解析対象データ中の前記第1目的変数データと前記第1予測値データとの差である第1予測誤差値と、前記第2解析対象データ中の前記第2目的変数データと前記第2予測値データとの差である第2予測誤差値と、を求める予測誤差値算出手段、
前記予測誤差値算出手段により算出された結果を出力手段へ出力する出力制御手段
として機能させることを特徴とする状態変動検出補助用プログラム。 Computer,
Regarding the first machine learning data including the first objective variable data and the first explanatory variable data, data different from the first objective variable data in the first machine learning data is set as the second objective variable data, and the first A second machine learning data creation means for creating a second machine learning data in which data different from the first explanatory variable data in the machine learning data is used as the second explanatory variable data.
First prediction model generated using teacher data, which predicts the objective variable from the explanatory variables,
A second prediction model that predicts the objective variable from the explanatory variables and is generated from the teacher data using the second teacher data created by the second machine learning data creation means.
A first predicted value acquisition means for obtaining first predicted value data by applying first analysis target data including first objective variable data and first explanatory variable data acquired from a plurality of acquisition sources to the first predicted model.
A second predicted value acquisition means for obtaining second predicted value data by applying the second analysis target data obtained by applying the first analysis target data to the second machine learning data creating means to the second prediction model. ,
The first prediction error value, which is the difference between the first objective variable data and the first predicted value data in the first analysis target data, the second objective variable data in the second analysis target data, and the first. 2 Prediction error value calculation means for obtaining the second prediction error value, which is the difference from the prediction value data.
A state fluctuation detection assisting program characterized in that it functions as an output control means for outputting the result calculated by the prediction error value calculation means to the output means.
前記第1予測モデルの第1重要度を求めて、この第1重要度で前記第1予測誤差値に伸長を図る処理を行うと共に、前記第2予測モデルの第2重要度を求めて、この第2重要度で前記第2予測誤差値を伸長する処理を行う誤差値伸長手段として機能させ、
前記コンピュータを前記出力制御手段として、前記誤差値伸長手段による誤差値伸長の結果を出力するように機能させることを特徴とする請求項8または9に記載の状態変動検出補助用プログラム。 Further, the computer
The first importance of the first prediction model is obtained, and the first importance is used to extend the first prediction error value, and the second importance of the second prediction model is obtained. It is made to function as an error value extension means for performing a process of extending the second prediction error value at the second importance.
The state change detection assisting program according to claim 8 or 9, wherein the computer functions as the output control means so as to output the result of the error value expansion by the error value expansion means.
第1目的変数データと第1説明変数データを含む第1機械学習用データについて、前記第1機械学習用データ中の前記第1目的変数データとは異なるデータを第2目的変数データとし、前記第1機械学習用データ中の前記第1説明変数データとは異なるデータを第2説明変数データとした第2機械学習用データを作成する第2機械学習用データ作成手段、
説明変数から目的変数を予測する、教師データを用いて生成された第1予測モデル、
説明変数から目的変数を予測する、前記教師データから前記第2機械学習用データ作成手段によって作成された第2教師データを用いて生成された第2予測モデル、
複数の取得元から取得された第1目的変数データと第1説明変数データを含む第1解析対象データを前記第1予測モデルに適用し、第1予測値データを得る第1予測値取得手段、
前記第1解析対象データを前記第2機械学習用データ作成手段に適用して得られる第2解析対象データを前記第2予測モデルに適用し、第2予測値データを得る第2予測値取得手段、
前記第1解析対象データ中の前記第1目的変数データと前記第1予測値データとの差である第1予測誤差値と、前記第2解析対象データ中の前記第2目的変数データと前記第2予測値データとの差である第2予測誤差値と、を求める予測誤差値算出手段、
前記予測誤差値算出手段により算出された結果を出力手段へ出力する出力制御手段、
前記出力制御手段から出力される結果について統計処理を行って、予め設定された判定条件情報に基づき状態変動の有無を判定する状態変動判定手段、
として機能させることを特徴とする状態変動検出用プログラム。 Computer,
Regarding the first machine learning data including the first objective variable data and the first explanatory variable data, data different from the first objective variable data in the first machine learning data is set as the second objective variable data, and the first A second machine learning data creation means for creating a second machine learning data in which data different from the first explanatory variable data in the machine learning data is used as the second explanatory variable data.
First prediction model generated using teacher data, which predicts the objective variable from the explanatory variables,
A second prediction model that predicts the objective variable from the explanatory variables and is generated from the teacher data using the second teacher data created by the second machine learning data creation means.
A first predicted value acquisition means for obtaining first predicted value data by applying first analysis target data including first objective variable data and first explanatory variable data acquired from a plurality of acquisition sources to the first predicted model.
A second predicted value acquisition means for obtaining second predicted value data by applying the second analysis target data obtained by applying the first analysis target data to the second machine learning data creating means to the second prediction model. ,
The first prediction error value, which is the difference between the first objective variable data and the first predicted value data in the first analysis target data, the second objective variable data in the second analysis target data, and the first. 2 Prediction error value calculation means for obtaining the second prediction error value, which is the difference from the prediction value data.
An output control means that outputs the result calculated by the prediction error value calculation means to the output means.
A state change determination means that performs statistical processing on the result output from the output control means and determines the presence or absence of a state change based on preset determination condition information.
A program for detecting state fluctuations, which is characterized by functioning as.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019160751A JP6961312B2 (en) | 2019-09-03 | 2019-09-03 | State change detection auxiliary device, state change detection device, state change detection auxiliary program, and state change detection program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019160751A JP6961312B2 (en) | 2019-09-03 | 2019-09-03 | State change detection auxiliary device, state change detection device, state change detection auxiliary program, and state change detection program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2021039565A JP2021039565A (en) | 2021-03-11 |
JP6961312B2 true JP6961312B2 (en) | 2021-11-05 |
Family
ID=74847151
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019160751A Active JP6961312B2 (en) | 2019-09-03 | 2019-09-03 | State change detection auxiliary device, state change detection device, state change detection auxiliary program, and state change detection program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6961312B2 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116724276B (en) * | 2021-03-18 | 2024-06-21 | 三菱电机株式会社 | Storage medium readable by computer system, operation improvement support device, and operation improvement support method |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010211684A (en) * | 2009-03-12 | 2010-09-24 | Toshiba Corp | Method, program and apparatus for processing data |
JP5674623B2 (en) * | 2011-05-31 | 2015-02-25 | 京セラドキュメントソリューションズ株式会社 | Image forming apparatus and image forming system |
EP3125057B1 (en) * | 2014-03-27 | 2019-07-03 | Nec Corporation | System-analyzing device, analysis-model generation method, system analysis method, and system-analyzing program |
JP6164311B1 (en) * | 2016-01-21 | 2017-07-19 | 日本電気株式会社 | Information processing apparatus, information processing method, and program |
JP6844204B2 (en) * | 2016-03-22 | 2021-03-17 | 日本製鉄株式会社 | Prediction table design equipment, methods and programs |
WO2017221856A1 (en) * | 2016-06-21 | 2017-12-28 | 日本電気株式会社 | Analysis device, analysis method, and recording medium |
JP2017228086A (en) * | 2016-06-22 | 2017-12-28 | 富士通株式会社 | Machine learning management program, machine learning management method, and machine learning management device |
WO2019111866A1 (en) * | 2017-12-07 | 2019-06-13 | 日本電気株式会社 | Analysis device, communication system, data processing method, and non-transitory computer readable medium |
-
2019
- 2019-09-03 JP JP2019160751A patent/JP6961312B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2021039565A (en) | 2021-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108628281B (en) | Abnormality detection system and abnormality detection method | |
EP2905665B1 (en) | Information processing apparatus, diagnosis method, and program | |
US11216741B2 (en) | Analysis apparatus, analysis method, and non-transitory computer readable medium | |
CN108780314B (en) | Model adaptation and online learning for unstable environments | |
US7882050B2 (en) | Data division apparatus, data division method and program | |
Kaya et al. | Process capability analyses based on fuzzy measurements and fuzzy control charts | |
JP5538597B2 (en) | Anomaly detection method and anomaly detection system | |
US11657121B2 (en) | Abnormality detection device, abnormality detection method and computer readable medium | |
WO2019026134A1 (en) | Information processing device and information processing method | |
WO2020188696A1 (en) | Abnormality detection device, and abnormality detection method | |
US10884841B2 (en) | Information processing device, information processing method, and recording medium | |
JP2019082943A (en) | Causal relationship model construction system and method | |
JP6961312B2 (en) | State change detection auxiliary device, state change detection device, state change detection auxiliary program, and state change detection program | |
de Moraes et al. | Data-driven failure prediction in brittle materials: A phase field-based machine learning framework | |
CN115617867A (en) | Time series prediction method, electronic device and storage medium | |
US20230229136A1 (en) | Abnormal irregularity cause identifying device, abnormal irregularity cause identifying method, and abnormal irregularity cause identifying program | |
JP7012928B2 (en) | State change detection device and state change detection program | |
JP6929260B2 (en) | Time-series feature extraction device, time-series feature extraction method and program | |
JP6910997B2 (en) | Information processing equipment, calculation method and calculation program | |
JP6765769B2 (en) | State change detection device and state change detection program | |
JP6233038B2 (en) | Assembly yield prediction apparatus, assembly yield prediction program, and assembly yield prediction method | |
CN111868715A (en) | Method for synchronizing multivariate time series data and multivariate time series data processing device | |
JP6347771B2 (en) | Abnormality diagnosis apparatus, abnormality diagnosis method, and abnormality diagnosis program | |
JP6115607B2 (en) | Abnormality diagnosis apparatus, abnormality diagnosis method, and abnormality diagnosis program | |
JP7233258B2 (en) | Abnormal diagnosis device |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200221 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210420 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20210713 |
|
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: 20211012 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20211012 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6961312 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313111 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |