特開2023-156850相関係数算出プログラム、情報処理装置および相関係数算出方法
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2023156850
(43)【公開日】2023-10-25
(54)【発明の名称】相関係数算出プログラム、情報処理装置および相関係数算出方法
(51)【国際特許分類】
G06F 17/15 20060101AFI20231018BHJP
【ＦＩ】
G06F17/15
【審査請求】未請求
【請求項の数】6
【出願形態】ＯＬ
(21)【出願番号】P 2022066461
(22)【出願日】2022-04-13
(71)【出願人】
【識別番号】000005223
【氏名又は名称】富士通株式会社
(74)【代理人】
【識別番号】110002147
【氏名又は名称】弁理士法人酒井国際特許事務所
(72)【発明者】
【氏名】松田 一仁
【テーマコード（参考）】
5B056
【Ｆターム（参考）】
5B056BB23
5B056FF03
(57)【要約】
【課題】相関が現れる条件を採用する場合に用いられる相関係数を効率的に算出する。
【解決手段】情報処理装置１は、観測値リスト２１から第１の属性の値を抜き出した第１の列データに条件Ｘをマスクして残った要素の平均値Ｘを算出し、第１の列データに条件Ｘを否定する条件Ｙをマスクして残った要素の平均値Ｙを算出し、第１の列データを第１のレジスタにロードし、平均値Ｘに条件Ｘをマスクした値と、平均値Ｙに条件Ｙをマスクした値とを第２のレジスタにロードし、演算器を用いて、第１のレジスタにロードされた値列と第２のレジスタにロードされた値列との第１の減算を行い、第２の列データについて、第１の減算を行う処理と同じ方法で第２の減算を行い、第１の減算によって得られた第１の値列と、第２の減算によって得られた第２の値列とを用いて、条件Ｘおよび条件Ｙにおける第１の列データと第２の列データとの間の相関係数を算出する。
【選択図】図１
【解決手段】情報処理装置１は、観測値リスト２１から第１の属性の値を抜き出した第１の列データに条件Ｘをマスクして残った要素の平均値Ｘを算出し、第１の列データに条件Ｘを否定する条件Ｙをマスクして残った要素の平均値Ｙを算出し、第１の列データを第１のレジスタにロードし、平均値Ｘに条件Ｘをマスクした値と、平均値Ｙに条件Ｙをマスクした値とを第２のレジスタにロードし、演算器を用いて、第１のレジスタにロードされた値列と第２のレジスタにロードされた値列との第１の減算を行い、第２の列データについて、第１の減算を行う処理と同じ方法で第２の減算を行い、第１の減算によって得られた第１の値列と、第２の減算によって得られた第２の値列とを用いて、条件Ｘおよび条件Ｙにおける第１の列データと第２の列データとの間の相関係数を算出する。
【選択図】図１
【特許請求の範囲】
【請求項1】
サンプル毎に各サンプルが持つ複数の属性の値を蓄積した表データから第１の属性の値を抜き出した第１の列データに第１の条件をマスクして残った要素の第１の平均値を算出し、
前記第１の列データに前記第１の条件を否定する第２の条件をマスクして残った要素の第２の平均値を算出し、
前記第１の列データを第１のレジスタにロードし、
前記第１の平均値に前記第１の条件をマスクした値と、前記第２の平均値に前記第２の条件をマスクした値とを第２のレジスタにロードし、
前記第１のレジスタと前記第２のレジスタを入力として演算器を用いて、前記第１のレジスタにロードされた値列と前記第２のレジスタにロードされた値列との第１の減算を行い、
前記第１の属性と異なる第２の属性の値を抜き出した第２の列データについて、前記第１の減算を行う処理と同じ方法で第２の減算を行い、
前記第１の減算によって得られた第１の値列と、前記第２の減算によって得られた第２の値列とを用いて、前記第１の条件および前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する
処理をコンピュータに実行させる相関係数算出プログラム。
前記第１の列データに前記第１の条件を否定する第２の条件をマスクして残った要素の第２の平均値を算出し、
前記第１の列データを第１のレジスタにロードし、
前記第１の平均値に前記第１の条件をマスクした値と、前記第２の平均値に前記第２の条件をマスクした値とを第２のレジスタにロードし、
前記第１のレジスタと前記第２のレジスタを入力として演算器を用いて、前記第１のレジスタにロードされた値列と前記第２のレジスタにロードされた値列との第１の減算を行い、
前記第１の属性と異なる第２の属性の値を抜き出した第２の列データについて、前記第１の減算を行う処理と同じ方法で第２の減算を行い、
前記第１の減算によって得られた第１の値列と、前記第２の減算によって得られた第２の値列とを用いて、前記第１の条件および前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する
処理をコンピュータに実行させる相関係数算出プログラム。
【請求項2】
前記相関係数を算出する処理は、
前記第１の値列を前記第１のレジスタおよび前記第２のレジスタにロードし、前記第１のレジスタと前記第２のレジスタを入力として前記演算器を用いて、前記第１の値列の２乗の第１の計算を算出し、
前記第２の値列を前記第１のレジスタおよび前記第２のレジスタにロードし、前記第１のレジスタと前記第２のレジスタを入力として前記演算器を用いて、前記第２の値列の２乗の第２の計算を算出し、
前記第１の値列を前記第１のレジスタにロードし、前記第２の値列を前記第２のレジスタにロードし、前記第１のレジスタと前記第２のレジスタを入力として前記演算器を用いて、前記第１の値列と前記第２の値列との乗算の第３の計算を算出し、
前記第１の計算の計算結果、前記第２の計算の計算結果、および前記第３の計算の計算結果を用いて、前記第１の条件および前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する
ことを特徴とする請求項１に記載の相関係数算出プログラム。
前記第１の値列を前記第１のレジスタおよび前記第２のレジスタにロードし、前記第１のレジスタと前記第２のレジスタを入力として前記演算器を用いて、前記第１の値列の２乗の第１の計算を算出し、
前記第２の値列を前記第１のレジスタおよび前記第２のレジスタにロードし、前記第１のレジスタと前記第２のレジスタを入力として前記演算器を用いて、前記第２の値列の２乗の第２の計算を算出し、
前記第１の値列を前記第１のレジスタにロードし、前記第２の値列を前記第２のレジスタにロードし、前記第１のレジスタと前記第２のレジスタを入力として前記演算器を用いて、前記第１の値列と前記第２の値列との乗算の第３の計算を算出し、
前記第１の計算の計算結果、前記第２の計算の計算結果、および前記第３の計算の計算結果を用いて、前記第１の条件および前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する
ことを特徴とする請求項１に記載の相関係数算出プログラム。
【請求項3】
前記第２のレジスタにロードする処理は、プレディケートレジスタを用いて、前記第１の平均値に前記第１の条件をマスクしてコピーした値と、前記第２の平均値に前記第２の条件をマスクしてコピーした値とを第２のレジスタにロードする
ことを特徴とする請求項１に記載の相関係数算出プログラム。
ことを特徴とする請求項１に記載の相関係数算出プログラム。
【請求項4】
前記相関係数を算出する処理は、
さらに、プリディケートレジスタを用いて、前記第１の計算の計算結果、前記第２の計算の計算結果、前記第３の計算の計算結果にそれぞれ前記第１の条件をマスクした値列を算出し、算出したそれぞれの値列を用いて、前記第１の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出し、
前記プリディケートレジスタを用いて、前記第１の計算の計算結果、前記第２の計算の計算結果、前記第３の計算の計算結果にそれぞれ前記第２の条件をマスクした値列を算出し、算出したそれぞれの値列を用いて、前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する
ことを特徴とする請求項２に記載の相関係数算出プログラム。
さらに、プリディケートレジスタを用いて、前記第１の計算の計算結果、前記第２の計算の計算結果、前記第３の計算の計算結果にそれぞれ前記第１の条件をマスクした値列を算出し、算出したそれぞれの値列を用いて、前記第１の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出し、
前記プリディケートレジスタを用いて、前記第１の計算の計算結果、前記第２の計算の計算結果、前記第３の計算の計算結果にそれぞれ前記第２の条件をマスクした値列を算出し、算出したそれぞれの値列を用いて、前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する
ことを特徴とする請求項２に記載の相関係数算出プログラム。
【請求項5】
サンプル毎に各サンプルが持つ複数の属性の値を蓄積した表データから第１の属性の値を抜き出した第１の列データに第１の条件をマスクして残った要素の第１の平均値を算出する第１の平均値計算部と、
前記第１の列データに前記第１の条件を否定する第２の条件をマスクして残った要素の第２の平均値を算出する第２の平均値計算部と、
前記第１の列データを第１のレジスタにロードし、前記第１の平均値に前記第１の条件をマスクした値と、前記第２の平均値に前記第２の条件をマスクした値とを第２のレジスタにロードする第１の偏差計算部と、
前記第１のレジスタと前記第２のレジスタを入力として演算器を用いて、前記第１のレジスタにロードされた値列と前記第２のレジスタにロードされた値列との第１の減算を行う第２の偏差計算部と、
前記第１の属性と異なる第２の属性の値を抜き出した第２の列データについて、前記第１の減算を行う処理と同じ方法で第２の減算を行う第３の偏差計算部と、
前記第１の減算によって得られた第１の値列と、前記第２の減算によって得られた第２の値列とを用いて、前記第１の条件および前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する相関係数計算部と、
を有することを特徴とする情報処理装置。
前記第１の列データに前記第１の条件を否定する第２の条件をマスクして残った要素の第２の平均値を算出する第２の平均値計算部と、
前記第１の列データを第１のレジスタにロードし、前記第１の平均値に前記第１の条件をマスクした値と、前記第２の平均値に前記第２の条件をマスクした値とを第２のレジスタにロードする第１の偏差計算部と、
前記第１のレジスタと前記第２のレジスタを入力として演算器を用いて、前記第１のレジスタにロードされた値列と前記第２のレジスタにロードされた値列との第１の減算を行う第２の偏差計算部と、
前記第１の属性と異なる第２の属性の値を抜き出した第２の列データについて、前記第１の減算を行う処理と同じ方法で第２の減算を行う第３の偏差計算部と、
前記第１の減算によって得られた第１の値列と、前記第２の減算によって得られた第２の値列とを用いて、前記第１の条件および前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する相関係数計算部と、
を有することを特徴とする情報処理装置。
【請求項6】
サンプル毎に各サンプルが持つ複数の属性の値を蓄積した表データから第１の属性の値を抜き出した第１の列データに第１の条件をマスクして残った要素の第１の平均値を算出し、
前記第１の列データに前記第１の条件を否定する第２の条件をマスクして残った要素の第２の平均値を算出し、
前記第１の列データを第１のレジスタにロードし、
前記第１の平均値に前記第１の条件をマスクした値と、前記第２の平均値に前記第２の条件をマスクした値とを第２のレジスタにロードし、
前記第１のレジスタと前記第２のレジスタを入力として演算器を用いて、前記第１のレジスタにロードされた値列と前記第２のレジスタにロードされた値列との第１の減算を行い、
前記第１の属性と異なる第２の属性の値を抜き出した第２の列データについて、前記第１の減算を行う処理と同じ方法で第２の減算を行い、
前記第１の減算によって得られた第１の値列と、前記第２の減算によって得られた第２の値列とを用いて、前記第１の条件および前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する
処理をコンピュータが実行する相関係数算出方法。
前記第１の列データに前記第１の条件を否定する第２の条件をマスクして残った要素の第２の平均値を算出し、
前記第１の列データを第１のレジスタにロードし、
前記第１の平均値に前記第１の条件をマスクした値と、前記第２の平均値に前記第２の条件をマスクした値とを第２のレジスタにロードし、
前記第１のレジスタと前記第２のレジスタを入力として演算器を用いて、前記第１のレジスタにロードされた値列と前記第２のレジスタにロードされた値列との第１の減算を行い、
前記第１の属性と異なる第２の属性の値を抜き出した第２の列データについて、前記第１の減算を行う処理と同じ方法で第２の減算を行い、
前記第１の減算によって得られた第１の値列と、前記第２の減算によって得られた第２の値列とを用いて、前記第１の条件および前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する
処理をコンピュータが実行する相関係数算出方法。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、相関係数算出プログラムなどに関する。
【背景技術】
【0002】
近年、相関関係を持つ条件の抽出により、因果探索すべき条件の数を効率的に絞り込む研究がなされている。図１７は、統計的因果探索のための条件抽出を示す参考図である。図１７に示すように、かかる技術では、ＡＩに学習させるための過去のデータから相関関係を持つデータの条件候補をすべて抽出する。そして、当該技術では、抽出したデータの条件候補から因果関係をもつデータの条件をすべて抽出する。ところが、当該技術では、すべてのデータの条件候補に対して因果関係を探索するが、計算量の観点から非現実的であるという問題がある。
【0003】
そこで、条件の探索対象を、因果関係から相関関係に緩和することで、因果探索すべき条件の数を効率的に絞り込む技術が開示されている（例えば、非特許文献）。図１８は、個々の特徴的な因果関係を発見する技術を示す参考図である。図１８に示すように、かかる技術では、まず、顕在パターン発見技術を用いて、過去のサンプル集合から、特定の条件下で目的変数と強い相関を持つ重要因子候補と、その時の条件との組を網羅的に求める。なお、過去のサンプル集合は、閾値に基づき２値化した後に使用される。
【0004】
その後、求めた条件それぞれに対して、因果探索技術を用いて、その条件下での重要因子候補が正しく重要因子であるかを判定する。例えば、「ｘ１∧ｘ３∧ｘ４→ｙ」（ｘ１＝ｘ３＝ｘ４＝１ならばｙ＝１）がある場合とする。かかる場合には、左辺から選んだ１変数を「重要因子の候補」、残りを「条件」とする。ここでは、ｘ４が「重要因子の候補」、残りの「ｘ１∧ｘ３」が「条件」を示すとする。かかる技術では、「条件」を満たす過去のサンプル集合の中で「重要因子の候補」と右辺のｙとに高い相関があれば、その「条件」を採用する。このように求めた条件と重要因子は、ＤＢ（DataBase）に保持される。そして、適用時には、因果関係を知りたいサンプルに対して、そのサンプルが満たす条件をＤＢから選択し、対応する重要因子を提示する。
【0005】
ここで、相関を取ろうとする２種類の信号を、ダイナミックレンジの中間値未満か否かに応じた１ビット信号に変換し、相関演算に関する演算量を低減させる技術が開示されている（例えば、特許文献１参照）。
【先行技術文献】
【特許文献】
【0006】
【特許文献1】特開２００８－１５８８５５号公報
【非特許文献】
【0007】
【非特許文献1】小柳 祐介，他４名，“個々の特徴的な因果関係を発見する技術の開発とマーケティングデータへの適用”，人工知能学会第１８回ビジネスインフォマティクス研究会，2021年3月，＜URL:http://sig-bi.jp/doc/18thSIG-BI2021/18thSIG-BI2021 paper13.pdf＞
【発明の概要】
【発明が解決しようとする課題】
【0008】
しかしながら、過去のサンプル集合を抽出するための条件を採用する場合に、「重要因子の候補」と右辺のｙとの間で相関係数を算出するが、相関係数を効率的に算出することが求められている。
【0009】
本発明は、１つの側面では、相関が現れる条件を採用する場合に用いられる相関係数を効率的に算出することを目的とする。
【課題を解決するための手段】
【0010】
１つの態様では、相関係数算出プログラムが、サンプル毎に各サンプルが持つ複数の属性の値を蓄積した表データから第１の属性の値を抜き出した第１の列データに第１の条件をマスクして残った要素の第１の平均値を算出し、前記第１の列データに前記第１の条件を否定する第２の条件をマスクして残った要素の第２の平均値を算出し、前記第１の列データを第１のレジスタにロードし、前記第１の平均値に前記第１の条件をマスクしてコピーした値と、前記第２の平均値に前記第２の条件をマスクしてコピーした値とを第２のレジスタにロードし、前記第１のレジスタと前記第２のレジスタを入力として演算器を用いて、前記第１のレジスタにロードされた値列と前記第２のレジスタにロードされた値列との第１の減算を行い、前記第１の属性と異なる第２の属性の値を抜き出した第２の列データについて、前記第１の減算を行う処理と同じ方法で第２の減算を行い、前記第１の減算によって得られた第１の値列と、前記第２の減算によって得られた第２の値列とをそれぞれ前記第１のレジスタまたは前記第２のレジスタにロードし、前記第１のレジスタと前記第２のレジスタを入力として前記演算器を用いて、前記第１の条件および前記第２の条件における前記第１の列データと前記第２の列データとの間の相関係数を算出する、処理をコンピュータに実行させる。
【発明の効果】
【0011】
１実施態様によれば、相関が現れる条件を採用する場合に用いられる相関係数を効率的に算出することができる。
【図面の簡単な説明】
【0012】
【発明を実施するための形態】
【0013】
以下に、本願の開示する相関係数算出プログラム、情報処理装置および相関係数算出方法の実施例を図面に基づいて詳細に説明する。なお、本発明は、実施例により限定されるものではない。
【実施例0014】
まず、個々のインスタンスに関して観測した複数事象の数値データ群に対して、「相関関係を持つ観測値ペアを有する」サンプル集合を抽出するための条件を列挙することを考える。図１４Ａ～図１４Ｃは、相関関係を持つサンプル群抽出のための条件列挙を説明する参考図である。
【0015】
図１４Ａ左図には、観測値リストが表わされている。観測値リストは、インスタンスｉｄに対して観測した複数の観測値（属性）の数値データ群を記憶したリストである。図１４Ａ右図には、条件リストが表わされている。条件リストは、「相関関係を持つ観測値ペアを有する」サンプル集合を抽出するための条件候補のリストである。条件は、例えば、観測値から生成されるが、これに限定されるものではない。条件は、網羅的な探索を目的とするため、１つの条件に対し、ｎｏｔ条件が含まれる。ここでは、１つの条件が「ａ＜５」である場合に、ｎｏｔ条件として「ｎｏｔ（ａ＜５）」が記憶される。
【0016】
図１４Ｂに示すように、条件式Ｘは、全条件からｋ個抜き出して作成される。ここでは、条件式Ｘは、「ａ＜５」と「ｎｏｔ（ｃ＝１）」との論理積を示す「（ａ＜５） ａｎｄ ｎｏｔ（ｃ＝１）」である。条件列挙処理は、条件リストから条件式Ｘを用いて論理積を計算し、サンプルのフィルタリングを行う。フィルタリングされた観測値リストが図１４Ｂ右下に表わされている。すなわち、条件式Ｘを用いて論理積を計算した結果が「１」を示すサンプル集合が表わされている。
【0017】
図１４Ｃには、条件式Ｘでフィルタリングされた観測値リストが表わされている。かかる観測値リストは、図１４Ｂ右下の図と同じである。条件列挙処理は、フィルタリングされた観測値リストから２つの観測値を抜き出して、抜き出した観測値ペアに対する相関係数を計算する。例えば、観測値ａと観測値ｂが抜き出され、観測値ａ－ｂ間の相関係数が計算される。観測値ａと観測値ｃが抜き出され、観測値ａ－ｃ間の相関係数が計算される。観測値ｂと観測値ｃが抜き出され、観測値ｂ－ｃ間の相関係数が計算される。そして、条件列挙処理は、相関係数が閾値を超える観測値ペアがｎ個以上ある等の判定を行い、合致する場合にフィルタリングに使用した条件を列挙する。ここでは、列挙された１つの条件は、「（ａ＜５） ａｎｄ ｎｏｔ（ｃ＝１）」となる。
【0018】
条件を網羅的に探索する場合、全条件数ｄｘ、抜き出される条件数ｋに対して、ｄｘＣｋパターンの条件式が判定されることになる。すなわち、それぞれの条件式が、「相関関係を持つ観測値ペアを有する」サンプル集合を抽出するための条件であるかどうかについて個別に判定されることになる。
【0019】
図１５は、観測値に対する相関係数演算を説明する参考図である。図１５に示すように、相関係数演算処理は、条件Ｘでフィルタリングされた観測値リストから２つの観測値を抜き出した行列Ｙ１、Ｙ２に関して、以下の式（１）～（４）を用いて相関係数ＲＹ１２を求める。なお、式（１）および式（２）に示されるＹ１ａｖｅは、行列Ｙ１の平均値である。Ｙ２ａｖｅは、行列Ｙ２の平均値である。
【数1】
【数2】
【数3】
【数4】
【0020】
すなわち、相関係数演算処理は、条件Ｘで行列Ｙ１および行列Ｙ２をマスクし、残った要素（Ｘ行列で１になっている行の要素）だけで相関係数を計算する。
【0021】
ところで、Ａｒｍ社のＳＶＥ（Scalable Vector Extension）では、ｐｒｅｄｉｃａｔｅレジスタに入力されたビット列によって、演算を有効にするか否かを制御できる。図１６は、ＳＩＭＤレジスタとｐｒｅｄｉｃａｔｅレジスタを用いた相関係数演算を説明する参考図である。なお、ここでは、相関係数演算の中の式（２）の演算について説明する。
【0022】
相関係数演算処理は、ＳＩＭＤレジスタａに行列Ｙ１をロードし、ＳＩＭＤレジスタｂに平均値Ｙ１ａｖｅをロードし、ｐｒｅｄｉｃａｔｅレジスタにＸ行列を与える。すると、相関係数演算処理は、ｐｒｅｄｉｃａｔｅレジスタに与えられたＸ行列によってＡＬＵ（演算器）がマスクされるので、図１４Ｂ右下に示されるフィルタリングされたＹを生成しないで、ＳＩＭＤ計算することができる。ここでは、Ｘ行列が「０１１０」である場合には、１ビット目と４ビット目に対するＡＬＵが不活性（ｉｎ－ａｃｔｉｖｅ）となり、式（２）が計算されない。そして、２ビット目と３ビット目に対するＡＬＵは活性（ａｃｔｉｖｅ）となり、式（２）が計算される。
【0023】
ところが、相関係数演算処理は、ｐｒｅｄｉｃａｔｅレジスタによってマスクをかけた部分のＡＬＵを有効活用できないので、ＣＰＵ利用率が低下するという問題がある。そこで、実施例では、ＡＬＵを有効活用できる相関係数演算処理について説明する。
【0024】
［情報処理装置の構成］
図１は、実施例に係る情報処理装置の構成を示す機能ブロック図である。
図１は、実施例に係る情報処理装置の構成を示す機能ブロック図である。
【0025】
情報処理装置１は、制御部１０と、記憶部２０とを有する。制御部１０は、平均値計算処理部１１と、偏差計算処理部１２と、相関係数計算部１３と、判定部１４とを有する。記憶部２０は、観測値リスト２１と、条件リスト２２とを有する。
【0026】
観測値リスト２１は、観測値リスト２１は、インスタンスｉｄに対して観測した複数の観測値の数値データ群を記憶したリストである。言い換えれば、観測値リスト２１は、各インスタンスｉｄが持つ複数の観測値（属性）の値を蓄積した表データである。ここでいうインスタンスｉｄとは、個々の人などを一意に識別する識別子である。観測値リスト２１の各列は、各観測値（属性）に対応する。ここで、観測値リスト２１の一例を、図２を参照して説明する。
【0027】
図２は、観測値リストの一例を示す図である。図２に示すように、観測値リスト２１は、インスタンスｉｄと、観測値ａと、観測値ｂと、観測値ｃと、・・・を対応付けたリストである。一例として、インスタンスｉｄが「１」である場合に、観測値ａとして「１．３」、観測値ｂとして「４．８」、観測値ｃとして「１」が記憶されている。インスタンスｉｄが「２」である場合に、観測値ａとして「２．１」、観測値ｂとして「３．７」、観測値ｃとして「０」が記憶されている。
【0028】
図１に戻って、条件リスト２２は、「相関関係を持つ観測値ペアを有する」インスタンスｉｄのサンプル集合を抽出するための条件候補のリストである。言い換えれば、条件リスト２２は、観測値リスト２１から、各インスタンスｉｄが持つ複数の観測値（属性）の値を条件候補に基づいて２値化した表データである。すなわち、条件リスト２２の列単位の配列は、条件に対するビット列となる。条件は、例えば、観測値から生成されるが、これに限定されるものではない。また、条件は、網羅的な探索を目的とするため、１つの条件に対し、ｎｏｔ条件が含まれる。ここで、条件リスト２２の一例を、図３を参照して説明する。
【0029】
図３は、条件リストの一例を示す図である。条件リスト２２は、インスタンスｉｄ、各種条件を対応付けたリストである。条件の一例として、「ａ＜５」、「ｂ＞４」、「ｃ＝１」が記憶されている。そして、これらの条件に対し、ｎｏｔ条件である「ｎｏｔ（ａ＜５）」、「ｎｏｔ（ｂ＞４）」、「ｎｏｔ（ｃ＝１）」が記憶されている。条件「ａ＜５」は、観測値ａが「５」より小さいという条件を示している。条件「ｎｏｔ（ａ＜５）」は、観測値ａが「５」より小さいという条件の否定条件、すなわち、観測値ａが「５」以上という条件を示している。
【0030】
そして、各種条件について、それぞれ条件を満たす場合には「１」が設定され、条件を満たさない場合には「０」が設定されている。値の一例として、インスタンスｉｄが「１」である場合には、条件「ａ＜５」について、観測値ａが「１．３」であり、観測値ａが「５」より小さいので、「１」が設定されている。条件「ｎｏｔ（ａ＜５）」について、観測値ａが「１．３」であり、観測値ａが「５」より小さいので、「０」が設定されている。そして、条件「ａ＜５」に対する列単位の配列は、「１１１・・・」のビット列となっている。条件「ｎｏｔ（ａ＜５）」に対する列単位の配列は、「０００・・・」のビット列となっている。
【0031】
図１に戻って、情報処理装置１は、観測値リスト２１の数値データ群に対して、「相関関係を持つ観測値ペアを有する」インスタンスｉｄのサンプル集合を抽出するための条件を列挙する。情報処理装置１は、条件を探索する際に、探索する複数の条件の中にある条件と対になる否定条件を含めることを利用して、ある条件とある条件の否定条件との相関係数の演算を一括して行う。
【0032】
ここで、相関係数の演算方法の原理を、図４Ａおよび図４Ｂを参照して説明する。図４Ａおよび図４Ｂは、実施例に係る演算方法の原理を示す図である。図４Ａに示すように、条件ｘ１に対し、必ず条件ｘ１の否定条件を含めることを利用する。例えば、ある条件式Ｘが（ｘ０∧ｘ１）である場合には、ある条件式Ｘに対し、ｘ１をｘ１の否定条件に置換した条件式Ｘ‘（ｘ０∧ｎｏｔ（ｘ１））も判定する条件として含める。そうすると、条件式Ｘに対応するビット列と条件式Ｘ´に対応するビット列とは、条件ｘ０を軸として排他的なビット列となる。
【0033】
つまり、条件式Ｘに対応するビット列は、条件ｘ０に対応するビット列と条件ｘ１に対応するビット列とのｂｉｔａｎｄ（論理積）で計算される。一方、条件式Ｘ´に対応するビット列は、条件ｘ０に対応するビット列と条件ｎｏｔ（ｘ１）に対応するビット列とのｂｉｔａｎｄ（論理積）で計算される。条件式Ｘに対応するビット列と条件式Ｘ´に対応するビット列との論理和は、条件ｘ０に対応するビット列となる。したがって、条件式Ｘと条件式Ｘ´とは、条件ｘ０を共通とした重複しない条件式であるので、一括して判定することができる。
【0034】
図４Ｂには、ＳＩＭＤレジスタとｐｒｅｄｉｃａｔｅレジスタを用いた式（２）に含まれる偏差計算を実行する模式図が表わされている。ここでは、条件式Ｘは、「０１１０」であり、条件式Ｘ´は、「１００１」であるとする。図４Ｂに示すように、ｐｒｅｄｉｃａｔｅレジスタに与えられた条件式Ｘに対応するビット列によってＡＬＵ（演算器）がマスクされるので、ＳＩＭＤレジスタにロードされた平均値Ｙ１Ｘａｖｅがマスクされないビットに対応するＡＬＵが活性化し、Ｃ［１］、Ｃ［２］に結果が算出される。また、ｐｒｅｄｉｃａｔｅレジスタに与えられた条件式Ｘ´のビット列によってＡＬＵ（演算器）がマスクされるので、ＳＩＭＤレジスタにロードされた平均値Ｙ１Ｘ´ａｖｅがマスクされないビットに対応するＡＬＵが活性化し、Ｃ´［１］、Ｃ´［２］に結果が算出される。このように、条件式Ｘに対応するビット列と条件式Ｘ´に対応するビット列とは、条件ｘ０を軸として排他的なビット列であるので、条件式Ｘと条件式Ｘ´とを一括して判定することができる。
【0035】
図１に戻って、平均値計算処理部１１は、観測値リスト２１から１つの観測値を抜き出した行列Ｙ１に関して、プレディケートレジスタおよびＳＩＭＤレジスタを利用して、条件式Ｘ（＝ｘ０∧ｘ１）をマスクして残った要素の平均値を計算する。また、平均値計算処理部１１は、同じ行列Ｙ１に関して、プレディケートレジスタおよびＳＩＭＤレジスタを利用して、条件式Ｘ´（＝ｘ０∧ｎｏｔ（ｘ１））をマスクして残った要素の平均値を計算する。条件式Ｘ´は、条件式Ｘを否定する条件式である。
【0036】
【0037】
図５に示すように、平均値計算処理部１１は、観測値リスト２１から第１の観測値の値を抜き出した行列Ｙ１をＳＩＭＤレジスタＡにロードする。平均値計算処理部１１は、行列Ｙ１の要素を条件式Ｘで２値化した第１のビット列をプレディケートレジスタにロードする。そして、平均値計算処理部１１は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、ＳＩＭＤレジスタＡの要素の第１の平均値Ｙ１Ｘａｖｅを算出する（符号ａ１）。すなわち、平均値計算処理部１１は、ＳＩＭＤレジスタＡの要素をプレディケートレジスタでマスクして残った要素の第１の平均値Ｙ１Ｘａｖｅを算出する。
【0038】
また、平均値計算処理部１１は、観測値リスト２１から第１の観測値の値を抜き出した行列Ｙ１の要素をＳＩＭＤレジスタＡにロードする。平均値計算処理部１１は、行列Ｙ１の要素を条件式Ｘ´で２値化した第２のビット列をプレディケートレジスタにロードする。そして、平均値計算処理部１１は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、ＳＩＭＤレジスタＡの要素の第２の平均値Ｙ１Ｘ´ａｖｅを算出する（符号ａ２）。すなわち、平均値計算処理部１１は、ＳＩＭＤレジスタＡの要素をプレディケートレジスタでマスクして残った要素の第２の平均値Ｙ１Ｘ´ａｖｅを算出する。
【0039】
また、平均値計算処理部１１は、条件式Ｘで２値化した第１のビット列をプレディケートレジスタにロードする。平均値計算処理部１１は、第１の平均値Ｙ１ＸａｖｅをＳＩＭＤレジスタＡにロードする。そして、平均値計算処理部１１は、第１の平均値Ｙ１Ｘａｖｅを、プレディケートレジスタでマスクしてＳＩＭＤレジスタＢにコピーする（符号ａ３）。また、平均値計算処理部１１は、条件式Ｘ´で２値化した第２のビット列をプレディケートレジスタにロードする。平均値計算処理部１１は、第２の平均値Ｙ１Ｘ´ａｖｅをＳＩＭＤレジスタＡにロードする。そして、平均値計算処理部１１は、第２の平均値Ｙ１Ｘ´ａｖｅを、プレディケートレジスタでマスクしてＳＩＭＤレジスタＢにコピーする（符号ａ４）。条件式Ｘで２値化した第１のビット列と条件式Ｘ´で２値化した第２のビット列とは、共通する条件ｘ０を軸として排他的なビット列であるので、ＳＩＭＤレジスタＢにコピーされた第１の平均値Ｙ１Ｘａｖｅおよび第２の平均値Ｙ１Ｘ´ａｖｅは同じビットにコピーされない。したがって、平均値計算処理部１１は、プレディケートレジスタを利用して、条件式Ｘと条件式Ｘ´のそれぞれでマスクされた第１の観測値の平均値Ｙ１Ｘａｖｅ、Ｙ１Ｘ´ａｖｅをＳＩＭＤレジスタＢにロードできる。
【0040】
また、平均値計算処理部１１は、行列Ｙ１の場合と同様に、観測値リスト２１から第２の観測値を抜き出した行列Ｙ２に関して、条件式Ｘ（＝ｘ０∧ｘ１）をマスクして残った要素の平均値Ｙ２Ｘａｖｅを計算すれば良い。また、平均値計算処理部１１は、行列Ｙ１の場合と同様に、同じ行列Ｙ２に関して、条件式Ｘ´（＝ｘ０∧ｎｏｔ（ｘ１））をマスクして残った要素の平均値Ｙ２Ｘ´ａｖｅを計算すれば良い。そして、平均値計算処理部１１は、プレディケートレジスタを利用して、条件式Ｘと条件式Ｘ´のそれぞれでマスクされた第２の観測値の平均値Ｙ２Ｘａｖｅ、Ｙ２Ｘ´ａｖｅをＳＩＭＤレジスタＢにロードすれば良い。
【0041】
図１に戻って、偏差計算処理部１２は、プレディケートレジスタおよびＳＩＭＤレジスタを利用して、観測値リスト２１から第１の観測値を抜き出した行列Ｙ１と、条件式ＸおよびＸ´それぞれでマスクされ得られた平均値の列との偏差を計算する。
【0042】
【0043】
図６に示すように、ＳＩＭＤレジスタＢには、平均値計算処理部１１によって条件式Ｘ、Ｘ´それぞれでマスクされた平均値（Ｙ１Ｘａｖｅ、Ｙ１Ｘ´ａｖｅ）の列が設定されている（符号ｂ２）。また、ＳＩＭＤレジスタＡには、第１の観測値の値を抜き出した行列Ｙ１の要素が設定されている（符号ｂ１）。偏差計算処理部１２は、条件式Ｘ、Ｘ´の共通の条件ｘ０をプレディケートレジスタにロードする（符号ｂ３）。そして、偏差計算処理部１２は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、ＳＩＭＤレジスタＡからＳＩＭＤレジスタＢを減算する演算を実行する。そして、偏差計算処理部１２は、減算結果をＳＩＭＤレジスタＣにロードする（符号ｂ４）。この結果、偏差計算処理部１２は、条件式Ｘおよび条件式Ｘ´に対する「Ｙ１－Ｙ１ａｖｅ」を計算できる。
【0044】
また、偏差計算処理部１２は、行列Ｙ１の場合と同様に、観測値リスト２１から第２の観測値を抜き出した行列Ｙ２についても、条件式Ｘおよび条件式Ｘ´に対する「Ｙ２－Ｙ２ａｖｅ」を計算すれば良い。このように、偏差計算処理部１２は、条件式Ｘ、Ｘ´の共通の条件ｘ０をプレディケートレジスタにロードして計算することで、条件式Ｘ、Ｘ´に対する偏差計算を一括して行うことができる。
【0045】
相関係数計算部１３は、プレディケートレジスタおよびＳＩＭＤレジスタを利用して、条件式Ｘに対する第１の観測値の行列Ｙ１と第２の観測値の行列Ｙ２との相関係数および条件式Ｘ´に対する行列Ｙ１と行列Ｙ２との相関係数を一括して計算する。
【0046】
【0047】
図７Ａに示すように、相関係数計算部１３は、第１の観測値の行列Ｙ１について、偏差計算処理部１２によって処理された条件式Ｘおよび条件式Ｘ´における「Ｙ１－Ｙ１ａｖｅ」の列値をＳＩＭＤレジスタＡおよびＳＩＭＤレジスタＢにロードする（符号ｃ１）。相関係数計算部１３は、条件式Ｘ、Ｘ´の共通の条件ｘ０をプレディケートレジスタにロードする（符号ｃ２）。そして、相関係数計算部１３は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、ＳＩＭＤレジスタＡとＳＩＭＤレジスタＢとを乗算する演算を実行する。そして、相関係数計算部１３は、乗算結果をＳＩＭＤレジスタＣにロードする（符号ｃ６）。この結果、相関係数計算部１３は、第１の観測値の行列Ｙ１について、条件式Ｘおよび条件式Ｘ´に対する式（２）を一括して計算して、Ｓｙ１を求める。
【0048】
また、相関係数計算部１３は、第２の観測値の行列Ｙ２について、偏差計算処理部１２によって処理された条件式Ｘおよび条件式Ｘ´における「Ｙ２－Ｙ２ａｖｅ」の列値をＳＩＭＤレジスタＡおよびＳＩＭＤレジスタＢにロードする（符号ｃ４）。そして、相関係数計算部１３は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、ＳＩＭＤレジスタＡとＳＩＭＤレジスタＢとを乗算する演算を実行する。そして、相関係数計算部１３は、乗算結果をＳＩＭＤレジスタＣにロードする（符号ｃ６）。この結果、相関係数計算部１３は、第２の観測値の行列Ｙ２について、条件式Ｘおよび条件式Ｘ´に対する式（３）を一括して計算して、Ｓｙ２を求める。
【0049】
また、相関係数計算部１３は、第１の観測値の行列Ｙ１について、条件式Ｘおよび条件式Ｘ´における「Ｙ１－Ｙ１ａｖｅ」の列値をＳＩＭＤレジスタＡにロードする（符号ｃ７）。相関係数計算部１３は、第２の観測値の行列Ｙ２について、条件式Ｘおよび条件式Ｘ´における「Ｙ２－Ｙ２ａｖｅ」の列値をＳＩＭＤレジスタＢにロードする（符号ｃ８）。そして、相関係数計算部１３は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、ＳＩＭＤレジスタＡとＳＩＭＤレジスタＢとを乗算する演算を実行する。そして、相関係数計算部１３は、乗算結果をＳＩＭＤレジスタＣにロードする（符号ｃ１０）。この結果、相関係数計算部１３は、第１の観測値の行列Ｙ１および第２の観測値の行列Ｙ２について、条件式Ｘおよび条件式Ｘ´に対する式（１）を一括して計算して、Ｓｙ１２を求める。
【0050】
図７Ｂに示すように、相関係数計算部１３は、プレディケートレジスタを利用して、条件式Ｘおよび条件式Ｘ´について一括して計算された結果を分ける。
【0051】
例えば、相関係数計算部１３は、プレディケートレジスタに、行列Ｙ１の要素を条件式Ｘで２値化した第１のビット列をロードする（符号ｄ１）。相関係数計算部１３は、式（２）によって計算された、条件式Ｘおよび条件式Ｘ´に対するＳｙ１をＳＩＭＤレジスタＡにロードする（符号ｄ２）。そして、相関係数計算部１３は、ＳＩＭＤレジスタＡをプレディケートレジスタでマスクして、マスクされない要素の和（ｓｕｍ）を計算して、「ｓｕｍＳｙ１Ｘ」を求める（符号ｄ３）。
【0052】
そして、相関係数計算部１３は、プレディケートレジスタに、行列Ｙ１の要素を条件式Ｘ´で２値化した第２のビット列をロードする（符号ｄ４）。そして、相関係数計算部１３は、ＳＩＭＤレジスタＡをプレディケートレジスタでマスクして、マスクされない要素の和（ｓｕｍ）を計算して、「ｓｕｍＳｙ１Ｘ´」を求める（符号ｄ５）。
【0053】
また、相関係数計算部１３は、同様に、式（２）によって計算された、条件式Ｘおよび条件式Ｘ´に対するＳｙ２から、「ｓｕｍＳｙ２Ｘ」を「ｓｕｍＳｙ２Ｘ´」を求める。相関係数計算部１３は、同様に、式（１）によって計算された、条件式および条件式Ｘ´に対するＳｙ１２から、「ｓｕｍＳｙ１２Ｘ」を「ｓｕｍＳｙ１２Ｘ´」を求める。
【0054】
そして、相関係数計算部１３は、図７Ｂ右図上段の式（４）に基づいて、行条件式Ｘに対する第１の観測値の行列Ｙ１と第２の観測値の行列Ｙ２との相関係数Ｒｙ１２Ｘを計算する。相関係数計算部１３は、図７Ｂ右図下段の式（４）に基づいて、条件式Ｘ´に対する第１の観測値の行列Ｙ１と第２の観測値の行列Ｙ２との相関係数Ｒｙ１２Ｘ´を計算する。
【0055】
同様にして、相関係数計算部１３は、条件式Ｘ、Ｘ´に対する別の観測値ペアの相関係数を計算する。これにより、相関係数計算部１３は、条件式Ｘに対し、否定条件Ｘ´を一括して計算可能であるため、相関係数演算のうち、式（２）で示されるＳｙ１，式（３）で示されるＳｙ２，式（１）で示されるＳｙ１２の計算に必要なサイクル数を半減させることができる。
【0056】
図１に戻って、判定部１４は、フィルタリングに利用した個々の条件に対し、相関係数が閾値を超える観測値ペアがｎ個以上あるか否かを判定する。ここでいうｎは、「１」でも良いし、予め定義されれば良い。そして、判定部１４は、相関係数が閾値を超える観測値ペアがｎ個以上あれば、フィルタリングに使用した条件を保存する。
【0057】
【0058】
図８Ａには、観測値リスト２１および条件リスト２２が表わされている。観測値リスト２１には、インスタンスｉｄに対応付けた観測値ａの値、観測値ｂの値、観測値ｃの値が記憶されている。すなわち、観測値ａの列単位の配列、観測値ｂの列単位の配列、観測値ｃの列単位の配列が記憶されている。そして、条件リスト２２には、インスタンスｉｄに対応付けて、各条件に対応する「０」または「１」が記憶されている。「０」は、条件に合致していないことを示し、「１」は、条件に合致していることを示している。また、条件リスト２２に含まれる条件には、対になるｎｏｔ条件が存在する。
【0059】
そして、条件リスト２２から条件式が生成される。図８Ｂでは、条件式は、全条件から３個抜き出して生成されるものとする。ここでは、条件式は、「ａ＜５」と「ｂ＞４」と「ｃ＝１」との論理積を示す「（ａ＜５） ａｎｄ （ｂ＞４） ａｎｄ （ｃ＝１）」である。この条件式をＸとする。また、３つ目の条件「ｃ＝１」を否定条件に変えた条件式「（ａ＜５） ａｎｄ （ｂ＞４） ａｎｄ ｎｏｔ（ｃ＝１）」をＸ´とする。そして、条件式Ｘを用いてｂｉｔａｎｄ（論理積）が算出される。条件式Ｘ´を用いてｂｉｔａｎｄ（論理積）が算出される。なお、ここでは、条件式Ｘ、Ｘ´の共通の条件であるｘ０「（ａ＜５） ａｎｄ （ｂ＞４）」のビット列は、「１１１１」となる。
【0060】
図８Ｃに示すように、平均値計算処理部１１は、相関係数を判定する観測値のペアを選択する。ここでは、観測値ａ、観測値ｂのペアが選択されるとする。平均値計算処理部１１は、条件式Ｘで２値化したビット列をプレディケートレジスタにロードする。そして、平均値計算処理部１１は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、観測値ａの平均値ＹａａｖｅＸを計算する。平均値計算処理部１１は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、観測値ｂの平均値ＹｂａｖｅＸを計算する。ここでは、（ＹａａｖｅＸ，ＹｂａｖｅＸ）は、（３．５（＝（２．１＋４．９）×０．５），３．６（＝（３．７＋３．５）×０．５））と計算される。
【0061】
同様に、平均値計算処理部１１は、条件式Ｘ´で２値化したビット列をプレディケートレジスタにロードする。そして、平均値計算処理部１１は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、観測値ａの平均値ＹａａｖｅＸ´を計算する。平均値計算処理部１１は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、観測値ｂの平均値ＹｂａｖｅＸ´を計算する。ここでは、（ＹａａｖｅＸ´，ＹｂａｖｅＸ´）は、（１．８５（＝（１．３＋２．４）×０．５），４．０（＝（４．８＋３．２）×０．５））と計算される。
【0062】
図８Ｄに示すように、平均値計算処理部１１は、条件式Ｘで２値化したビット列「０１０１」をプレディケートレジスタにロードする。そして、平均値計算処理部１１は、平均値ＹａａｖｅＸを、プレディケートレジスタでマスクしてＳＩＭＤレジスタｚ０にコピーする。ここでは、ＳＩＭＤレジスタｚ０の２ビット目と４ビット目に平均値ＹａａｖｅＸ「３．５」がコピーされている。
【0063】
また、平均値計算処理部１１は、条件式Ｘ´で２値化したビット列「１０１０」をプレディケートレジスタにロードする。そして、平均値計算処理部１１は、平均値ＹａａｖｅＸ´を、プレディケートレジスタでマスクしてＳＩＭＤレジスタｚ０にコピーする。ここでは、ＳＩＭＤレジスタｚ０の１ビット目と３ビット目に平均値ＹａａｖｅＸ´「１．８５」がコピーされている。
【0064】
図８Ｅに示すように、偏差計算処理部１２は、条件式Ｘ、Ｘ´の共通の条件ｘ０で２値化したビット列「１１１１」をプレディケートレジスタにロードする。そして、偏差計算処理部１２は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、観測値ａの値を抜き出した行列ＹａからＳＩＭＤレジスタｚ０を減算する演算を実行する。そして、偏差計算処理部１２は、減算結果をＳＩＭＤレジスタＳａｔｍｐにロードする。すなわち、偏差計算処理部１２は、条件式Ｘおよび条件式Ｘ´に対する「Ｙａ－Ｙａａｖｅ」を一括して計算できる。
【0065】
観測値ｂについても、観測値ａの場合と同様に、偏差計算処理部１２は、演算器（ＡＬＵ）をプレディケートレジスタでマスクして、観測値ｂの値を抜き出した行列ＹｂからＳＩＭＤレジスタｚ０を減算する演算を実行する。そして、偏差計算処理部１２は、減算結果をＳＩＭＤレジスタＳｂｔｍｐにロードする。すなわち、偏差計算処理部１２は、条件式Ｘおよび条件式Ｘ´に対する「Ｙｂ－Ｙｂａｖｅ」を一括して計算できる。
【0066】
図８Ｆに示すように、相関係数計算部１３は、プレディケートレジスタを利用して、ＳＩＭＤレジスタＳａｔｍｐの値列の２乗を計算して、ＳＩＭＤレジスタＳａにロードする。かかる場合のプレディケートレジスタには、条件式Ｘ、Ｘ´の共通の条件ｘ０で２値化したビット列「１１１１」がロードされる。この結果、相関係数計算部１３は、観測値ａの行列Ｙａについて、条件式Ｘおよび条件式Ｘ´に対する式（２）を一括して計算して、Ｓａを求める。
【0067】
観測値ｂについても、観測値ａの場合と同様に、相関係数計算部１３は、プレディケートレジスタを利用して、ＳＩＭＤレジスタＳｂｔｍｐの値列の２乗を計算して、ＳＩＭＤレジスタＳｂにロードする。かかる場合のプレディケートレジスタには、条件式Ｘ、Ｘ´の共通の条件ｘ０で２値化したビット列「１１１１」がロードされる。この結果、相関係数計算部１３は、観測値ｂの行列Ｙｂについて、条件式Ｘおよび条件式Ｘ´に対する式（２）を一括して計算して、Ｓｂを求める。
【0068】
また、相関係数計算部１３は、プレディケートレジスタを利用して、ＳＩＭＤレジスタＳａｔｍｐの値列とＳＩＭＤレジスタＳｂｔｍｐの値列とを乗算して、ＳＩＭＤレジスタＳａｂにロードする。かかる場合のプレディケートレジスタには、条件式Ｘ、Ｘ´の共通の条件ｘ０で２値化したビット列「１１１１」がロードされる。この結果、相関係数計算部１３は、観測値ａの行列Ｙａおよび観測値ｂの行列Ｙｂについて、条件式Ｘおよび条件式Ｘ´に対する式（１）を一括して計算して、Ｓａｂを求める。
【0069】
そして、相関係数計算部１３は、プレディケートレジスタに条件式Ｘで２値化したビット列「０１０１」をロードする。そして、相関係数計算部１３は、条件式Ｘおよび条件式Ｘ´に対するＳａをプレディケートレジスタでマスクして、マスクされない要素の和（ｓｕｍ）を計算して、「ｓｕｍＳａＸ」を求める。相関係数計算部１３は、条件式Ｘおよび条件式Ｘ´に対するＳｂをプレディケートレジスタでマスクして、マスクされない要素の和（ｓｕｍ）を計算して、「ｓｕｍＳｂＸ」を求める。相関係数計算部１３は、条件式Ｘおよび条件式Ｘ´に対するＳａｂをプレディケートレジスタでマスクして、マスクされない要素の和（ｓｕｍ）を計算して、「ｓｕｍＳａｂＸ」を求める。
【0070】
そして、相関係数計算部１３は、プレディケートレジスタに条件式Ｘ´で２値化したビット列「１０１０」をロードする。そして、相関係数計算部１３は、条件式Ｘおよび条件式Ｘ´に対するＳａをプレディケートレジスタでマスクして、マスクされない要素の和（ｓｕｍ）を計算して、「ｓｕｍＳａＸ´」を求める。相関係数計算部１３は、条件式Ｘおよび条件式Ｘ´に対するＳｂをプレディケートレジスタでマスクして、マスクされない要素の和（ｓｕｍ）を計算して、「ｓｕｍＳｂＸ´」を求める。相関係数計算部１３は、条件式Ｘおよび条件式Ｘ´に対するＳａｂをプレディケートレジスタでマスクして、マスクされない要素の和（ｓｕｍ）を計算して、「ｓｕｍＳａｂＸ´」を求める。
【0071】
そして、相関係数計算部１３は、（ＳａＸ，ＳｂＸ，ＳａｂＸ）を用いて、条件式Ｘに対する、観測値ａと観測値ｂとの相関係数ＲｙａｂＸを計算する。相関係数計算部１３は、（ＳａＸ´，ＳｂＸ´，ＳａｂＸ´）を用いて、条件式Ｘ´に対する観測値ａと観測値ｂとの相関係数ＲｙａｂＸ´を計算する。
【0072】
［相関係数計算処理のフローチャート］
図９は、実施例に係る相関係数計算処理の全体のフローチャートの一例を示す図である。なお、観測値リスト２１および条件リスト２２は、記憶部２０に格納されている。そして、情報処理装置１は、観測値リスト２１から選択された観測値行列Ｙ１、Ｙ２および条件リスト２２から選択された条件ｘ０、ｘ１を受け付けたとする。観測値行列Ｙ１、Ｙ２は、観測値リスト２１から２つの観測値の行列を抜き出したものである。
図９は、実施例に係る相関係数計算処理の全体のフローチャートの一例を示す図である。なお、観測値リスト２１および条件リスト２２は、記憶部２０に格納されている。そして、情報処理装置１は、観測値リスト２１から選択された観測値行列Ｙ１、Ｙ２および条件リスト２２から選択された条件ｘ０、ｘ１を受け付けたとする。観測値行列Ｙ１、Ｙ２は、観測値リスト２１から２つの観測値の行列を抜き出したものである。
【0073】
情報処理装置１は、条件ｘ０、ｘ１を用いて、観測値行列Ｙ１、Ｙ２における平均値計算処理を実行する（ステップＳ１１）。なお、平均値計算処理のフローチャートについては、後述する。
【0074】
情報処理装置１は、条件ｘ０、ｘ１を用いて、観測値行列Ｙ１、Ｙ２における偏差計算処理を実行する（ステップＳ１２）。なお、偏差計算処理のフローチャートについては、後述する。
【0075】
情報処理装置１は、条件ｘ０、ｘ１を用いて、観測値行列Ｙ１、Ｙ２における相関係数処理および判定処理を実行する（ステップＳ１３）。なお、相関係数処理／判定処理のフローチャートについては、後述する。そして、情報処理装置１は、相関係数計算処理のプロセスを終了する。
【0076】
［平均値計算処理のフローチャート］
図１０は、実施例に係る平均値計算処理のフローチャートの一例を示す図である。なお、フローチャートで示すｐレジスタは、プレディケートレジスタである。
図１０は、実施例に係る平均値計算処理のフローチャートの一例を示す図である。なお、フローチャートで示すｐレジスタは、プレディケートレジスタである。
【0077】
まず、平均値計算処理部１１は、観測値行列Ｙ１、Ｙ２の全要素を順番に処理するまで、ステップＳ２１～Ｓ２８をループする。平均値計算処理部１１は、観測値行列Ｙ１をＳＩＭＤレジスタｚ０に最大数ロードする（ステップＳ２１）。平均値計算処理部１１は、観測値行列Ｙ２をＳＩＭＤレジスタｚ１に最大数ロードする（ステップＳ２２）。
【0078】
そして、平均値計算処理部１１は、条件式Ｘ「ｘ０∧ｘ１」に対応するビット列をｐレジスタに最大数ロードする(ステップＳ２３）。そして、平均値計算処理部１１は、ＳＩＭＤレジスタｚ０の処理対象の位置の要素をｐレジスタでマスクした上でＹ１Ｘａｖｅに加算する（ステップＳ２４）。平均値計算処理部１１は、ＳＩＭＤレジスタｚ１の処理対象の同じ位置の要素をｐレジスタでマスクした上でＹ２Ｘａｖｅに加算する（ステップＳ２５）。
【0079】
そして、平均値計算処理部１１は、条件式Ｘ´「ｘ０∧ｎｏｔ（ｘ１）」に対応するビット列をｐレジスタに最大数ロードする(ステップＳ２６）。そして、平均値計算処理部１１は、ＳＩＭＤレジスタｚ０の処理対象の位置の要素をｐレジスタでマスクした上でＹ１Ｘ´ａｖｅに加算する（ステップＳ２７）。平均値計算処理部１１は、ＳＩＭＤレジスタｚ１の処理対象の同じ位置の要素をｐレジスタでマスクした上でＹ２Ｘ´ａｖｅに加算する（ステップＳ２８）。
【0080】
続いて、平均値計算処理部１１は、観測値行列Ｙ１、Ｙ２の全要素を処理すると、Ｙ１Ｘａｖｅ、Ｙ１Ｘ´ａｖｅを、観測値行列Ｙ１の要素数で除算する（ステップＳ２９）。すなわち、平均値計算処理部１１は、観測値行列Ｙ１の要素を条件式Ｘに対応するビット列でマスクして残った要素の平均値Ｙ１Ｘａｖｅを算出する。平均値計算処理部１１は、観測値行列Ｙ１の要素を条件式Ｘ´に対応するビット列でマスクして残った要素の平均値Ｙ１Ｘ´ａｖｅを算出する。
【0081】
加えて、平均値計算処理部１１は、Ｙ２Ｘａｖｅ、Ｙ２Ｘ´ａｖｅを、観測値行列Ｙ２の要素数で除算する（ステップＳ３０）。平均値計算処理部１１は、観測値行列Ｙ２の要素を条件式Ｘに対応するビット列でマスクして残った要素の平均値Ｙ２Ｘａｖｅを算出する。平均値計算処理部１１は、観測値行列Ｙ２の要素を条件式Ｘ´に対応するビット列でマスクして残った要素の平均値Ｙ２Ｘ´ａｖｅを算出する。
【0082】
そして、平均値計算処理部１１は、平均値計算処理を終了する。
【0083】
［偏差計算処理のフローチャート］
図１１は、実施例に係る偏差計算処理のフローチャートの一例を示す図である。
図１１は、実施例に係る偏差計算処理のフローチャートの一例を示す図である。
【0084】
偏差計算処理部１２は、観測値行列Ｙ１の全要素を順番に処理するまで、ステップＳ４１～Ｓ４７をループする。偏差計算処理部１２は、条件式Ｘ「ｘ０∧ｘ１」に対応するビット列をｐレジスタに最大数ロードする(ステップＳ４１）。偏差計算処理部１２は、平均値Ｙ１Ｘａｖｅをｐレジスタでマスクした上でＳＩＭＤレジスタｚ０に最大数コピーする（ステップＳ４２）。
【0085】
そして、偏差計算処理部１２は、条件式Ｘ´「ｘ０∧ｎｏｔ（ｘ１）」に対応するビット列をｐレジスタに最大数ロードする(ステップＳ４３）。偏差計算処理部１２は、平均値Ｙ１Ｘ´ａｖｅをｐレジスタでマスクした上でＳＩＭＤレジスタｚ０に最大数コピーする（ステップＳ４４）。
【0086】
そして、偏差計算処理部１２は、条件式ＸとＸ´の共通の条件ｘ０に対応するビット列をｐレジスタに最大数ロードする（ステップＳ４５）。偏差計算処理部１２は、観測値行列Ｙ１をＳＩＭＤレジスタｚ１に最大数ロードする（ステップＳ４６）。そして、偏差計算処理部１２は、ｐレジスタでマスクした上でＳＩＭＤレジスタｚ１からＳＩＭＤレジスタｚ０を減算し、一次配列Ｓ１ｔｍｐに格納する（ステップＳ４７）。すなわち、偏差計算処理部１２は、条件式Ｘおよび条件式Ｘ´に対する「観測値行列Ｙ１－Ｙ１ａｖｅ」を計算して、一次配列Ｓ１ｔｍｐに格納する。
【0087】
続いて、偏差計算処理部１２は、観測値行列Ｙ１の全要素を処理すると、観測値行列Ｙ２の全要素を順番に処理するまで、ステップＳ４８～Ｓ５４をループする。偏差計算処理部１２は、条件式Ｘ「ｘ０∧ｘ１」に対応するビット列をｐレジスタに最大数ロードする(ステップＳ４８）。偏差計算処理部１２は、平均値Ｙ２Ｘａｖｅをｐレジスタでマスクした上でＳＩＭＤレジスタｚ０に最大数コピーする（ステップＳ４９）。
【0088】
そして、偏差計算処理部１２は、条件式Ｘ´「ｘ０∧ｎｏｔ（ｘ１）」に対応するビット列をｐレジスタに最大数ロードする(ステップＳ５０）。偏差計算処理部１２は、平均値Ｙ２Ｘ´ａｖｅをｐレジスタでマスクした上でＳＩＭＤレジスタｚ０に最大数コピーする（ステップＳ５１）。
【0089】
そして、偏差計算処理部１２は、条件式ＸとＸ´の共通の条件ｘ０に対応するビット列をｐレジスタに最大数ロードする（ステップＳ５２）。偏差計算処理部１２は、観測値行列Ｙ２をＳＩＭＤレジスタｚ１に最大数ロードする（ステップＳ５３）。そして、偏差計算処理部１２は、ｐレジスタでマスクした上でＳＩＭＤレジスタｚ１からＳＩＭＤレジスタｚ０を減算し、一次配列Ｓ２ｔｍｐに格納する（ステップＳ５４）。すなわち、偏差計算処理部１２は、条件式Ｘおよび条件式Ｘ´に対する「観測値行列Ｙ２－Ｙ２ａｖｅ」を計算して、一次配列Ｓ２ｔｍｐに格納する。
【0090】
［相関係数計算／判定処理のフローチャート］
図１２は、実施例に係る相関係数計算／判定処理のフローチャートの一例を示す図である。
図１２は、実施例に係る相関係数計算／判定処理のフローチャートの一例を示す図である。
【0091】
図１２に示すように、相関係数計算部１３は、「観測値行列Ｙ１－Ｙ１ａｖｅ」を格納した一次配列Ｓ１ｔｍｐ，「観測値行列Ｙ２－Ｙ２ａｖｅ」を格納した一次配列Ｓ２ｔｍｐの全要素を順番に処理するまで、ステップＳ６１～Ｓ７０をループする。相関係数計算部１３は、「観測値行列Ｙ１－Ｙ１ａｖｅ」を格納した一次配列Ｓ１ｔｍｐをＳＩＭＤレジスタｚ０に最大数ロードする（ステップＳ６１）。相関係数計算部１３は、「観測値行列Ｙ２－Ｙ２ａｖｅ」を格納した一次配列Ｓ２ｔｍｐをＳＩＭＤレジスタｚ１に最大数ロードする（ステップＳ６２）。
【0092】
相関係数計算部１３は、条件式Ｘ「ｘ０∧ｘ１」に対応するビット列をｐレジスタに最大数ロードする(ステップＳ６３）。相関係数計算部１３は、ＳＩＭＤレジスタｚ０をｐレジスタでマスクした上で自乗し、Ｓ１Ｘに加算する（ステップＳ６４）。すなわち、相関係数計算部１３は、条件式Ｘに対応する「ｓｕｍ（観測値行列Ｙ１－Ｙ１ａｖｅ）２」を計算する。
【0093】
相関係数計算部１３は、ＳＩＭＤレジスタｚ１をｐレジスタでマスクした上で自乗し、Ｓ２Ｘに加算する（ステップＳ６５）。すなわち、相関係数計算部１３は、条件式Ｘに対応する「ｓｕｍ（観測値行列Ｙ２－Ｙ２ａｖｅ）２」を計算する。相関係数計算部１３は、ｐレジスタでマスクした上で、ＳＩＭＤレジスタｚ０とＳＩＭＤレジスタｚ１とを乗算し、Ｓ１２Ｘに加算する（ステップＳ６６）。すなわち、相関係数計算部１３は、条件式Ｘに対応する「ｓｕｍ｛（観測値行列Ｙ１－Ｙ１ａｖｅ）×（観測値行列Ｙ２－Ｙ２ａｖｅ）｝」を計算する。
【0094】
相関係数計算部１３は、条件式Ｘ「ｘ０∧ｎｏｔ（ｘ１）」に対応するビット列をｐレジスタに最大数ロードする(ステップＳ６７）。相関係数計算部１３は、ＳＩＭＤレジスタｚ０をｐレジスタでマスクした上で自乗し、Ｓ１Ｘ´に加算する（ステップＳ６８）。すなわち、相関係数計算部１３は、条件式Ｘ´に対応する「ｓｕｍ（観測値行列Ｙ１－Ｙ１ａｖｅ）２」を計算する。
【0095】
相関係数計算部１３は、ＳＩＭＤレジスタｚ１をｐレジスタでマスクした上で自乗し、Ｓ２Ｘ´に加算する（ステップＳ６９）。すなわち、相関係数計算部１３は、条件式Ｘ´に対応する「ｓｕｍ（観測値行列Ｙ２－Ｙ２ａｖｅ）２」を計算する。相関係数計算部１３は、ｐレジスタでマスクした上で、ＳＩＭＤレジスタｚ０とＳＩＭＤレジスタｚ１とを乗算し、Ｓ１２Ｘ´に加算する（ステップＳ７０）。すなわち、相関係数計算部１３は、条件式Ｘ´に対応する「ｓｕｍ｛（観測値行列Ｙ１－Ｙ１ａｖｅ）×（観測値行列Ｙ２－Ｙ２ａｖｅ）｝」を計算する。
【0096】
続いて、相関係数計算部１３は、一次配列Ｓ１ｔｍｐ，一次配列Ｓ２ｔｍｐの全要素を処理すると、条件式Ｘに対応するＳ１Ｘ、Ｓ２Ｘ、Ｓ１２Ｘから相関係数Ｒ１２Ｘを計算し、判定部１４は、閾値の判定処理を行う（ステップＳ７１）。例えば、判定部１４は、式（１）に基づいて相関係数Ｒ１２Ｘを計算する。そして、判定部１４は、相関係数Ｒ１２Ｘが閾値を超えるか否かを判定し、判定結果を保持する。また、相関係数計算部１３は、条件式Ｘ´に対応するＳ１Ｘ´、Ｓ２Ｘ´、Ｓ１２Ｘ´から相関係数Ｒ１２Ｘ´を計算し、判定部１４は、閾値の判定処理を行う（ステップＳ７２）。例えば、判定部１４は、式（１）に基づいて相関係数Ｒ１２Ｘ´を計算する。そして、判定部１４は、相関係数Ｒ１２Ｘ´が閾値を超えるか否かを判定し、判定結果を保持する。
【0097】
この後、相関係数計算部１３は、観測値行列Ｙ１、Ｙ２のペアを変えて、条件式Ｘ，Ｘ´に対応する相関係数を計算する。そして、判定部１４は、相関係数が閾値を超える観測値ペアが予め定められたｎ個以上あるか否かを判定し、判定条件に合致する場合には、抽出に使用した条件式を保存する。
【0098】
このようにして、相関係数計算部１３は、条件Ｘに対応するＳ１２Ｘ，Ｓ１Ｘ、Ｓ２Ｘの計算を、条件Ｘ´に対応するＳ１２Ｘ´，Ｓ１Ｘ´、Ｓ２Ｘ´の計算と一括して行うことができるので、条件Ｘおよび条件Ｘ´における相関係数を効率的に計算できる。この結果、相関係数計算部１３は、条件式の判定に至る工程を従来に比較して迅速に行うことができる。
【0099】
［実施例の効果］
上記実施例によれば、情報処理装置１は、サンプル毎に各サンプルが持つ複数の属性の値を蓄積した表データから第１の属性の値を抜き出した第１の列データに第１の条件をマスクして残った要素の第１の平均値を算出する。情報処理装置１は、第１の列データに第１の条件を否定する第２の条件をマスクして残った要素の第２の平均値を算出する。情報処理装置１は、第１の列データを第１のレジスタにロードする。情報処理装置１は、第１の平均値に第１の条件をマスクした値と、第２の平均値に第２の条件をマスクした値とを第２のレジスタにロードする。情報処理装置１は、第１のレジスタと第２のレジスタを入力として演算器を用いて、第１のレジスタにロードされた値列と第２のレジスタにロードされた値列との第１の減算を行う。情報処理装置１は、第１の属性と異なる第２の属性の値を抜き出した第２の列データについて、第１の減算を行う処理と同じ方法で第２の減算を行う。そして、情報処理装置１は、第１の減算によって得られた第１の値列と、第２の減算によって得られた第２の値列とを用いて、第１の条件および第２の条件における第１の列データと第２の列データとの間の相関係数を算出する。かかる構成によれば、情報処理装置１は、相関関係を持つ属性のペアを有するサンプル群を抽出するための条件を見つける場合に、第１の条件および第１の条件を否定する第２の条件における属性のペアの相関係数を、効率良く演算器を活用して計算できる。すなわち、情報処理装置１は、第１の条件と第２の条件における属性のペアの相関係数の減算の計算について、個々の条件の演算では不活性であった演算器を活用できることとなり、効率良く演算器を活用して計算できる。この結果、情報処理装置１は、ＣＰＵ利用率を向上させることができる。
上記実施例によれば、情報処理装置１は、サンプル毎に各サンプルが持つ複数の属性の値を蓄積した表データから第１の属性の値を抜き出した第１の列データに第１の条件をマスクして残った要素の第１の平均値を算出する。情報処理装置１は、第１の列データに第１の条件を否定する第２の条件をマスクして残った要素の第２の平均値を算出する。情報処理装置１は、第１の列データを第１のレジスタにロードする。情報処理装置１は、第１の平均値に第１の条件をマスクした値と、第２の平均値に第２の条件をマスクした値とを第２のレジスタにロードする。情報処理装置１は、第１のレジスタと第２のレジスタを入力として演算器を用いて、第１のレジスタにロードされた値列と第２のレジスタにロードされた値列との第１の減算を行う。情報処理装置１は、第１の属性と異なる第２の属性の値を抜き出した第２の列データについて、第１の減算を行う処理と同じ方法で第２の減算を行う。そして、情報処理装置１は、第１の減算によって得られた第１の値列と、第２の減算によって得られた第２の値列とを用いて、第１の条件および第２の条件における第１の列データと第２の列データとの間の相関係数を算出する。かかる構成によれば、情報処理装置１は、相関関係を持つ属性のペアを有するサンプル群を抽出するための条件を見つける場合に、第１の条件および第１の条件を否定する第２の条件における属性のペアの相関係数を、効率良く演算器を活用して計算できる。すなわち、情報処理装置１は、第１の条件と第２の条件における属性のペアの相関係数の減算の計算について、個々の条件の演算では不活性であった演算器を活用できることとなり、効率良く演算器を活用して計算できる。この結果、情報処理装置１は、ＣＰＵ利用率を向上させることができる。
【0100】
また、上記実施例によれば、情報処理装置１は、第１の値列を第１のレジスタおよび第２のレジスタにロードし、第１のレジスタと第２のレジスタを入力として演算器を用いて、第１の値列の２乗の第１の計算を算出する。情報処理装置１は、第２の値列を第１のレジスタおよび第２のレジスタにロードし、第１のレジスタと第２のレジスタを入力として演算器を用いて、第２の値列の２乗の第２の計算を算出する。情報処理装置１は、第１の値列を第１のレジスタにロードし、第２の値列を第２のレジスタにロードし、第１のレジスタと第２のレジスタを入力として演算器を用いて、第１の値列と第２の値列との乗算の第３の計算を算出する。そして、情報処理装置１は、第１の計算結果、第２の計算結果、および第３の計算結果を用いて、第１の条件および第２の条件における第１の列データと第２の列データとの間の相関係数を算出する。かかる構成によれば、情報処理装置１は、第１の条件と第２の条件における属性のペアの相関係数の減算結果を乗算する計算について、個々の条件の演算では不活性であった演算器を活用できることとなり、効率良く演算器を活用して計算できる。
【0101】
また、上記実施例によれば、情報処理装置１は、プレディケートレジスタを用いて、第１の平均値に第１の条件をマスクしてコピーした値と、第２の平均値に第２の条件をマスクしてコピーした値とを第２のレジスタにロードする。かかる構成によれば、情報処理装置１は、プレディケートレジスタを用いることで、第１の条件および第２の条件をそれぞれマスクしてそれぞれの平均値を同一のレジスタにロードすることができる。
【0102】
また、上記実施例によれば、情報処理装置１は、さらに、プリディケートレジスタを用いて、第１の計算結果、第２の計算結果、第３の計算結果にそれぞれ第１の条件をマスクした値列を算出し、算出したそれぞれの値列を用いて、第１の条件における第１の列データと第２の列データとの間の相関係数を算出する。そして、情報処理装置１は、プリディケートレジスタを用いて、第１の計算結果、第２の計算結果、第３の計算結果にそれぞれ第２の条件をマスクした値列を算出し、算出したそれぞれの値列を用いて、第２の条件における第１の列データと第２の列データとの間の相関係数を算出する。かかる構成によれば、情報処理装置１は、プレディケートレジスタを用いることで、第１の条件と第２の条件をそれぞれマスクして計算結果を分離することができる。この結果、情報処理装置１は、第１の条件と第２の条件における属性のペアの相関係数を効率良く計算することができる。
【0103】
なお、図示した情報処理装置１の各構成要素は、必ずしも物理的に図示の如く構成されていることを要しない。すなわち、情報処理装置１の分散・統合の具体的態様は図示のものに限られず、その全部または一部を、各種の負荷や使用状況などに応じて、任意の単位で機能的または物理的に分散・統合して構成することができる。また、記憶部２０を情報処理装置１の外部装置としてネットワーク経由で接続するようにしても良い。
【0104】
また、上記実施例で説明した各種の処理は、予め用意されたプログラムをパーソナルコンピュータやワークステーションなどのコンピュータで実行することによって実現することができる。そこで、以下では、図１に示した情報処理装置１と同様の機能を実現する情報処理プログラムを実行するコンピュータの一例を説明する。ここでは、情報処理装置１と同様の機能を実現する情報処理プログラムを一例として説明する。図１３は、情報処理プログラムを実行するコンピュータの一例を示す図である。
【0105】
図１３に示すように、コンピュータ２００は、各種演算処理を実行するＣＰＵ２０３と、ユーザからのデータの入力を受け付ける入力装置２１５と、表示装置２０９を制御する表示制御部２０７とを有する。また、コンピュータ２００は、記憶媒体からプログラムなどを読取るドライブ装置２１３と、ネットワークを介して他のコンピュータとの間でデータの授受を行う通信制御部２１７とを有する。また、コンピュータ２００は、各種情報を一時記憶するメモリ２０１と、ＨＤＤ（Hard Disk Drive）２０５を有する。そして、メモリ２０１、ＣＰＵ２０３、ＨＤＤ２０５、表示制御部２０７、ドライブ装置２１３、入力装置２１５、通信制御部２１７は、バス２１９で接続されている。
【0106】
ドライブ装置２１３は、例えばリムーバブルディスク２１０用の装置である。ＨＤＤ２０５は、情報処理プログラム２０５ａおよび情報処理関連情報２０５ｂを記憶する。
【0107】
ＣＰＵ２０３は、情報処理プログラム２０５ａを読み出して、メモリ２０１に展開し、プロセスとして実行する。かかるプロセスは、情報処理装置１の各機能部に対応する。情報処理関連情報２０５ｂは、観測値リスト２１および条件リスト２２に対応する。そして、例えばリムーバブルディスク２１０が、情報処理プログラム２０５ａなどの各情報を記憶する。
【0108】
なお、情報処理プログラム２０５ａについては、必ずしも最初からＨＤＤ２０５に記憶させておかなくても良い。例えば、コンピュータ２００に挿入されるフレキシブルディスク（ＦＤ）、ＣＤ－ＲＯＭ（Compact Disk Read Only Memory）、ＤＶＤ（Digital Versatile Disk）、光磁気ディスク、ＩＣ（Integrated Circuit）カードなどの「可搬用の物理媒体」に当該プログラムを記憶させておく。そして、コンピュータ２００がこれらから情報処理プログラム２０５ａを読み出して実行するようにしても良い。
【符号の説明】
【0109】
１ 情報処理装置
１１ 平均値計算処理部
１２ 偏差計算処理部
１３ 相関係数計算部
１４ 判定部
２０ 記憶部
２１ 観測値リスト
２２ 条件リスト
１１ 平均値計算処理部
１２ 偏差計算処理部
１３ 相関係数計算部
１４ 判定部
２０ 記憶部
２１ 観測値リスト
２２ 条件リスト
【図1】
【図2】
【図3】
【図4A】
【図4B】
【図5】
【図6】
【図7A】
【図7B】
【図8A】
【図8B】
【図8C】
【図8D】
【図8E】
【図8F】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14A】
【図14B】
【図14C】
【図15】
【図16】
【図17】
【図18】