SEMなどの際の推定法の中にWeighted Least Square（重み付き最小2乗法）というものがある。通常の（重みをつけない）最小2乗法は標本共分散行列とモデルの共分散行列の残差が最も小さくなるようにモデルの母数を推定するが, WLSでは残差に対して異なる重みを…

ポリコリック相関係数とは, 2つの順序付きのカテゴリー変数の間に2変量正規分布を仮定して母相関を最尤推定する方法である。Olsson(1979)の方法が有名である。Rのパッケージでも計算できるがお勉強ついでに自作してみる。解説はこの論文を参考にした。 https…

はじめに 『心理統計学の基礎』のp.114に相関係数の標本分布が図示されている。確率密度関数は数学的にかなり複雑だからとの理由で省略される代わりに近似的にその期待値を求める式が載っている。 その標本誤差についてもやはり近似的に求める次のような式が…

barplot()はベクトルもしくは行列データを与えるといい感じで棒グラフや積み上げの棒グラフを作ってくれる。 その際のデフォルトの色について。 まず、ベクトルを与えると単純な棒グラフが出力される。デフォルトの色はcol="grey"である。 x <- c(4,5,7) bar…

Rでプロットするといい感じに枠と範囲を決めてくれるが、諸般の事情から枠の隅の座標が知りたくて調べたことの覚書。 結論から言えば、デフォルトではplotした範囲（xlim,ylimで指定した範囲）の4%分の長さを外にとっているっぽい。 プロットしたときに出る…

結構昔に、シングルケースのグラフをRで描く記事を書いた。 シングルケース研究のグラフをRで描く 続・シングルケース研究のグラフをRで描く なんやかんやでR歴が2年近くになった今見ると、すごい無駄な事しているなぁと思ったりする。 今回はより洗練された…

確証的因子分析を行うと、モデルがどの程度データに当てはまっているかの指標が複数算出される。その解釈の仕方については検索するとよく出てくるが、その算出については調べてもあまり出てこないので手元で色々と試してた記録を残しておく。 確証的因子分析…

因子分析シリーズ記事。 その1（とりあえず試した編） その2（因子数の決定編） その3（様々な推定法編） その4（因子軸の回転編） その5（高次因子モデルと双因子モデル編） その6（確証的因子分析編） ずっと、同じデータセットというのも飽きてきたので違…

シリーズものの記事。過去の記事はこちら。 その1（とりあえず試した編） その2（因子数の決定編） その3（様々な推定法編） その4（因子軸の回転編） その5（高次因子モデルと双因子モデル編） 前回までは基本的に探索的因子分析（Exploratory Factor Analy…

因子分析についてのシリーズものの記事。 その1（とりあえず試した編） その2（因子数の決定編） その3（様々な推定法編） その4（因子軸の回転編） 今回はもう少し複雑なモデルということで高次因子分析モデルと階層因子分析モデルのうち双因子モデルを試し…

まだまだ因子分析の続き。 その1（とりあえず試した編） その2（因子数の決定編） その3（様々な推定法編） 今回は因子の回転について。 様々な因子軸の回転法 このシリーズものの記事では、Rのpsychというパッケージのfa()という関数を使っているのだが、こ…

以前からの続きで因子分析を試してみるシリーズ。 その1 その2 今回はいろいろな推定方法について。あいもかわらず清水先生のスライドを頼りにしながら試してみる。 Rで因子分析 商用ソフトで実行できない因子分析のあれこれ データセットは前回、前々回と同…

Rは大変便利な言語なのだけれど、大規模なデータの扱いは苦手と言われている。その理由は、データをメモリに全部載せて作業するので、メモリの量によって扱える操作が制限されるとういもの。最近やや大規模なデータを扱う必要があってその処理の仕方を勉強し…

以前に書いた記事の続き。 nekomosyakushimo.hatenablog.com 前回は固有値を求めスクリープロットを書いて因子数を決定したが、因子数の決め方にも色々あるらしいので、それを試してみる。参考にしたのは関西学院大の清水先生の次のスライド。この記事の内容…

思いつきで何かをやりたくなるということがあって、今は因子分析な気分なのでやってみる。ちなみに多変量解析の経験値はほぼゼロ。 データセットの準備 ここのところマイブームがアセスメントなので、知能検査のデータセットを使うこととする。Rのパッケージ…

csvからデータを読み込むときにカンマが入った数字（29,000）みたいなものが因子型（factor）のベクトルとして読み込まれることがある。そういうときには、read.csvのオプションにstringsAsFactors=FALSEというオプションをつけると良い。こうするととりあえ…

Rで美しいグラフを描きたい欲が急に湧いたため次の本を買った。Rグラフィックスクックブック ―ggplot2によるグラフ作成のレシピ集作者: Winston Chang,石井弓美子,河内崇,瀬戸山雅人,古畠敦出版社/メーカー: オライリージャパン発売日: 2013/11/30メディア: …

実験を実施した結果、一つのフォルダに複数の.csvファイルが保存されたりする。サンプルサイズが小さかったり、処理が複雑でなければそれぞれを開いて手作業でやっても良いが、数が多くなってくると大変に時間がかかるので単純作業は機械にやらせると良い。…