不適解が出てしまった時の対処。不適解とは、 ・推定値の絶対値が1以上 ・誤差分散が撹乱項の分散が0やマイナス など。これらは一般に論文には出せない。ただし「サンプル変動による不適解」と判断できれば、 論文に載せることも不可能ではない。 参考。 htt…

昨日に引き続き、共分散構造分析の勉強。Rでsemパッケージがうまくいかないので、今日はlavaanパッケージと格闘。 2．lavaanパッケージ使い方は、semパッケージとあんまり変わらない。 ただ、モデルの記述方法が違う。モデルの記述方法 =~ 潜在変数の定義 ~ …

共分散構造分析をやってみた。ツールはRとAMOS AMOS student版（AMOS5.0）が無料でダウンロードできるため使用。 使い方は以下を参照。 心理データ解析第10回(2)だけど、student版はパスが54本までしか解析できない。 あと観測変数が8変数までしか使えない。…

最初に書いた 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は今入っていないらしい最適なモ…

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

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