最初に書いた mclustBICとMclustの出力は 本来一緒のはずらしいけど， 今回は異なってしまった理由は不明...なぜなんだろう(´Д` ) あとクラスターの分布が違うのは 使った混合分布の型(球、楕円球)、体積、形と軸の向きが 異なるから． （最初"EEE"，二回目"…

いろいろ調べてたら もう一つやり方見つけたただこれやったら全然結果変わってしまった(´Д` )とりあえずメモ。 > res > plot(res) これよりEEVのクラスター数7が最適 > res > cl > cl[1:10] 1 2 3 4 5 6 7 8 9 10 6 7 6 7 6 6 7 3 6 6 > table(cl) cl 1 2 3 …

モデルに基づいたクラスター分析http://mjin.doshisha.ac.jp/R/29/29.html を参考にやってみたmclustパッケージをダウンロードまず書かれているように plot(EMclust) ってやったけどエラーどうやらmclustパッケージに EMclustは今入っていないらしい最適なモ…

apply関数を使って 行列処理をしてた時の エラーというか困ったこと func list2 for(i in 1:3){ if(list[i]==list[7]) list2[2,i] if(list[i+3]==list[8]) list2[3,i] } if(sum(list2[2,])==1) return(list2[1,list2[2,]==1]) else return(list2[1,list2[3,]…

原因が本当にこれかは不明だけど「はい」「いいえ」みたいな二値データは 基本的には因子分析では扱われなくて、因子分析したいなら カテゴリカル因子分析ってやつを使うらしいhttp://home.hiroshima-u.ac.jp/keiroh/maeda/statsarekore/gedoufa.htmlテトラ…

アンケート結果を因子分析して 因子を抽出しようとしているのだけど、なかなかうまくいかない 使っているのは、 factanal関数 factanal2関数（青木先生） fa関数（psychパッケージ） とかいろいろデータ形式は 「いいものは長く使いたいですか？」：「はい」…