"ウェハのプロセス処理、段取り処理を実現するための構成を備えた半導体製造装置において、半導体製造装置の各部の駆動、プロセスを制御する制御装置、および半導体製造装置内の各物理量、または各化学反応量の状態をモニタするセンサと、製造中の半導体製造装置の前記制御装置、センサの出力信号を所定の間隔でサンプリングした装置データを記録するデータベースと、評価対象の半導体製造装置の装置データを前記データベースより検索して、解析対象の装置データを読み出し、比較対象の複数の信号データより時点間の相関行列を計算し、該相関行列より固有値、固有ベクトルを算出して、主成分分析により主成分スコアを算出する処理と、各主成分の固有値の大きさを比較して、降順に並べて一覧をユーザインタフェース画面に表示して、寄与率がある固有値を判定可能とする処理と、前記寄与率がある固有値に対応する主成分を選択した特徴空間に各信号の主成分スコアをプロットした散布図をユーザインタフェース画面に表示する処理とを実行する演算部とを備えたことを特徴とする半導体製造装置。前記演算部は、更に、前記寄与率がある固有値に対応する主成分ベクトルを、横軸をベクトルの要素番号順に処理時間として、縦軸をベクトルの値としたグラフでユーザインタフェース画面に表示する処理を実行することを特徴とする請求項1に記載の半導体製造装置。前記演算部は、更に、予め入力された異常検知判定基準を閾値として、前記信号間の主成分スコアの値の差が前記閾値を越えるものがある場合に異常と判定し、異常発生の警告を出力することを特徴とする請求項1に記載の半導体製造装置。前記演算部は、主成分スコアである変換値を用いて、各信号間の類似性により信号の変化を検知する際には、主成分分析の寄与率固有値の大きさの順に、指定された累積寄与率に達するまでの主成分数nを求め、第1主成分から第n主成分の主成分スコアの値を用いて類似性を評価する、ことを特徴とする請求項1に記載の半導体製造装置。前記演算部は、更に、前記寄与率がある固有値に対応する主成分iを順次選択して、該主成分iに対する全ての信号の主成分スコアxijより平均xiと標準偏差σiを算出し、次数式を満たさない信号を異常な信号と判定し、ここで、kは正数であり、異常と判定した信号を特定する情報を出力することを特徴とする請求項1に記載の半導体製造装置。ウェハのプロセス処理、段取り処理を実現するための構成を備えた半導体製造装置において、半導体製造装置の各部の駆動、プロセスを制御する制御装置、および半導体製造装置内の各物理量、または各化学反応量の状態をモニタするセンサと、製造中の半導体製造装置の前記制御装置、センサの出力信号を所定の間隔でサンプリングした装置データを記録するデータベースと、評価対象の半導体製造装置の装置データを前記データベースより検索して、解析対象の装置データを読み出し、比較対象の複数の信号データ間のそれぞれにおいて、サンプリング各時点に亘って信号強度の差の2乗の総和で信号間の距離と定義して、距離行列、内積行列を求め、該内積行列の固有値と固有ベクトルを算出し、多次元尺度構成法によりMDS写像における各信号の座標値を算出する処理と、各信号に対応する固有値の大きさを比較して、降順に並べて一覧をユーザインタフェース画面に表示して、寄与率がある固有値を判定可能とする処理と、前記寄与率がある固有値に対応するMDS写像座標軸を選択した特徴空間に各信号の座標値をプロットした散布図をユーザインタフェース画面に表示する処理とを実行する演算部とを備えたことを特徴とする半導体製造装置。前記演算部は、更に、予め入力された異常検知判定基準を閾値として、前記信号間のMDS写像座標値の差が前記閾値を越えるものがある場合に異常と判定し、異常発生の警告を出力することを特徴とする請求項6に記載の半導体製造装置。前記演算部は、座標値である変換値を用いて、各信号間の類似性により信号の変化を検知する際には、多次元尺度構成法の寄与率固有値の大きさの順に、指定された累積寄与率に達するまでの次元数nを求め、第1軸第1次元から第n軸第n次元の座標値を用いて類似性を評価する、ことを特徴とする請求項6に記載の半導体製造装置。前記演算部は、更に、前記寄与率がある固有値に対応するMDS写像の座標軸iを順次選択して、該MDS写像座標軸iに対する全ての信号のMDS写像値xijより平均xiと標準偏差σiを算出し、次数式を満たさない信号を異常な信号と判定し、ここで、kは正数であり、異常と判定した信号を特定する情報を出力することを特徴とする請求項6に記載の半導体製造装置。ウェハのプロセス処理、段取り処理を実現するための構成を備えた半導体製造装置において、半導体製造装置の各部の駆動、プロセスを制御する制御装置、および半導体製造装置内の各物理量、または各化学反応量の状態をモニタするセンサと、製造中の半導体製造装置の前記制御装置、センサの出力信号を所定の間隔でサンプリングした装置データを記録するデータベースと、評価対象の半導体製造装置の装置データを前記データベースより検索して、解析対象の装置データを読み出し、比較対象の複数の信号データよりサンプリング時点間で同一信号の差を取り2乗して、全信号について総和として時点間の距離を定義し、時点間の距離関係に基づき距離行列、内積行列を算出し、該内積行列の固有値と固有ベクトルを算出し、多次元尺度構成法により特徴空間の座標値行列と距離行列によりMDSスコアを算出する処理と、各時点に対応する固有値の大きさを比較して、降順に並べて一覧をユーザインタフェース画面に表示して、寄与率がある固有値を判定可能とする処理と、前記寄与率がある固有値に対応するMDSスコアを、横軸を固有ベクトルの要素番号順に処理時間として、縦軸をMDSスコアの値とした波形グラフでユーザインタフェース画面に表示する処理を実行する演算部とを備えたことを特徴とする半導体製造装置。前記演算部は、特徴空間の各軸に対する信号の特徴であるMDSスコアを用いて、各信号間の類似性により信号の変化を検知する際には、時点間の多次元尺度構成法の寄与率固有値の大きさの順に、指定された累積寄与率に達するまでの次元数nを求め、第1軸第1次元から第n軸次元に対する信号のMDSスコアを用いて類似性を評価する、ことを特徴とする請求項10に記載の半導体製造装置。前記演算部は、半導体製造装置内に備えられた第1の演算部と、ネットワークを介して接続された装置データ監視装置に備えられた第2の演算部に分けられ、前記データベースは、ネットワークを介して半導体製造装置と接続されるように設置され、第1の演算部は、半導体製造装置の前記制御装置、センサの出力信号を所定の間隔でサンプリングした装置データを、前記ネットワーク、前記第2の演算部を介して前記データベースへ記録を行い、第2の演算部は、評価対象の半導体製造装置の装置データを前記データベースより検索して、解析対象の装置データを読み出し、比較対象の複数の信号データの間の特徴的な変化の違いを定量値に変換して、変換値を用いて各信号間の類似性から異常と識別しうる信号を特定する処理を実行し、複数の信号間の定量値をグラフ表現した出力、および異常検知した警告出力を第1の演算部へ送信を行い、第1の演算部は、更に、第2の演算部から受信した前記グラフ情報、および異常検知警告情報をユーザインタフェース画面へ出力を行うことを特徴とする請求項1乃至請求項11のいずれかの請求項に記載の半導体製造装置。本発明は、製造処理中に、短いサンプリング間隔での装置データをモニタできる半導体製造装置において、装置データを解析することで、装置の異常発生を監視する方法に関する。監視対象の装置データには、流量、圧力、電流、電圧、温度といった、複数の信号間で相関関係が無いものも、化学反応中の信号値に変化があるOESデータといった、複数の信号間で相関関係があるものも含まれる。複数の信号とは、信号項目の異なる複数の信号でも、プロセスの処理の繰り返しにおいて同一信号項目の信号を過去から現在までの複数の信号であっても、解析の対象となる。また信号の異常を表す波形パターンの例としては、例えば図2に示したような、脈動、ハンチングといったノイズ、遅れやシフトといった変化であるが、さらにそれらに限定されるものではなく、事前に予測できない想定外の波形パターンも含まれる。このため実際にサンプリングした複数の信号を、信号の時系列での相関関係の有無に関わらず解析して、信号間の差異を表わすことが課題となる。特許文献1には、異常検知のための信号の処理として相関関係を利用した処理、また周波成分を取得する処理を利用することが記載される。これでは、例えば図2のシフトといった相関関係が無く、周波成分も無い信号の差異を検出できない。特許文献2では、相関の有無に関わらず、サンプリングした信号の時系列の強度変化から異常を検出できない。特許文献3で挙げられた各種の信号データ変換方法は、基本的には信号間の相関関係に基づく方法である。独立成分分析ICAは、信号を、白色雑音ではない信号の和に分解する方法、非負行列因子分解NMFは信号を非負の行列の積に分解する方法と記載される。これらの処理は音声信号処理、画像信号処理で用いられ、信号にノイズが混入したデータから特徴を抽出する解析法である。そのため、例えば図2に示した、複数の信号変化の違いを解析する方法ではない。特許文献4は、周波成分のある信号の処理に限定される。特許文献5は、相関関係のあるデータに限定される。また変化が無い信号の識別については、複数の信号の違いを解析するものではないので、図2のような信号の識別はできない。本発明は、複数の信号を対象に、信号間の時系列での相関関係の有無に関わらず、また事前に波形パターンを仮定せずに信号間の差異を表すことを目的とする。本発明によりどのような波形パターンの異常を示す信号変化があっても、異常を検知可能とする。また信号間の相関関係の検出を前提とする従来手法と異なり、相関関係を求めずに、信号間の微妙な変化の差異を検知可能とする。得られた複数の信号のみを利用して差異を表すことができるので、事前のパラメータ設定や波形パターン設定も不要とできるので、活用も容易となる。上記課題を解決するために本発明では、ウェハのプロセス処理、段取り処理を実現するための構成を備えた半導体製造装置を、半導体製造装置の各部の駆動、プロセスを制御する制御装置、および半導体製造装置内の各物理量、または各化学反応量の状態をモニタするセンサと、製造中の半導体製造装置の前記制御装置、センサの出力信号を所定の間隔でサンプリングした装置データを記録するデータベースと、評価対象の半導体製造装置の装置データを前記データベースより検索して、解析対象の装置データを読み出し、比較対象の複数の信号データより時点間の相関行列を計算し、該相関行列より固有値、固有ベクトルを算出して、主成分分析により主成分スコアを算出する処理と、各主成分の固有値の大きさを比較して、降順に並べて一覧をユーザインタフェース画面に表示して、寄与率がある固有値を判定可能とする処理と、前記寄与率がある固有値に対応する主成分を選択した特徴空間に各信号の主成分スコアをプロットした散布図をユーザインタフェース画面に表示する処理とを実行する演算部とを備えて構成した。また、上記課題を解決するために本発明では、ウェハのプロセス処理、段取り処理を実現するための構成を備えた半導体製造装置を、半導体製造装置の各部の駆動、プロセスを制御する制御装置、および半導体製造装置内の各物理量、または各化学反応量の状態をモニタするセンサと、製造中の半導体製造装置の前記制御装置、センサの出力信号を所定の間隔でサンプリングした装置データを記録するデータベースと、評価対象の半導体製造装置の装置データを前記データベースより検索して、解析対象の装置データを読み出し、比較対象の複数の信号データ間のそれぞれにおいて、サンプリング各時点に亘って信号強度の差の2乗の総和で信号間の距離と定義して、距離行列、内積行列を求め、該内積行列の固有値と固有ベクトルを算出し、多次元尺度構成法によりMDS写像における各信号の座標値を算出する処理と、各信号に対応する固有値の大きさを比較して、降順に並べて一覧をユーザインタフェース画面に表示して、寄与率がある固有値を判定可能とする処理と、前記寄与率がある固有値に対応するMDS写像座標軸を選択した特徴空間に各信号の座標値をプロットした散布図をユーザインタフェース画面に表示する処理とを実行する演算部とを備えて構成した。また、上記課題を解決するために本発明では、ウェハのプロセス処理、段取り処理を実現するための構成を備えた半導体製造装置を、半導体製造装置の各部の駆動、プロセスを制御する制御装置、および半導体製造装置内の各物理量、または各化学反応量の状態をモニタするセンサと、製造中の半導体製造装置の前記制御装置、センサの出力信号を所定の間隔でサンプリングした装置データを記録するデータベースと、評価対象の半導体製造装置の装置データを前記データベースより検索して、解析対象の装置データを読み出し、比較対象の複数の信号データよりサンプリング時点間で同一信号の差を取り2乗して、全信号について総和として時点間の距離を定義し、時点間の距離関係に基づき距離行列、内積行列を算出し、該内積行列の固有値と固有ベクトルを算出し、多次元尺度構成法により特徴空間の座標値行列と距離行列によりMDSスコアを算出する処理と、各時点に対応する固有値の大きさを比較して、降順に並べて一覧をユーザインタフェース画面に表示して、寄与率がある固有値を判定可能とする処理と、前記寄与率がある固有値に対応するMDSスコアを、横軸を固有ベクトルの要素番号順に処理時間として、縦軸をMDSスコアの値とした波形グラフでユーザインタフェース画面に表示する処理を実行する演算部とを備えて構成した。本発明によると、半導体製造装置や処理の異状検知、診断が容易となり、また対策の迅速化が図れる。特に、異常検知、診断に用いる、短いサンプリング時間間隔で取得した装置データはどのような信号変化でも良く、信号間の相関の有無、バラツキ、ノイズの有無と大小、信号強度差の大小に応じて計算方法を分ける必要が無いので、作業の手続きを簡便にできるとともに、本計算処理の装置組み込みと処理自動化も容易となる。信号の変化パターンが、例えばハンチングとシフトが混ざったような、複合的な変化であっても、解析処理後の特徴空間の軸に応じて差異を取りだせる。本発明の計算処理結果によれば、異常検知の自動化ができるとともに、複数の信号中の共通の信号変化と特殊な信号変化を特徴的な波形パターンとして抽出できるので、異常の調査、診断を効率的に進めることができる。これにより対策作業を迅速化、容易化できる。さらに本発明の計算結果によれば、処理中の装置データを特徴空間の座標値として、信号を特徴量のベクトルとして定量化できる。特に固有値の大きさによれば、信号の差異をあらわすための特徴量を絞れる。よって、処理毎に、ひとつ、もしくは少数の特徴量の変化を、処理の順にとっていくことで、処理による装置の経時変動をモニタできる。この情報を活用すればメンテナンスのタイミングを決定でき、さらには半導体ウェハプロセスで利用されるRun-to-Run制御法といったプロセス制御にも利用できる。添付の図面を参照しながら、本発明の実施形態について以下に説明する。本発明が対象とする半導体製造装置は、構成要素である各装置をモニタして、短いサンプリング間隔で信号を取得できる構成となっており、またサンプリングした装置データを記憶するデータベースと接続されていなければならない。半導体製造装置の構成要素は、装置を動作させる部分的な装置、部品であるが、OESといったプロセス状態をモニタするための付加的な装置を含んでいても良い。また取得した信号を解析するための計算機を備えている。主成分分析PCA複数の信号の違いは、相関関係の有無、バラツキの有無と大小、信号強度差の大小という違いである。それぞれ個別には、主成分分析PCAといった相関解析、フーリエ変換といった周波数解析、平均値の差の検定といった解析方法により差異を検出できるが、解析方法を適用するためには、信号に現れる波形パターンの特徴を事前に知っておかなければならない。これらの内、特に複数の信号を一括して計算処理し、自動的に差異を検出できる手法は主成分分析PCAに限られるが、信号間の時間軸に亘っての相関関係の有無を解析する方法であるため、バラツキや信号強度差を検出、識別できない。信号間に相関が無いということがわかるのみである。この主成分分析PCAは以下の式で求める相関行列Rの主成分を求めることで、信号の類似性を識別する方法である。相関行列Rは以下の式で求める。ここでxは信号データである。ボルド体xは信号のデータ数nを要素数とするベクトルであり、ntime行1列の行列とも解釈できる。変数上のバー-は平均を意味する。また添え字のi、jは各信号を意味する。右肩添え字の括弧の右肩添え字のTは行列の転置を表す。sは分散、共分散であり、rは相関係数である。主成分分析PCAは相関行列の固有値、固有ベクトルを求めることで、項目間の共変関係と、各関係の全体の変化に占める大きさを解析する方法である。元の各時点の信号データと固有ベクトルとの内積をとることで、その時点における主成分の大きさを、主成分スコアとして求めることができる。固有値、固有ベクトルと相関行列の関係は次となる。ここでλiは第i主成分の固有値、viが第i主成分ベクトルである。行列Vは主成分ベクトルviを列方向に並べた行列である。diagは引数を対角成分に並べて、それ以外の要素をゼロとした正方行列である。なお固有値は大きい順に、降順に並び、また主成分ベクトルも対応した並びとなっているとする。主成分スコアPcは次式で求める。Pcは行に時点毎、列に主成分毎の主成分スコアが並ぶ。Xは信号データxを列に並べた行列である。図3に示す信号に対して、この主成分分析を行った結果を説明する。図3には、処理ステップがStep1311とStep2312の2つのステップの連続処理である、4つの同じ処理の繰り返しの信号のデータを示す。処理においては、信号強度は処理条件の設定値となることが想定されるが、実際にモニタされた信号強度は、Step2312において、信号1303、信号2304、信号3305、信号4306に亘って、処理毎に縦軸でマイナス方向にシフトしている。また細かなバラツキを含む。これらの信号をStep1311、Step2312の、時点t0321から時点t2326の範囲で相関行列Rを求めると、rij1と全要素がほぼ1となる。主成分分析PCAによれば、全ての波形は同様に変化しているということがわかる。またStep1311、もしくはStep2312のみの範囲で相関行列Rを求めると、rii=1、rij0と、ほぼ単位行列となり、主成分分析PCAによれば、各波形はそれぞれ独立であるということがわかる。バラツキの大きさが変わったとしても、そのバラツキ成分の平均がゼロで、信号間で独立であれば、相関行列Rはほぼ単位行列となり、同様の結果となる。以上より複数の信号のばらつきや信号強度差を一括して検出、識別できないことがわかる。時点間主成分分析IT-PCA複数の信号を一括処理し、差異を検出、識別するために、主成分分析PCAで解析するならば、時間軸に沿った信号間の相関関係ではなく、他の評価軸に沿った信号間の相関関係を用いなければならない。もしくは、バラツキであるにしろ平均の差であるにしろ、信号間の強度差を直接に一括処理し、差異を検出、識別できる方法を利用しなければならない。図3のStep1311、Step2312とも、それぞれの範囲で信号間に相関は無い。信号毎に信号強度差はある。そこで横軸に信号、縦軸に信号強度をとったグラフを図4に示す。プロットは図3における時点t11323、t12324、t13325の信号強度である。各時点の信号強度を信号に亘って線分で結べば、時点t11、t12、t13は、それぞれ折れ線411、412、413となる。このように時点間では信号強度に相関がある。信号1421、信号2422、信号3423、信号4424に沿って負の相関である。図3のStep1311の範囲に着目すると、信号強度は昇順で、信号1303、信号2304、信号3305、信号4306の順に大きくなっているので、やはり時点間で信号強度に相関がある。このように時点間で相関を取れば、信号強度差に相関関係を求めることができる。また特定の信号に図2に示したようなハンチングがあったとしても、ハンチングのある時点で、信号間で信号強度差があるため、時点間で信号強度に相関がでる。時点間で相関関係が無い複数の信号とは、複数の信号が一致している場合、もしくは複数の信号について各時点での平均との差が正負に均等にばらついて、時点間の積和が相殺されてゼロになる場合である。時点間の相関行列Rを求める計算は、数式1、数式2を用いればよい。この際、xは各時点における信号データであり、ボルド体のxは信号の数nsignalを要素数とするベクトルである。また添え字のi、jは各時点を意味する。なお、複数の信号について時点間で相関関係が無くなる場合とは、数式1がゼロとなることである。主成分を軸とした特徴空間に、主成分スコアPcを配置すれば、信号の差異を検出、識別できる。強度変化が類似した複数の信号は近い位置に配置され、強度変化の差異が大きい信号は離れた位置に配置される。また主成分ベクトルには特徴的な信号が現れる。例えば複数の信号で共通した変化、特定の信号にしか存在しない変化が現れ、またこれらの変化は、主成分の固有値の大きさ寄与に応じた順番の主成分ベクトルに現れる。この時点間の相関行列を用いた主成分分析PCAを、時点間主成分分析IT-PCAと呼ぶ。主成分分析PCAの計算方法自体は数式3から数式5と同じである。ひとつの信号について、時点数ntimeの主成分スコアが求まる。主成分スコアを配置した空間において、各信号の差異を検出、識別できる。多次元尺度構成法MDS信号間の強度差を直接に一括処理し、差異を検出、識別するためには、信号間の非類似性dissimilarity、もしくは距離を定義できれば、多次元尺度構成法MDSにより信号間の関係を特徴量による空間に配置して、差異を検出、識別できる。多次元尺度構成法MDSは、空間に配置されている各サンプル間の距離のデータから、サンプルの座標値を求める方法である。今、ndimension次元空間中にnsample個のサンプルが配置されているとする。サンプル間の距離がユークリッド距離で定まる場合には、次式がサンプル間の距離となる。ここでi、jは各サンプルのインデクスであり、kは次元のインデクスである。nsample×nsampleの距離行列Tは次式となる。なお距離行列Tを、非類似性行列とも呼ぶ。中心化行列Gnを用いて、内積行列Bcを求める。ここでInはn×nの単位行列である。1nは全ての要素を1とする要素数nのベクトルである。サンプル数はnsampleであるので、n=nsampleである。内積行列Bcを固有値分解する。Λは対角要素を固有地、他要素をゼロとしたnsample×nsampleの行列、Vは固有ベクトルを列に並べたnsample×nsampleの行列である。内積行列Bcとndimension次元空間中の各サンプルの座標値の行列Xcとの間にはBc=XcXcTの関係がある。よって座標値の行列Xcは次式で得られる。2次元空間中のサンプルである場合、行列Xcの1列目と2列目が第1軸、第2軸の座標値となる。つまり数式11の座標値xcijのインデクスiはサンプル、jは多次元尺度構成法MDSの写像である特徴空間の軸を意味する。XcをMDS写像と呼ぶ。MDS写像の座標値はnsample数分、算出される。信号間多次元尺度構成法IS-MDS複数の信号の強度差を用いて、差異を検出、識別するために、信号間に距離を定義する。各時点に亘って、信号強度の差の2乗の総和で距離を定義できる。つまり式のi、jを信号についてのインデクス、kを時点についてのインデクスとすることである。ndimension=ntime、nsample=nsignalである。この距離の定義により多次元尺度構成法MDSで解析すれば、数式11で算出される特徴空間の座標値で、信号の差異が表現される。この距離の定義に基づき信号を解析する多次元尺度構成法MDSを信号間多次元尺度構成法IS-MDSと呼ぶ。時点間多次元尺度構成法IT-MDS距離の定義について、時点間主成分分析IT-PCAと同様に時点間の信号強度の差に着目すれば、時点間での信号強度の差を距離とすることもできる。すなわち、時点間で同一信号の差を取り2乗して、全信号について総和として時点間の距離を定義する。図5に、図3の信号について時点t0321とt2326の信号強度の、信号に亘るグラフを示す。図5中の両矢印の差を取り2乗して、全信号について総和として時点間の距離とする。数式6のi、jを時点についてのインデクス、kを信号についてのインデクスとすることである。ndimension=nsignal、nsample=ntimeである。この距離の定義により多次元尺度構成法で解析すれば、数式11で算出される特徴空間の座標値は、信号に含まれる特徴的な信号の波形となる。つまり座標値行列Xcの行を横軸として、またデータ値を縦軸として、各列のデータ値を波形としてグラフ化すれば、グラフには特徴的な信号が示される。この信号の波形を、実際の信号の差異の大きさを反映して表すためには、座標値行列Xcに距離行列Tを掛け合わせる。これを多次元尺度構成法MDSスコアMDSscoreと本明細書では呼ぶ。MDSスコアMDSscoreはntime×ntimeの行列である。列毎に第i座標軸に対応した、要素数ntimeのMDSスコアmdsscoreiとなる。時点間での信号強度の差による距離の定義に基づき、信号を解析する多次元尺度構成法MDSを時点間多次元尺度構成法IT-MDSと呼ぶ。以上の方法によれば、処理時間に亘って、信号間に相関関係が無くても、バラツキの有無と大小、信号強度差の大小という違いがあれば、複数の信号を一括して計算し、特徴空間において、信号の差異を検出、識別可能となる。また事前に信号の変化パターンを用意する必要もない。半導体製造装置の構成例図6および図7により、本発明の一実施の形態に係る半導体製造装置の構成および動作について説明する。図6は本発明の一実施の形態に係る半導体製造装置の構成例を示す構成図、図7は本発明の一実施の形態に係る半導体製造装置601に備えられた計算機111において実行される、またはネットワーク131を介して接続される装置データ監視装置133において実行される半導体製造装置監視処理演算部の機能構成を示す構成図である。図6に示す通り、半導体製造装置601は、処理室602、プロセス処理用の材料の供給系603、プロセス処理後の廃材の排出系604、装置コントローラ、外部通信装置109、計算機システムである計算機、記憶装置111、端末である画面、ユーザーインターフェース112、所定の処理を実現するための装置の各部位を動作変化させる装置制御装置1611から装置制御装置n615、プロセス処理の状況を検出するための処理中センサ1621から処理中センサm623で構成される。処理室602とは、図1に示したエッチング装置101の例においては、チャンバ102のことである。薄膜を堆積するプラズマCVD装置やPVD装置でも処理室のことをチャンバと呼ぶ。熱酸化装置、LPCVD装置では、複数のウェハをボートに装填して処理を行い、この処理室のことを炉furnaceと呼ぶ。洗浄装置やめっき装置では、特に薬液にウェハを入れるので、処理室は浴槽となる。CMP装置は、研磨ヘッドでウェハを持ち、回転定盤上のパッドに押し付けて、回転運動することで研磨を行う。この研磨のための装置構成部位が処理室である。ウェハのプロセス処理を直接に行う部位が処理室である。供給系603、排出系604は、図1に示した例においては、ガス供給系108、排気系107となる。他の種類の装置においても、プロセス処理実現のための材料を供給し、処理後の材料を排出する構成を有する。装置制御装置1611から装置制御装置n615は、図1に示した例においては、流量調整装置113、圧力調整装置114、電力調整装置115、温度調整装置116に該当する。他の種類の装置においても、装置を動作させ、プロセスを調整するための装置を備える。また半導体製造装置は、ウェハを処理室に搬送する構造を備えており、この搬送装置を動作させるための制御装置も、装置制御装置として扱ってよい。処理中センサ1621から処理中センサm623は、図1に示した例においては、OESである。各装置制御装置611-615、各処理中センサ621-623の出力は、図示してはいないが装置コントローラ、外部通信装置109を経由して、計算機、記憶装置111へ入力される。また、製造ラインに設置された複数の半導体製造装置を総合的に監視する装置データ監視装置133が設置されている場合には、装置コントローラ、外部通信装置109より各装置制御装置611-615、各処理中センサ621-623の出力をネットワーク131を介して、装置データ監視装置133へ送信する。装置データ監視装置133により、装置データDB132へ記録される。半導体製造装置601はネットワーク131を介して、装置データDB132と接続されており、またデータ共用の便を図り、計算機システムである装置データ監視装置133もネットワーク131に接続されている。計算機、記憶装置111は、製造処理中においては、前記各装置制御装置、各処理中センサの出力を取得して、装置データとして装置コントローラ、外部通信装置109を介して装置データDB132へ記録し、半導体製造装置の監視処理においては、前記装置データDB132より装置コントローラ、外部通信装置109を介して該当する装置データを読み出して、解析を行うことで装置監視処理を実行する。また計算機、記憶装置111は、画面、ユーザーインターフェース112と接続されており、装置監視の状況や異常発生の警告、解析結果である主成分スコア、MDS写像のグラフや、信号波形のグラフ、固有値一覧といった情報をユーザーに表示する。これらの処理の中で前記装置データDB132より該当する装置データを読み出して、解析を行う装置監視処理は、ネットワーク131を介して接続された装置データ監視装置133が処理しても良い。もしくは、装置データ監視装置133の機能は、半導体製造装置601の計算機、記憶装置111により実行されるソフトウェアとして、および装置データDB132が計算機、記憶装置111内に構成されて半導体製造装置601に備えられているという構成でも良く、この場合、半導体製造装置601をネットワークに接続しなくても監視処理システムは装置を監視できる。半導体製造装置監視処理演算部の機能構成計算機、記憶装置111において実行する半導体製造装置監視処理演算部701は、ロット、ウェハ毎処理来歴取得部702、装置データ取得部703、相関行列算出部704、距離行列、内積行列算出部705、固有値計算部706、主成分スコア算出部707、MDS写像算出部708、MDSスコア算出部709、固有値一覧表示部710、特徴空間表示部711、特徴信号波形表示部712、異常検知判定基準取得部713、異常検知部714とを有する。ロット、ウェハ毎処理来歴取得部702は、各装置制御装置611-615、各処理中センサ621-623からの出力信号を、所定のサンプリング周期において取得して、それらの装置データを装置データDB132に記憶する。装置データDB132のデータ構成装置データDB132のデータ構成例を図29に示す。装置データDB132をネットワークに接続して、複数の半導体製造装置の装置データを記憶する共通のデータベースとして構成する場合には、記憶する装置データを分類するためのディレクトリとして、半導体製造装置ID2901を有する。さらに、当該半導体装置において処理されたロット又はウェハID2902、製造工程ID2903、製造レシピID2904、製造日時2905のデータレコードより成るディレクトリによって、装置データを分類している。そして、当該半導体製造装置において、当該ロットまたはウェハの製造処理が実施されている間、信号ごとに所定のサンプリン周期で、装置データを格納する。装置データは、信号ID2906、観測開始日時2907、観測終了日時2908、サンプリング間隔2909、時系列装置データ2910のデータ項目より成るデータレコード形式に記憶される。例えば、信号IDがsig-1の信号データは観測開始日時が10月18日09時20分10秒00からサンプリングを開始して、観測終了日時が10月18日09時26分25秒50まで、サンプリング間隔が0.50秒で取得された信号強度データを、時系列装置データの欄に格納する。半導体製造装置601をネットワーク131に接続して、ネットワーク上の装置データ監視装置より装置監視処理を実行しない場合には、上記した装置データDB132を計算機、記憶装置111内に構成する。半導体製造装置601の監視処理の概要半導体製造装置601の監視処理は、例えば同一工程の、同一レシピの複数回の製造処理履歴を検索して、該当装置データの実値を取得して解析を行うことで、半導体製造装置状態を監視し、異常を検知する。なおこの処理は計算機記憶装置111が行う場合と、ネットワーク上の装置データ監視装置133が行う場合とが考えられる。まず解析の対象とする装置データを特定し、装置データ取得部703により、装置データDB132から装置データを取得する。取得した装置データを対象として、時点間主成分分析IT-PCA、信号間多次元尺度構成法IS-MDS、または時点間多次元尺度構成法IT-MDSにより解析処理を行う。時点間主成分分析IT-PCAでは、まず相関行列算出部704で時点間の相関行列を計算する。そして固有値計算部706で固有値と固有ベクトルを算出し、主成分スコア算出部707で主成分スコアを算出する。固有値一覧表示部710で固有値の大きさを比較表示する。これにより信号中の共通の変化パターン、もしくは特異な変化パターンがあるかを、固有値の大きさにより判別可能とする。異常検知判定基準取得処理713で信号が異常であると判定する基準をユーザ入力より取得する。この基準とは、信号間の主成分スコアの値の差の大きさについての閾値である。閾値を超えた場合に異常であるとする。異常検知部714で異常の判定をし、異常発生と判定すれば、画面ユーザーインターフェース112に警告を表示する。計算処理の結果については、特徴空間表示部711で画面、ユーザーインターフェース112に主成分スコアのグラフを表示し、特徴信号波形表示部712で、主成分ベクトルを、横軸をベクトルの要素番号順に処理時間とし、縦軸をベクトルの値としたグラフで表示する。また固有値一覧も表示する。これによりユーザは発生した信号の、過去の信号との差異や波形の特徴を確認できる。これらのユーザに対する出力においては、ネットワークを介して外部に情報を通知しても良く、電子メールや音声処理など出力の形態は様々とできる。信号間多次元尺度構成法IS-MDSでは、距離行列、内積行列算出部705で波形間の距離を計算し、内積行列を求める。固有値計算部706で内積行列の固有値と固有ベクトルを求め、MDS写像算出部708でMDS写像の座標値を求める。異常検知判定基準取得部713で異常発生の判定のための基準を取得し、異常検知部714で異常を判定し、警告を出力する。特徴空間表示部711で画面、ユーザーインターフェース112にMDS写像のグラフを表示する。また固有値一覧表示部710で固有値一覧も表示する。時点間多次元尺度構成法IT-MDSでは、距離行列、内積行列算出部705で時点間の距離を計算し、内積行列を求める。固有値計算部706で内積行列の固有値と固有ベクトルを求め、MDS写像算出部708でMDS写像の座標値行列Xcを得る。MDSスコア算出部709で座標値行列Xcと距離行列TよりMDSスコアを求める。MDSスコアは処理時間に対する信号の特徴的な波形パターンである。特徴信号波形表示部712で画面、ユーザーインターフェース112にMDSスコアの処理時間に対するグラフを表示する。また固有値一覧表示部710で固有値一覧も表示する。以上のシステムの構成により、複数の装置データより、信号の違いを解析して、自動的に異常を検知可能となるとともに、解析結果を表示することによりユーザが診断や対策を進めていくことができる。まず半導体製造装置における、本発明での半導体製造装置の監視方法の概要を、図8を用いて説明する。半導体製造装置601の構成は図6に示したものと同一である。装置監視処理を実行するのは、本実施例では装置の計算機、記憶装置111が行うものとしている。半導体製造装置601で、複数回のウェハの処理が実施されて、処理中の各装置制御装置611、各処理中センサ621の出力信号は、装置コントローラ、外部通信装置109を介して計算機、記憶装置111に入力されて装置データDB132に記憶される。記憶された装置データである各種信号を、処理回数の分、すなわち各種信号×処理回数801のデータを、装置データDB132から計算機、記憶装置111が取得し、解析処理を実行する。本実施例では、例えばエッチング装置で処理が行われた12回の処理の信号を対象とする。信号変化はグラフ802に示される。12回の処理で信号強度はシフトし、7回目のひとつの信号でハンチングが発生している。これらの信号の、サンプリング時点間の相関803を算出し、主成分分析PCA804を実行すると、特徴量写像807が得られる。またサンプリング信号間の差、時点間の差805を算出し、多次元尺度構成法MDS806を実行しても、特徴量写像807が得られる。特徴量写像807においては、各信号は、シフトの大きさの変化はその大きさに応じてpc1xc1軸方向に並び、またpc2xc2軸方向にハンチングがある信号が分離される。以上の処理は自動的に計算機で実行でき、これらの特徴量の差をとることで、異常といった信号間の違いを自動で検知できる。また主成分分析PCA804ならば主成分ベクトル、また時点間の差を用いた多次元尺度構成法MDS806ならば特徴量で、処理時間に亘っての特徴的な信号変化のデータを得ることができる。これらの計算結果を、グラフとして画面、ユーザーインターフェース112に表示すれば、ユーザは異常の有無や信号変化の特徴を確認できる。図9に示す前記した12回の処理の信号を対象に、各信号を特徴空間に配置し、また信号に含まれる特長的な波形パターンを抽出する方法を以下に示す。グラフの縦軸Signal902の単位系A.U.は任意単位系を意味する。図9では、各信号はプロセス処理の毎に信号001911から012922に亘って信号強度がシフトしていき、また信号007917にはハンチングが入っている。まず時点間主成分分析IT-PCAによる方法を示す。相関行列算出部704は時点間の信号の相関行列Rを算出する。信号のデータを、行に各信号、列に各時点を取り、列間の平均の縦ベクトルを列ごとに差し引いた行列を今Xtmpとおくならば、として相関行列は算出される。そして固有値計算部706が固有値λiと固有ベクトルviを計算し、主成分スコア算出部707が主成分スコアpciを算出する。iは1からntimeまでの整数であり、固有ベクトルviは要素数ntimeのベクトル、pciは要素数nsignalのベクトルとなる。pciの要素の順番は、信号のデータの各行の信号に対応する。図10に、固有値一覧表示部710により表示される固有値、すなわち寄与率CR1002を大きさの順に、降順に示した表を示す。CCR1003は累積寄与率である。これによれば第1主成分の寄与が95とほとんどを占め、第2主成分の寄与が5、それ以降は寄与が無いことがわかる。なお、固有値、固有ベクトルはntime個、算出される。図11に特徴空間表示部711による特徴空間における各信号の主成分スコアのプロットを示す。は第1主成分と第2主成分に対する散布図であり、は第3主成分と第4主成分に対する散布図である。これよりpc11101で信号の強度変化が表現され、またpc21102で信号007のハンチングによる違いが表現されているとわかる。図では各信号がばらついているだけであり、固有値の情報の寄与率が0であることから、これらの主成分は無意味であることがわかる。よって第1主成分と第2主成分の値を信号間で比較し、差を判定すれば自動的に異常が発生したか判別できることとなる。図12に、特徴信号波形表示部712が算出した特徴信号波形である、処理時間に亘っての主成分ベクトルを示す。図12固有ベクトル1では、平均値は0.07程度となっており、すなわち信号強度に信号間で差があることがわかる。また図12固有ベクトル2より信号007に含まれるハンチングが特徴信号波形として得られていることがわかる。図12との主成分ベクトルはランダムにばらついているだけであり、固有値の情報の寄与率が0であることから、これらの主成分は無意味であることがわかる。次に信号間多次元尺度構成法IS-MDSによる方法を示す。距離行列、内積行列算出部705は、信号間の時点毎の信号強度の差に基づき、数式7の距離行列Tを算出する。数式8から数式10の計算により内積行列Bcを算出してから、固有値計算部706が固有値λi、固有ベクトルviを計算する。MDS写像算出部708が数式11の計算で特徴空間であるMDS写像における各信号の座標値xciを算出する。iは1からnsignalまでの整数であり、固有ベクトルviは要素数nsignalのベクトルである。iの順番は、距離行列Tの行、列に対応する信号の順番に対応する。なお、固有値、固有ベクトルはnsignal個、算出される。図13が算出した距離行列Tである。対角成分はゼロとなる。また処理の繰り返し毎の信号001から012に向かって第1行、第1列ともに距離が増加していくが、信号007で距離が増加しており、ハンチングによる影響が入っていることがわかる。図14に信号間多次元尺度構成法IS-MDSで算出した、固有値の一覧を示す。第1主成分が90とほとんどを占め、第2主成分が10、それ以降は0で寄与は無い。図15に、MDS写像算出部708、特徴空間表示部711が算出した特徴空間である、各信号のMDS写像におけるプロットを示す。図15が第1軸xc1、第2軸xc2に対するプロットであり、これより第1軸xc11501で信号強度差、第2軸xc21502で信号007のハンチングによる違いが表現されているとわかる。図15は第3軸xc3、第4軸xc4に対するプロットであり、固有値の情報の寄与率が0であることから、これらの軸は無意味であることがわかる。時点間多次元尺度構成法IT-MDSによる方法を示す。この方法は信号中の特徴信号波形を抽出するための方法である。距離行列、内積行列算出部705は、各信号における時点間の差を取り、数式7の距離行列Tを算出する。信号間多次元尺度構成法IS-MDSと同様に内積行列Bcを算出し、固有値計算部706が固有値λi、固有ベクトルviを計算する。iは1からntimeまでの整数であり、固有ベクトルviは要素数ntimeのベクトルである。iの順番は、距離行列Tの行、列に対応し、すなわち処理時間に亘る並びである。そして、MDSスコア算出部709が数式12でMDSスコアmdsscoreiを求める。mdsscoreiはMDSscoreのi列をとったベクトルである。図16に固有値一覧表示部710が画面、ユーザーインターフェース112に表示した固有値の一覧を示す。第1固有値の寄与が95とほとんどを示す。累積寄与率CRR1603を見ることで、第10固有値までは小数点2桁の精度で寄与率があり、それ以降は寄与率0である。なお、固有値、固有ベクトルはntime個、算出される。図17に、特徴信号波形表示部712が画面、ユーザーインターフェース112に表示した処理時間に亘るMDSスコアのグラフを示す。第1から第12までのMDSスコアを示す。図17はMDSスコアの全波形を示したグラフであり、図17はMDSスコアの値の範囲をゼロの付近に拡大したグラフである。第1MDSスコアscore01によりハンチングと信号強度シフトの特徴があることがわかる。また第2から第11までのMDSスコアには信号強度シフトの特徴があることがわかる。第12以降のMDSスコアはゼロとなる。第1MDSスコアにハンチングと信号強度シフトの特徴があるので、このMDSスコアだけで、複数の信号に含まれる特徴的な波形パターンを確認できる。信号間多次元尺度構成法IS-MDSと時点間多次元尺度構成法IT-MDSによれば、MDS写像による特徴空間で異常を自動判定でき、また信号に含まれる特徴信号波形を抽出できることがわかる。図18に示す12回の処理の信号を対象に、各信号を特徴空間に配置し、また信号に含まれる特長的な波形パターンを抽出する方法を示す。各信号はプロセス処理の毎に信号強度がプラス方向にシフトしていき、さらに信号0011811、0041814、0071817、0101820には、それぞれ15秒、35秒、55秒、75秒近傍でハンチングが発生している。また0011811、0041814、0071817、0101820は0秒の解析対象の開始時点から立上り遅れが発生している。0011811、0041814、0071817、0101820の順に遅れは大きくなっている。この遅れのため4つの信号には信号間で、相関関係もあることになる。本実施例では、時点に亘る信号間の相関関係の有無に関わらず、本発明の方法である時点間主成分分析IT-PCA、信号間多次元尺度構成法IS-MDS、時点間多次元尺度構成法IT-MDSにより、特徴空間において異常を自動判定でき、また信号に含まれる特徴信号波形を抽出できることを示す。時点間主成分分析IT-PCAによる、複数信号の特徴空間への変換方法と波形パターンの抽出方法を示す。計算方法は実施例1に示した方法と同じである。図19に寄与率CR1902を降順に示した表を、累積寄与率CCR1903とともに示す。第1主成分の寄与率が80と大きいが、第2主成分は12、また第3主成分から第5主成分までも数の寄与率となっている。すなわち第5主成分までは寄与がある信号中の波形パターンであることがわかる。それ以降は寄与が無いこともわかる。図20に特徴空間表示部711が画面、ユーザーインターフェース112に表示した各信号の主成分スコアのプロットを示す。は第1主成分と第2主成分の散布図であり、は第3主成分と第4主成分の散布図、は第5主成分と第6主成分の散布図である。寄与率のある第5主成分までが入るようにグラフを示した。これらより、まずpc1でシフトによる強度変化が表現され、pc2には遅れとハンチングによる違いが表現されていることがわかる。特にpc2の負の方向に、001、004、007、010が並んでいるため、特に遅れの違いをよく反映していると考えられる。pc3ではハンチングの有無が表現され、pc4では特に004、007の違いが出ている。またpc5ではシフトだけの信号に対し、正負にハンチングのある信号が散っている。第6主成分の寄与率はゼロなので、pc6の信号の配置は無意味である。よって第1主成分から第5主成分までの主成分の値を信号間で比較し、差を判定すれば自動的に異常が発生したか判別できる。図21に、特徴信号波形表示部712画面、ユーザーインターフェース112に表示した処理時間に亘っての主成分ベクトルを示す。図21からは、第1主成分から第6主成分の固有ベクトルをそれぞれ示している。固有ベクトル1では平均がゼロではなく、およそ0.06程度となっているので、信号強度に信号間に差があることがわかる。固有ベクトル2には、ハンチングの波形が入り込んでいるが、立上り遅れの波形変化が強く現れて、また010のハンチングの強度が強く現れている。固有ベクトル3は001のハンチングが、固有ベクトル4は004のハンチング、固有ベクトル5には007のハンチングが強く現れている。寄与がゼロの固有ベクトル6はランダムにばらついている。固有ベクトル1から5までを観察することで、信号の特徴的な波形パターンがわかった。次に信号間多次元尺度構成法IS-MDSによる複数信号の特徴空間への変換方法を示す。計算機における計算内容は実施例1と同じである。図22に信号間多次元尺度構成法IS-MDSで算出した固有値の一覧を示す。第1主成分から第5主成分まで順に60、30、5、3、2である。以降で寄与は無い。図23に特徴空間であるMDS写像を示す。図23が第1軸xc1、第2軸xc2、図23第3軸xc3、第4軸xc4、図23が第5軸xc5、第6軸xc6に対するプロットである。第1軸についてはxc1=2にハンチングと立上り遅れの信号が並んでいるが、他の信号は負の方向に並んでいる。また第2軸については正の方向に信号の強度のシフトが現れる。但し、第2軸でハンチングと立上り遅れのある信号と無い信号とではxc2の値が対応していない。以上から第1軸、第2軸には、いずれにも信号の強度のシフト、立上り遅れとハンチングの影響が含まれていることとなる。第3軸、第4軸は信号001、004、010の信号の違いが現れ、また第5軸は主に信号007の違いが現れる。立上り遅れとハンチングのある信号と無い信号との、MDS写像における差、もしくは距離を判定すれば自動的に異常が発生したか判別できる。時点間多次元尺度構成法IT-MDSによる、複数の信号に含まれる波形パターンの抽出方法を示す。計算方法は実施例1に示した方法と同じである。図24に固有値の一覧を示す。第1固有値の寄与が83と大きく、第2固有値は10、第3固有値、第4固有値はそれぞれ5、3である。以降、寄与は無い。図25に、処理時間に亘るMDSスコアのグラフを示す。第1から第13までのMDSスコアを示す。図25はMDSスコアの全波形を示したグラフであり、図25は縦軸Score2512を-4,12としたグラフ、図25は縦軸Score2522を-0.04,0.04としたグラフである。図25ではscore01を、また図25ではscore01、02、03、04を描画していない。第1MDSスコアscore01には、強度シフトと立上り遅れ、ハンチングの特徴が現れる。第2MDSスコアscore02、第4MDSスコアscore04には主に強度シフトとハンチングが現れる。第3MDSスコアscore03にはハンチングが現れる。第5MDSスコアscore05から第12MDSスコアscore12は、非常に小さなMDSスコアの値であるが、強度シフトの特徴が現れる。寄与率の大きい第1MDSスコアだけでも立上り遅れ、ハンチング、強度シフトの特徴的な波形パターンを確認できる。もしくは寄与率がゼロではない第4MDSスコアまでをグラフとして、特徴的な波形パターンを確認しても良い。本実施例においても、信号間多次元尺度構成法IS-MDSと時点間多次元尺度構成法IT-MDSによれば、MDS写像による特徴空間で異常を自動判定でき、また信号に含まれる特徴信号波形を抽出できることがわかる。本実施例では、異常検知する方法、及び解析結果の画面、ユーザーインターフェースへの表示の例を示す。時点間主成分分析IT-PCAでは主成分スコア、信号間多次元尺度構成法IS-MDSではMDS写像を特徴空間とし、特徴空間における複数信号の配置関係を判定することにより、異常な信号を識別することが、異常検知処理である。特徴空間における複数信号の配置関係とは、例えば図11、図20に示した主成分スコアのプロット、また図15、図23に示したMDS写像のプロットである。ここで、一つの主成分スコアに対する異常の判定方法を説明する。図26に実施例1の時点間主成分分析IT-PCAにより計算された第2主成分pc2の信号の分布を示す。また図26に実施例2の時点間主成分分析IT-PCAにより計算された第3主成分の信号の分布を示す。図26の分布において、全ての信号の平均average2603の1σ2602と-1σ2604の範囲に信号007を除く全ての信号が含まれる。つまり信号007が他の信号に対して特異であることになる。そこで異常が発生したと判定するためには、数式14が満たされるならば異常無しとして、異常な信号を求めればよい。ここで、xは主成分スコア、もしくはMDS写像値であり、iは主成分番号、もしくはMDS写像の軸番号である。jは信号を意味するインデクスである。変数の上のバー-は平均を意味する。またσは標準偏差である。範囲を1σに限定する必要は無く、一般的には正数kにより範囲を定めておけばよい。図26においては、数式14によれば、信号001、004、007が異常と判定される。そこで図20を確認すれば、第2主成分pc22002において信号010が最も平均から解離していることがわかる。このように主成分、もしくはMDS写像の軸を変えれば異常な信号を判定できる。なお評価すべき軸は、計算で得られた固有値で判断すればよく、例えば実施例1の時点間主成分分析IT-PCAならば図10より寄与のある第2主成分までを評価すればよいことがわかる。実施例2の時点間主成分分析IT-PCAならば図19より第5主成分までを評価すればよい。異常の評価対象とする主成分数、もしくはMDS写像の軸数をlまでとし、多変量の主成分スコア、もしくはMDS写像値のベクトルにより、複数の主成分、もしくは複数のMDS写像の軸において、異常を判定することもできる。この場合、数式16を満たす場合に異常無しとして、異常な信号を求めればよい。ここでxは要素数がlの主成分スコア、もしくはMDS写像値のベクトルである。ベクトルの上のバー-は要素毎に信号間で平均を取ったベクトルであることを意味する。σは主成分、もしくはMDS写像の軸毎の、信号間の標準偏差のベクトルである。Sは信号の主成分スコア、もしくはMDS写像値の、l×lの標本分散共分散行列である。特徴空間において信号の異常を判定する方法として、予め異常と正常の判別の基準を過去の事例から決めることができるならば、判別分析やサポートベクターマシンSVMを用いてもよい。また特徴空間における信号の異常を判定するために、グループ分類を利用する方法もある。図27に図26に示した時点間主成分分析IT-PCAの第3主成分を用いた、クラスター分析による、信号の樹形図を示す。これによれば第1の分岐1stdivision2701により、Group12702として信号001、004、007、010が、他の信号のグループGroup22703と分離される。信号の数を比較し、信号の数が少ないグループを異常とすればよい。このように上位階層2704での分岐によるグループ分けで異常を判定できる。グループ分類の方法としては、クラスター分析の他にk-means法も利用できる。次に分析結果の画面、ユーザーインターフェースへの表示内容を示す。この表示は、異常の有無に関わらずに処理しても良く、プロセス処理の度に表示しても良い。また異常の発生は、必ずしも画面でユーザーが確認しなければならないものではなく、メールやサイレンなどの音で警告を出しても良い。画面表示の例を図28に示す。画面2801には、異常判定結果2802、固有値による寄与率一覧2803、特徴空間における信号のプロット2805と信号のプロットの凡例2804、そして信号の特徴波形2806が表示される。これらの情報によりユーザーは異常の発生とそのときの信号に含まれる波形の特徴を確認できる。これらの情報は、時点間主成分分析IT-PCAならば、固有値と主成分スコア、主成分ベクトルにより、それぞれ寄与率一覧2803、特徴空間における信号のプロット2805、信号の特徴波形2806を表示できる。多次元尺度構成法で分析を行うならば、信号間多次元尺度構成法IS-MDSによる固有値、MDS写像で、それぞれ寄与率一覧2803、特徴空間における信号のプロット2805を表示し、時点間多次元尺度構成法IT-MDSによるMDSスコアで信号の特徴波形2806を表示できる。時点間多次元尺度構成法IT-MDSによる固有値も寄与率一覧2803に含めても良い。または、時点間主成分分析IT-PCA、信号間多次元尺度構成法IS-MDS、時点間多次元尺度構成法IT-MDSの全ての情報を画面に表示しても良い。各情報の表示の形式は図28の各画面を要素のものとすればよい。以上が本発明による監視方法を備えた半導体製造装置の処理内容である。なお、本発明は半導体製造装置に関するものであるが、時点間主成分分析IT-PCA、信号間多次元尺度構成法IS-MDS、時点間多次元尺度構成法IT-MDSの計算方法自体は、あらゆる信号、更に一般的にはデータ点数が同じ複数のデータ項目の解析に適用可能である。これらの計算方法自体は半導体製造装置に限定されるものではない。本実施例では、半導体製造装置601の装置監視処理を実行するのは、図6に示すネットワーク131を介して接続される装置データ監視装置133である場合の構成を説明する。実施例1においては、半導体製造装置601に備えられた計算機、記憶装置111にある半導体製造装置監視処理演算部701で装置監視処理の全ての機能を実行する例を示した。また、装置データDB132が計算機、記憶装置111内に構成される例を説明した。本実施例の装置データ監視装置133は、図7に示す半導体製造装置監視処理演算部701を備え、ネットワーク131に接続された複数の半導体製造装置601の装置監視処理を実行する。また、ネットワーク131に接続された装置データDB132は、複数の半導体製造装置601の装置データを共通に格納する。各半導体製造装置601の計算機、記憶装置111のロット、ウェハ毎処理来歴取得部702は、各装置制御装置611-615、各処理中センサ621-623の出力信号を、所定のサンプリング間隔の装置データとして取得して、ネットワーク131を介して、装置データDB132へ記録する。装置データ監視装置133は、各半導体製造装置601ごとに記憶された装置データを装置データDB132から検索して読み出し、各半導体製造装置601の装置監視処理を、実施例1,2,3に記載した処理と同様に実行する。そして、処理結果を、図6,7には記載していない出力部へ出力すると共に、または該当する半導体製造装置601の計算機、記憶装置111へ処理結果を全て送信する。装置監視処理結果のデータを受信した計算機、記憶装置111は、各表示処理部、異常検知部714により、画面、ユーザーインターフェース112へ表示、出力する。なお装置データ監視装置133は、例えば装置の近くに設置されている必要はなく、データセンターといった建物に設置されていてもよい。処理結果の出力は装置の画面、ユーザーインターフェース112でなくても、オフィスのPC端末画面、ユーザーインターフェースであってもよい。これならば複数の装置を同時に、同一の画面で管理できる。"