「統計学関連なんでもあり」の過去ログ--- 040

No.03121　Ｒによる欠測値の推定について　　【太原】　2007/04/04(Wed) 19:26

初歩的な質問で申し訳ありません。

反復数が3の乱塊法による8水準×6水準の2要因分散分析について，欠測（反復0）が20あります。この後に交互作用効果の推定値から分散共分散行列を計算し，固有値と固有ベクトルを計算したいため，EMアルゴリズムによる欠測データの補完を考えています。

上記の計算をできればフリーソフトのRで行いたいのですが，そのための関数または反復計算を行うような関数はあるのでしょうか？Web上のRに関する説明を調べてみましたが，プログラミング初心者ということもあり，分かりませんでした。

申し訳ありませんがご助言いただけたらと思います。

No.03124　Re: Rによる欠測値の推定について　　【青木繁伸】　2007/04/04(Wed) 21:32

とても初歩的とはいえない難しい問題だと思いますが

質問の本題とは外れますが，反復0のが20ということは，20/(8*6)≒42% が欠測なんでしょうか？それを補完するというのは，大変でしょうね。

No.03129　Re: Rによる欠測値の推定について　　【太原】　2007/04/05(Thu) 08:51

ご返信ありがとうございます。

欠測42%ではどのような方法でもやはり補完は難しいのでしょうか？実験データの取り直しはできないため，6/(6*4)≒25%でもう一度考えてみたいと思います。