diagonizable.統計解析に行列を使うとき、データセットが単純に理解できるような基底の取り直しをすることがよくある。基底の取り直しは線形変換で行うので、線形変換で互いに行きあえる行列同市は「相似(似ている)」関係にあると言う。相互に相似な関係にある行列の中で、最も単純な行列を代表とすることがあり、それを使うことでデータ解釈が最もわかりやすくなる。相似な行列の中で最も単純なのは、ジョルダン標準形のそれである。対角化可能な行列なら、ジョルダン標準形は対角化行列である。実際、統計解析で現れる行列の多くは、対角化可能である。任意の正方行列$M$はジョルダン標準形$J$と逆行列を持つ行列$R$とを用いて、次のように表せる。R'JR'^{-1}$となるとき、$M$と$M'$は相似であると言う。以下では、ジョルダン標準形が対角行列であるようなもの(対角化可能行列)に絞って話を進めるが、特に断らなければ、提示している性質等は対角化可能行列に限るものではない。相互に相似な行列は見かけは違うが、「本質的に同じ」なので、その「本質」に関する性質は共通である。1-1逆行列を持つ行列Rを使って、対角行列Jを変換した行列M($M=RJR^{-1}$)がJと相似であることを、固有値が一致していることを確認して確かめよ。ランク、行列式、トレースが同一であることを確かめよ。1-2行列が相互に相似かどうかは、(すべての固有値が異なるならば)固有値のセットが等しいことで確認できる。したがって相似な行列に共通する特徴は、固有値と紐づけて算出できるはずである。1-3$d$個の異なる固有値を持つ行列を作り、それと相似な行列を複数作れ。それらの固有値のセットが一致すること、トレース・行列式が一致することを確かめよ。トレースの計算は、対角成分の和を取る方法と、固有値の和を取る方法の両方を実施せよ。行列式の計算もdet()関数を使用する方法と固有値の積を取る方法の両方を実施せよ。対称行列$M=M^T$である正方行列を対称行列と言う。I$$と対角化で現れる行列$R$が特別な性質(直交行列)を持つ。2-1対角行列と直交行列によって相似な関係にある行列を複数作り、それらが対称行列であることを確かめよ。直交行列をランダムに発生するには、GPRrotation::Rnadom.Start()関数を用いるとよい。行列には相似という間柄があり、それは固有値で決まることが分かった。相似な行列の間で保存される性質があることもわかった。R^{-1}$という変換が相似な行列を作ることもわかった。R^T$と書け、Rが直交行列であることもわかった。これらは、一般的な行列の話であったが、この節では、統計解析でよく登場する特別な行列を取り上げ、その場合に固有値・相似行列の特徴がどのような性質を持つかを確認する。距離行列$n$個の要素があって、それらの要素ペアの距離を行列とすることがある。要素間の関連を調べるときの基本的なデータ構造であり距離の場合、自身から自身への距離は0である。自身以外への距離は正である。AからBへの距離とBからAへの距離は同じなので、距離行列は、次の性質を持つ。非対角成分は正このような行列では、以下のような特徴がある。k$の分散共分散行列が算出できる。対角成分は$k$個の変数の分散である。トレースは$k$個の変数の分散の和である。$k$個の変数の分散の和が固有値の和である。これは距離行列の場合と同じ。t(eigen.out[[2]]),8)eigen.out[[1]]```分散共分散行列は、対象行列であって、対角成分はすべて正。以下に示すように、固有値もすべて正。固有値がすべて正であるような対称行列は正定値行列を呼ばれている。正定値行列の特殊例でもある、単位行列の場合は、以下のように、「ベクトルの長さの二乗」である。これが0ベクトルなら0となり、そうでなければ正になる、というのは、ベクトルの「長さ」を定義するのに適切な性質を持つ行列を正定値行列と言う、というようにも理解できる。逆に言えば、統計解析で登場する行列は、この性質を満たすものが多いことともなる。v$の値を決めるという意味では、「計量」を定めるとも言える。対角成分がすべて正である対角行列と回転行列から作る行列は正定値行列であることになる。vのみを取り出し、そのrangeを確認。統計解析でよく使われる。$$N(\mu,\Sigma)$$とパラメタ表現される。変数の数$n$個の多変量正規分布の場合、$\mu$は長さ$n$のベクトルであり、各変数の平均値を表す。n$行列であり、正定値行列であることになっている。RJR^{-1}$と対角化したときに、$R$が直交行列($R^{-1}=R^T$)であるようなものであった。これを利用して、多変量正規分布からの乱数発生をしてみる。2$の単位行列であるようなもののことである。3-12変量正規分布において、各変数の分散を(1,1)からずらして、乱点を発生せよ。分散共分散行列は対角行列のままとする。3-24-1で用いたSigmaと相似な行列で対称なものは正定値行列であるから、mvtnorm::rmvnorm()のsigma引数として適当なはずである。そのようなsigmaを作り、乱点を発生せよ。変換に用いた直交行列$R$を利用して、乱点を変換すると、楕円の軸を正すことができる。以下のコードを読み、変換した分散共分散行列に従う2変量正規分布の軸が補正されていることを、コードの各行の処理にコメントを加えることで確認せよ。4-1xの各変数の標本平均を求めよ。xの分散共分散行列を求めよ。求めた平均値ベクトルと分散共分散行列とに基づく、多変量正規乱数を発生し、ペアワイズにプロットせよ。4-3回転行列を用いて、xを回転し(軸を取り直し)、その標本平均と分散共分散行列とを求め、4-1のそれと比較せよ。各変数の平均も分散も異なることを確かめよ。4-4分散は異なるが、分散の和は変わっていないことを、分散共分散行列のトレースを計算することで確かめよ。同様に、分散共分散行列の固有値の和としてトレースを計算することでも確かめよ。4-5データセットxの軸の取り方を回転することによって変えても、各軸の分散の和は変わらない。分散の和は変わらないが、分散の総和の分配具合は変わる。そのことを、多数の回転行列をランダムに発生させ、データセットを回転した上で各軸方向の分散を計算し、その分散の値を大小順序でソートしてプロットし、分配具合が変わることを具体的に確認せよ。4-6xの分散共分散行列の固有値を求めることで、分散共分散行列を対角化するような回転をしたときの、各軸の分散を求めることができる。その分散を4-4のプロットに重ねてプロットせよ。4-74-4,4-5と同様のプロットを変数の数を増やして実施せよ。GWAS/GWASメタアナリシスとその周辺解析用のアプリケーションへのリンク複合遺伝性疾患とは----------News----------20万アクセス2008/07/2040万アクセス2010/09/12-------------------------はてなブログをはじめよう!ryamada22さんは、はてなブログを使っています。ブログを報告するスターをつけました引用をストックしましたストック一覧を見る閉じる引用するにはまずログインしてくださいログイン閉じる引用をストックできませんでした。再度お試しください閉じる限定公開記事のため引用できません。