S-Plus と R での固有ベクトル計算

　例の CWRES の計算でちょっとはまってしまっているのです．S-Plus と R とで計算結果が異なるようだ，ということでいろいろ調べてみると，どうやら，固有ベクトル計算のところで違いが出てしまっているらしい．はっきりしたらまた別のところに書いておこうと思いますが，一言で言うと，S-Plus では「普通に」分散共分散行列を作っただけでは，それは厳密には対称行列になってくれなかった，ということ．そのために，うまく直行行列を構成できなかったのです．
　今検証中．