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 PDF

Info

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
Application number
JP2019160751A
Other languages
Japanese (ja)
Other versions
JP2021039565A (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.)
Toshiba Information Systems Japan Corp
Original Assignee
Toshiba Information Systems Japan Corp
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 Toshiba Information Systems Japan Corp filed Critical Toshiba Information Systems Japan Corp
Priority to JP2019160751A priority Critical patent/JP6961312B2/en
Publication of JP2021039565A publication Critical patent/JP2021039565A/en
Application granted granted Critical
Publication of JP6961312B2 publication Critical patent/JP6961312B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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つを、気象データとして抽出する気象データ抽出手段を備える。更に、上記気象データ抽出手段で上記気象データが抽出された日にちの電波到来に関し、予め蓄積された測定データに基づき算出された上記電波の異常伝搬を予測するための予測関数を蓄積しており、蓄積されている予測関数に、上記気象データ抽出手段で抽出した気象データを入力して、上記電波の異常伝搬が発生する発生確率である予測値を決定するものである。 Patent Document 1 discloses a system for predicting the arrival status of the radio wave when an abnormal propagation of the radio wave occurs. This system uses the meteorological information importing means for capturing meteorological information and the meteorological information captured by this meteorological information capturing means to at least information on temperature at each altitude, pressure at each altitude, and humidity at each altitude. One is provided with a meteorological data extraction means for extracting the meteorological data. Furthermore, regarding the arrival of radio waves on the day when the weather data was extracted by the weather data extraction means, a prediction function for predicting abnormal propagation of the radio waves calculated based on the measurement data accumulated in advance is accumulated. Meteorological data extracted by the meteorological data extracting means is input to the accumulated prediction function to determine a predicted value which is the probability of occurrence of abnormal propagation of the radio wave.

特許文献2には、生産物を生産可能な機械設備から、上記生産物の所定の生産単位ごとに上記機械設備の有するセンサの時系列データを取得する取得工程、前記時系列データを複数の分割データに分割した分割データ群を生成する分割工程、上記分割データ群に含まれる所定の正常品質分に対し、故障がないことを示す付加情報である第1付加情報を対応付けると共に、上記分割データ群に含まれる直近生産分に対し、故障があることを示す第2付加情報を無条件に対応付ける設定工程を備える品質異常検知方法が開示されている。 Patent Document 2 describes an acquisition step of acquiring time-series data of a sensor possessed by the machine-equipment for each predetermined production unit of the product from a machine-equipment capable of producing the product, and a plurality of divisions of the time-series data. In the division step of generating the divided data group divided into data, the first additional information which is additional information indicating that there is no failure is associated with the predetermined normal quality portion included in the divided data group, and the divided data group is associated with the first additional information. A quality abnormality detection method including a setting process for unconditionally associating a second additional information indicating that there is a failure with respect to the latest production amount included in the above is disclosed.

この品質異常検知方法は、更に、上記正常品質分及び上記直近生産分を混在させた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, Patent Document 3 discloses a data prediction method for predicting the future based on past achievements. In this data prediction method, the prediction model includes a first prediction model that outputs a prediction value for at least one feature quantity on the prediction target date and a prediction value output from this first prediction model as input factors, and is a prediction target. It is composed of a second prediction model that outputs a predicted value for each predetermined time of the day.

そして、特許文献3のデータ予測方法では、収集された至近実績データ及び過去実績データを用いて予測モデルを構築する予測モデル構築手段と、構築された予測モデルに予測用入力データを入力して予測を実行し、予測値を得る予測実行手段と、収集された至近実績データと予測値とから予測誤差またはモデル誤差を計算する予測誤差計算手段と、予測誤差またはモデル誤差に基づいて前記予測値を補正する補正係数又は補正量を算出し、軸方向補正予測値を得るものである。 Then, in the data prediction method of Patent Document 3, a prediction model construction means for constructing a prediction model using the collected closest actual data and past actual data, and a prediction input data for input to the constructed prediction model for prediction. The prediction execution means for obtaining the prediction value, the prediction error calculation means for calculating the prediction error or the model error from the collected nearest actual data and the prediction value, and the prediction value based on the prediction error or the model error. The correction coefficient or the correction amount to be corrected is calculated, and the axial correction predicted value is obtained.

図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 vertical axis value 400 and the like as shown by the arrow, and it is clearly anomalous "2" for analysis target data. Yes, there is a misclassification.

また、図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 Patent Document 4, a two-dimensional plane is considered as a space for acquiring signals, and the data used for machine learning and identification is multidimensional information having a three-dimensional structure in total of distribution information and spectral information in the plane. The procedure for machine learning and identification when using this multidimensional information is essentially the same as when using spectral data. Further, in this case, it is possible to acquire a plurality of feature quantities suitable for describing a data pattern and then use them as feature vectors for machine learning and identification processing, as a typical example of feature quantities. Is stated to have volume, curvature, spatial gradient, HLAC (higher-order local autocorrelation), and the like.

そして、機械学習に用いる特徴量を事前に選別することも可能であり、この場合、例えば、マハラノビス距離を算出して、識別に用いる特徴量を選別すれば良いことが述べられている。マハラノビス距離が大きければ、識別がより容易になるため、注目する各群間のマハラノビス距離が大きくなる様な特徴量を優先的に選別することもできるとしている。 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 Patent Document 4 states that when selecting a feature amount that facilitates machine learning and identification processing, a feature amount that increases the Mahalanobis distance is selected.

更に、特許文献5の第0117欄には、「上記所定のアルゴリズムとして決定木を利用した場合には、ランダムフォレスト法等により、各特徴量(説明変数)の関連度(重要度又は寄与度とも称しても良い)を算出することができる。具体的には、決定木を作成する際には、ジニ係数、交差エントロピー等で示される不純度が小さくなるように、各特徴量(説明変数)からノード(分岐)の選択が行われる。そのため、各特徴量(説明変数)をノードに選択した際の不純度の減少量をその特徴量(説明変数)の重要度として利用することができる。上記ステップS102では、この重要度を利用して、複数種類の特徴量の中から異常発生の予測に有効な特徴量を選択するようにしてもよい。」と記載がなされており、基本的には図2を用いて説明したような重要度による選択が示されている。 Further, in column 0117 of Patent Document 5, "When a decision tree is used as the above-mentioned predetermined algorithm, the degree of relevance (importance or contribution) of each feature amount (explanatory variable) is determined by a random forest method or the like. (It may be called) can be calculated. Specifically, when creating a decision tree, each feature amount (explanatory variable) so as to reduce the impureness indicated by the Gini coefficient, cross-entropy, etc. Therefore, the node (branch) is selected from the above. Therefore, the amount of decrease in purity when each feature amount (explanatory variable) is selected as the node can be used as the importance of the feature amount (explanatory variable). In step S102, the importance may be used to select a feature amount effective for predicting the occurrence of an abnormality from a plurality of types of feature amounts. " Shows selection by importance as described with reference to FIG.

特許文献6には、計算されたマハラノビスの距離D2iの計算値が"1+3*σ"以内であれば、危険率0.3% で製造装置の状態に変化無し、と判断され、一方、"1+3*σ"以上であれば、危険率0.3% で製造装置状態に変化有り或いはその動作が異常の兆候有りと判断することが開示されている。このシステムは、D2INIの値を基準空間の標準偏差σより3σから3.5σの範囲とするものである。 According to Patent Document 6, if the calculated Mahalanobis distance D2i is within "1 + 3 * σ", it is determined that there is no change in the state of the manufacturing apparatus with a risk rate of 0.3%, while "1 + 3". If it is * σ "or more, it is disclosed that it is judged that there is a change in the state of the manufacturing equipment or there is a sign of abnormality in the operation at a risk rate of 0.3%. In this system, the value of D2INI is in the range of 3σ to 3.5σ from the standard deviation σ of the reference space.

特開2005−315753号公報Japanese Unexamined Patent Publication No. 2005-315753 特開2018−147406号公報Japanese Unexamined Patent Publication No. 2018-147406 特開2004−94437号公報Japanese Unexamined Patent Publication No. 2004-94437 特開2015−52581号公報Japanese Unexamined Patent Publication No. 2015-52581 特開2018−116545号公報JP-A-2018-116545 特開2000−252179号公報Japanese Unexamined Patent Publication No. 2000-252179

本実施形態は、解析対象データを取得する取得元の数が少ない場合にも、取得元の数が多いときに匹敵するような精度で状態変動を検出することを可能とする状態変動検出補助装置、状態変動検出装置、状態変動検出補助用プログラム、状態変動検出用プログラムを提供することを目的とする。 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.

ランダムフォレストにより状態変動検出を行う場合の手法を示すための説明図。Explanatory drawing for showing the method in the case of performing state change detection by a random forest. ランダムフォレストにより状態変動検出を行う場合の手法の問題点を示すための説明図。Explanatory drawing to show the problem of the method when the state change detection is performed by a random forest. 本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bを実現するコンピュータシステムの構成図。The block diagram of the computer system which realizes the state change detection auxiliary device 1A and the state change detection device 1B which concerns on embodiment of this invention. 本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bの機能ブロック図。The functional block diagram of the state change detection auxiliary device 1A and the state change detection device 1B which concerns on embodiment of this invention. 本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる第2機械学習用データ作成手段による処理を説明するための図。The figure for demonstrating the process by the 2nd machine learning data creation means by the state change detection auxiliary device 1A and the state change detection device 1B which concerns on embodiment of this invention. 本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる第2機械学習用データ作成手段による別の処理手法を説明するための図。The figure for demonstrating another processing method by the 2nd machine learning data creation means by the state change detection auxiliary device 1A and the state change detection device 1B which concerns on embodiment of this invention. 本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる誤差伸長を行った結果を示す図。The figure which shows the result of having performed the error extension by the state change detection auxiliary device 1A and the state change detection device 1B which concerns on embodiment of this invention. 本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる3次元の出力を行った結果を示す図。The figure which shows the result of having performed the three-dimensional output by the state change detection auxiliary device 1A and the state change detection device 1B which concerns on embodiment of this invention. 本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる第1予測モデルと第2予測モデルとの生成に係る処理及び第1重要度と第2重要度を求める処理を示すフローチャート。The process related to the generation of the first prediction model and the second prediction model by the state change detection auxiliary device 1A and the state change detection device 1B and the process of obtaining the first importance and the second importance according to the embodiment of the present invention are shown. flowchart. 本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bによる処理を示すフローチャート。The flowchart which shows the processing by the state change detection auxiliary apparatus 1A and the state change detection apparatus 1B which concerns on embodiment of this invention.

以下添付図面を参照して本発明に係る状態変動検出補助装置1A、状態変動検出装置1B、状態変動検出補助用プログラム、及び状態変動検出用プログラムの実施形態を説明する。各図において同一の構成要素には、同一の符号を付して重複する説明を省略する。図3は、本発明の実施形態に係る状態変動検出補助装置1Aと状態変動検出装置1Bを実現するコンピュータシステムの構成図である。本発明の実施形態に係る状態変動検出補助装置1A及び状態変動検出装置1Bは、例えば図3に示されるようなパーソナルコンピュータやワークステーション、その他のコンピュータシステムにより構成することができる。このコンピュータシステムは、CPU10が主メモリ11に記憶されている或いは主メモリ11に読み込んだプログラムやデータに基づき各部を制御し、必要な処理を実行することにより状態変動検出補助装置1Aや状態変動検出装置として動作を行うものである。 Hereinafter, embodiments of the state change detection auxiliary device 1A, the state change detection device 1B, the state change detection auxiliary program, and the state change detection program according to the present invention will be described with reference to the accompanying drawings. In each figure, the same components are designated by the same reference numerals, and duplicate description will be omitted. FIG. 3 is a configuration diagram of a computer system that realizes the state change detection auxiliary device 1A and the state change detection device 1B according to the embodiment of the present invention. The state change detection auxiliary device 1A and the state change detection device 1B according to the embodiment of the present invention can be configured by, for example, a personal computer, a workstation, or another computer system as shown in FIG. In this computer system, the CPU 10 controls each part based on the program or data stored in the main memory 11 or read into the main memory 11 and executes necessary processing to execute the state change detection auxiliary device 1A or the state change detection. It operates as a device.

CPU10には、バス12を介して外部記憶インタフェース13、入力インタフェース14、表示インタフェース15、データ入力インタフェース16が接続されている。外部記憶インタフェース13には、状態変動検出用プログラム等のプログラムと必要なデータ等が記憶されている外部記憶装置23が接続されている。入力インタフェース14には、コマンドやデータを入力するための入力装置としてのキーボードなどの入力装置24とポインティングデバイスとしてのマウス22が接続されている。 An external storage interface 13, an input interface 14, a display interface 15, and a data input interface 16 are connected to the CPU 10 via a bus 12. A program such as a state change detection program and an external storage device 23 in which necessary data and the like are stored are connected to the external storage interface 13. An input device 24 such as a keyboard as an input device for inputting commands and data and a mouse 22 as a pointing device are connected to the input interface 14.

表示インタフェース15には、LEDやLCDなどの表示画面を有する表示装置25が接続されている。データ入力インタフェース16には、測定データを得るためのセンサ26−1、26−2、・・・、26−mが接続されている。このコンピュータシステムには、他の構成が備えられていても良く、また、図3の構成は一例に過ぎない。 A display device 25 having a display screen such as an LED or an LCD is connected to the display interface 15. Sensors 26-1, 26-2, ..., 26-m for obtaining measurement data are connected to the data input interface 16. This computer system may be provided with other configurations, and the configuration of FIG. 3 is only an example.

図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 auxiliary device 1A and the state change detection device 1B according to the embodiment of the present invention. In the above, in the CPU 10, each means and the like shown in FIG. 4 are realized by the state change detection assisting program in the external storage device 23. That is, the prediction model creation means 31, the second machine learning data creation means 32, the first prediction value acquisition means 33, the second prediction value acquisition means 34, the prediction error value calculation means 35, the output control means 36, and the error value extension means. 37, the dummy numerical value changing means 38 is realized. Further, the teacher data T is stored in the external storage device 23. The prediction model creating means 31 creates the first prediction model 30A using the teacher data T. Further, in the present embodiment, the prediction model creating means 31 includes the prediction model creating means 31 for creating the first prediction model 30A, but uses the first prediction model 30A created by another device or program. In this case, the prediction model creating means 31 may not be provided.

第1予測モデル30Aは、機械学習により説明変数から目的変数を予測するものである。ここに、機械学習のアルゴリズムとしては、パターンマイニングのランダムフォレストを挙げることができるが、これ以外に、分類木や回帰木などのように分類器により(例えばツリー構造で)分岐を行って予測を行う機械学習によるアリゴリズムを採用することができる。また、第1予測モデル30Aは、重回帰分析による機械学習を行うものであっても良い。 The first prediction model 30A 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 first prediction model 30A may perform machine learning by multiple regression analysis.

第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 Steps 1 to 4 in chronological order. A sensor 1, a sensor 2, a sensor 3, and a sensor 4 are provided in this manufacturing apparatus, and the data values detected at each step are acquired. The data acquired by the sensors 1 to 4 is usually a physical quantity that can be acquired from the manufacturing apparatus at the time of production of "something" such as temperature, humidity, voltage, current, and pressure. Of course, the sensor may be the five senses of a person, and may be sensory information such as hot, lukewarm, or cold.

以上のようなセンサ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 sensors 1 to 4 as described above. FIG. 5A is a diagram for explaining processing by the second machine learning data creating means by the state change detection auxiliary device 1A and the state change detection device 1B according to the embodiment of the present invention. In FIG. 5A, a plurality of acquired first analysis target data “Step1 to Step4” which are title information indicating the passage of time are arranged in the vertical direction, and “Sensor 1” which is the acquisition source information indicating the acquisition source. ~ 4 ”are arranged in the horizontal direction to form a first matrix in which the acquired value information is arranged at the intersection of each title information and each acquisition source information.

予測モデル作成手段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 data creating means 32. Further, in the present embodiment, the prediction model creating means 31 is provided for creating the second prediction model 30B, but the second prediction model 30B created by another device or program may be used. In this case, the prediction model creating means 31 is unnecessary.

第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 auxiliary device 1A and the state change detection device 1B according to the embodiment of the present invention, and is a diagram for explaining the transposed matrix. Is shown. In the first matrix, "Step1 to Step4" are the objective variables, whereas in the transposed matrix, "sensors 1 to 4" are the objective variables.

図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 auxiliary device 1A and the state change detection device 1B according to the embodiment of the present invention. It is a figure of. It is assumed that the second machine learning data creating means 32 obtains the first analysis target data as shown in FIG. 6A in addition to creating the transposed matrix as described above. In this case, the second machine learning data creating means 32 considers the data related to each vertical column of "sensors 1 to 4" as one block, for example, the leftmost vertical column of data "sensor 1". The data of "" is used as the objective variable, and the data of each one column in the vertical direction of "sensors 2 to 4" is used as the explanatory variable to create the second analysis target data. That is, one column is taken out from the explanatory variables existing in three columns, and the position is exchanged with the data of one column which is the objective variable. As a result, for example, the second analysis target data as shown in FIG. 6B is obtained.

更に、説明変数の位置を変えて図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 sensor 1 and a vertical row of the sensor 2 alternated, FIG. 6 (C) shows a row of the sensor 3 in the state of FIG. 6 (B). And the position of one row of the sensor 4 are changed. The second analysis target data obtained in this way is suitable for performing machine learning by multiple regression analysis.

第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 first prediction model 30A, and the first predicted value. It is for getting data. The second prediction value acquisition means 34 applies the second analysis target data obtained by applying the first analysis target data to the second machine learning data creation means 32 to the second prediction model 30B, and second. Predicted value data is obtained. The prediction error value calculating means 35 has a first prediction error value which is a difference between the first objective variable data and the first predicted value data in the first analysis target data, and the above second analysis target data. The second prediction error value, which is the difference between the second objective variable data and the second predicted value data, is obtained.

出力制御手段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 display device 25 already described. In addition to the display device 25, a printer or the like may be used. As described above, 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 data created by the second machine learning data creation means 32 is obtained. Since the second prediction error value 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. Moreover, the second analysis target data includes the second objective variable data and the second explanatory 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 prediction model. It can be said that the above expectations are sufficiently high because they are unreasonable.

誤差値伸長手段37は、上記第1予測モデル30Aの第1重要度を求めて、この第1重要度で上記第1予測誤差値に伸長を図る処理を行う共に、上記第2予測モデル30Bの第2重要度を求めて、この第2重要度で上記第2予測誤差値を伸長する処理を行うものである。重要度の求め方は、ランダムフォレストの予測モデルを作成する場合に重要度を求めて適切な決定木を作成するなどに用いる手法など、公知の手法を用いることができる。また、「伸長を図る処理」は、重要度の値を掛ける(乗算する)手法を挙げることができるが、重要度の値によっては除算を行うことも可能である。 The error value extending means 37 obtains the first importance of the first prediction model 30A, performs a process of extending the first importance to the first prediction error value at the first importance, and extends the first prediction error value of the second prediction model 30B. The second importance is obtained, and the process of extending the second prediction error value at this second importance is performed. As the method of determining the importance, a known method such as a method used for obtaining an appropriate decision tree when creating a prediction model of a random forest can be used. Further, as the "processing for stretching", a method of multiplying (multiplying) the value of importance can be mentioned, but division can also be performed depending on the value of importance.

図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 auxiliary device 1A and the state change detection device 1B according to the embodiment of the present invention. That is, FIG. 7 (A) shows the result of performing the process of extending the first error value obtained from the first analysis target data shown in FIG. 5 (A), and FIG. 7 (B) shows the result. , The result of performing the process of extending the second error value obtained from the second analysis target data shown in FIG. 5B is shown. In the portion indicated by the ellipse DA in FIG. 7 (A), an error value in which the difference is greatly emphasized is seen as compared with the other portions, which can be considered as the result of elongation. Further, in the portion indicated by the ellipse DB in FIG. 7B, an error value in which the difference is greatly emphasized is seen as compared with the other portions, which can be considered as the result of elongation. Moreover, this is the result of performing a process of extending the second error value obtained from the second analysis target data shown in FIG. 5 (B), and there is a difference as compared with FIG. 7 (A). Largely, "Since the second prediction error value 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." The effects mentioned above have been proven.

第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 auxiliary device 1A and the state change detection device 1B according to the embodiment of the present invention. In FIG. 8A, the horizontal direction is the time axis, the vertical direction is the axis of the title information (first axis) designated as “step” in the previous example, and the axis extending from the front to the back of the figure is assigned the acquisition source information. It is the second axis. In this display, the part having a larger difference than the other parts (P1 shown by enclosing it in an ellipse) is not identified, but by moving the three axes as shown in FIG. 8 (B), the difference is larger than the other parts. The large portion (P1 shown by enclosing it in an ellipse) becomes clear, which can contribute to the assistance when the operator visually detects the state change.

ダミー数値変更手段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 CPU 10, each means and the like shown in FIG. 4 is realized by a state change detection program in the external storage device 23. What is realized only by the state change detection program is the state change determination means 41.

状態変動判定手段41は、上記状態変動検出補助装置1A上記出力制御手段から出力される結果について統計処理を行って、予め設定された判定条件情報に基づき状態変動の有無を判定するものである。「統計処理を行って、予め設定された判定条件情報に基づき状態変動の有無を判定する」ことの具体的一例として、上記状態変動判定手段41は、上記出力制御手段から出力される結果について単位ブレ距離を求め、求めた単位ブレ距離の平均値と偏差値を求め、偏差値に基づき状態変動の有無を判定する基礎情報を求めるという手法を挙げることができる。 The state change determination means 41 performs statistical processing on the result output from the state change detection auxiliary device 1A output control means, and determines the presence or absence of the state change based on the determination condition information set in advance. As a specific example of "performing statistical processing and determining the presence or absence of a state change based on preset determination condition information", the state change determination means 41 is a unit for a result output from the output control means. An example is a method in which the blur distance is obtained, the average value and the deviation value of the obtained unit blur distance are obtained, and the basic information for determining the presence or absence of a state change is obtained based on the deviation value.

更に具体的には、第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 auxiliary device 1A and the state change detection device 1B configured by the above means and the like execute the processing operation according to the flowcharts shown in FIGS. 9 and 10, the operation will be described with reference to the flowcharts. .. The flowchart of FIG. 9 shows a process related to the generation of the first prediction model 30A and the second prediction model 30B, and a process of obtaining the first importance and the second importance, and is a process performed by the CPU 10. Using the given teacher data T, the second teacher data is generated by the second machine learning data creating means 32 (S11). The second teacher data generation can be performed by a process of obtaining a transposed matrix. Next, the first prediction model 30A is generated using the teacher data T (S12). In step S12, 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.

更に、第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 first prediction model 30A (S14). In addition, the second importance data is obtained using the second teacher data and the second prediction model 30B (S15). In the processing of steps S14 and S15, a known method such as a method used for obtaining an appropriate decision tree when creating a prediction model for a random forest is used as a method for obtaining the importance. Can be done.

図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 CPU 10. The first analysis target data including the first objective variable data and the first explanatory variable data is acquired from a plurality of acquisition sources, and the second analysis target data is processed from the first analysis target data by the second machine learning data creation means 32. (S21). The first analysis target data is applied to the first prediction model 30A to obtain the first prediction value (S22), and the second analysis target data is applied to the second prediction model 30B to obtain the second prediction value (S23). In steps S22 and S23, the numerical value of the objective variable and / or the numerical value of the explanatory variable can be changed to a dummy numerical value by the dummy numerical value changing means 38, if necessary.

第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目的変数データと第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予測誤差値と、を求める予測誤差値算出手段と、
前記予測誤差値算出手段により算出された結果を出力手段へ出力する出力制御手段と
を具備することを特徴とする状態変動検出補助装置。
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.
前記第2機械学習用データ作成手段は、前記取得された第1解析対象データを、時間経過を示す表題情報を縦方向に複数並べ、取得元を示す取得元情報を横方向に複数並べて、各表題情報と各取得元情報の交差する位置に取得された値情報を配置した第1の行列として、この第1の行列を転置して転置行列を作成し第2解析対象データとすることを特徴とする請求項1に記載の状態変動検出補助装置。 The second machine learning data creating means arranges the acquired first analysis target data in a plurality of title information indicating the passage of time in the vertical direction, and arranges a plurality of acquisition source information indicating the acquisition source in the horizontal direction. As the first matrix in which the acquired value information is arranged at the intersection of the title information and each acquisition source information, the first matrix is transposed to create a transposed matrix and used as the second analysis target data. The state change detection assisting device according to claim 1. 前記第1予測モデルの第1重要度を求めて、この第1重要度で前記第1予測誤差値に伸長を図る処理を行うと共に、前記第2予測モデルの第2重要度を求めて、この第2重要度で前記第2予測誤差値を伸長する処理を行う誤差値伸長手段を具備し、
前記出力制御手段は、前記誤差値伸長手段による誤差値伸長の結果を出力することを特徴とする請求項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乃至3のいずれか1項に記載の状態変動検出補助装置。 Claims 1 to 1, wherein the numerical value of the objective variable and / or the numerical value of the explanatory variable is changed to a dummy numerical value, and a dummy numerical value changing means for bringing the numerical value of the objective variable and the numerical value of the explanatory variable close to each other is provided. 3. The state change detection assisting device according to any one of 3. 第1解析対象データにおいて、時間経過を示す表題情報が所定時間間隔で繰り返される情報であり、
前記出力手段は、前記表題情報を第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.
請求項1乃至4のいずれか1項に記載の状態変動検出補助装置と、
前記状態変動検出補助装置の前記出力制御手段から出力される結果について統計処理を行って、予め設定された判定条件情報に基づき状態変動の有無を判定する状態変動判定手段と、
を具備する特徴とする状態変動検出装置。
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.
前記状態変動判定手段は、前記出力制御手段から出力される結果について単位ブレ距離を求め、求めた単位ブレ距離の平均値と偏差値を求め、偏差値に基づき状態変動の有無を判定することを特徴とする請求項6に記載の状態変動検出装置。 The state change determining means obtains a unit blur distance from the result output from the output control means, obtains an average value and a deviation value of the obtained unit blur distance, and determines the presence or absence of a state change based on the deviation value. The state change detection device according to claim 6. コンピュータを、
第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.
前記コンピュータを前記第2機械学習用データ作成手段として、前記取得された第1解析対象データを、時間経過を示す表題情報を縦方向に複数並べ、取得元を示す取得元情報を横方向に複数並べて、各表題情報と各取得元情報の交差する位置に取得された値情報を配置した第1の行列として、この第1の行列を転置して転置行列を作成し第2解析対象データとするように機能させることを特徴とする請求項8に記載の状態変動検出補助用プログラム。 Using the computer as the second machine learning data creation means, a plurality of the acquired first analysis target data are arranged in the vertical direction with title information indicating the passage of time, and a plurality of acquisition source information indicating the acquisition source is arranged in the horizontal direction. As the first matrix in which the acquired value information is arranged at the intersection of each title information and each acquisition source information, the first matrix is transposed to create a transposed matrix and used as the second analysis target data. The state change detection assisting program according to claim 8, wherein the program functions in such a manner. 前記コンピュータを更に、
前記第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.
前記コンピュータを更に、前記目的変数の数値及び/または前記説明変数の数値をダミー数値へ変更し、前記目的変数の数値と前記説明変数の数値を近接させるダミー数値変更手段として機能させることを特徴とする請求項8乃至10のいずれか1項に記載の状態変動検出補助用プログラム。 The computer is further characterized by changing the numerical value of the objective variable and / or the numerical value of the explanatory variable to a dummy numerical value, and functioning as a dummy numerical value changing means for bringing the numerical value of the objective variable and the numerical value of the explanatory variable close to each other. The state change detection assisting program according to any one of claims 8 to 10. 第1解析対象データにおいて、時間経過を示す表題情報が所定時間間隔で繰り返される情報であり、前記コンピュータを前記出力手段として、前記表題情報を第1の軸とし、前記取得元を示す取得元情報を第2の軸とし、時間を第3の軸として、3次元の出力を行うように機能させることを特徴とする請求項8乃至11のいずれか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 acquisition source information indicating the acquisition source with the computer as the output means and the title information as the first axis. The state change detection assisting program according to any one of claims 8 to 11, wherein the program is made to function so as to perform three-dimensional output with time as the second axis and time as the third axis. コンピュータを、
第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.
前記コンピュータを前記状態変動判定手段として、前記出力制御手段から出力される結果について単位ブレ距離を求め、求めた単位ブレ距離の平均値と偏差値を求め、偏差値に基づき状態変動の有無を判定するように機能させることを特徴とする請求項13に記載の状態変動検出用プログラム。 Using the computer as the state change determining means, the unit blur distance is obtained for the result output from the output control means, the average value and the deviation value of the obtained unit shake distance are obtained, and the presence or absence of the state change is determined based on the deviation value. The state change detection program according to claim 13, wherein the program is made to function as such.
JP2019160751A 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 Active JP6961312B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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