2006-02-13から1日間の記事一覧

偏相関係数の計算 メモ

R

Rの該当ページ 日本語のページ Rのページに従って計算。 8x8 相関係数行列の要素(要素数 64)をカンマ区切りで代入して、上記Rの該当ページに沿って実行(Matrix化、逆行列化、対角化、行列演算) > corMatlist<-c(1,0.014612303,0.026175719,0.002868667,0.023…

四捨五入とNaNの扱い メモ

相関係数・偏相関係数の計算ソースを書くと、それ以外の数値の扱いとも同様だが、分母がゼロになって割り算が不能だったり、切りのよい結果(1や0など)が返るべきところで、0.9999999999999999 などが出て困る。 NaN サイト "Not a Number"の略 NaNか否か…