A Primer on Partial Least Squares Structural Equation Modeling (PLS-SEM)作者:Hair Jr., Dr. Joe,Hult, G. Tomas M.,Ringle, Dr. Christian M.,Sarstedt, MarkoSAGE Publications, IncAmazon A Primer on Partial Least Squares Structural Equation Mod…

収束妥当性

Mplusで推定にWLSMV(adjusted diagonally weighted least squares)を使った測定の不変性の確認の方法について。 最尤法は正規分布と連続変数という2つの強力な仮定をすることから、順序尺度や正規分布にならないものの測定の不変性について確認する場合には…

WLSMV を用いた測定の不変性の計算[Mplus]

Structural Equation Modeling: Applications Using Mplus (Wiley Series in Probability and Statistics Book 9) (English Edition)作者:Wang, Jichuan,Wang, XiaoqianWileyAmazon 2.9章から。 従来のCFAモデルでは、各指標・項目は理論的に仮定された1つの…

ベイジアンCFAモデル

WLSMV(adjusted diagonally weighted least squares)での測定の不変性の方法が確立しているらしい。 https://www.tandfonline.com/doi/abs/10.1080/10705511.2019.1602776 最尤法での測定の不変性についてはこちらを参照のこと。 ides.hatenablog.com サンプ…

lavaanを用いてWLSMVによる測定の不変性を計算する[R]

Kaiser-Meyer-Olkin(KMO)の標本妥当性の測度とMeasures of Sampling Adequacy(MSA)について。 中澤港さんの資料 https://minato.sip21c.org/factor-in-R-j.pdf KMO と MSA KMO とは，Kaiser-Meyer-Olkin が提唱した因子分析全体についてのサンプリング適切性…

Kaiser-Meyer-Olkinの標本妥当性の測度とMeasures of Sampling Adequacy

因子スコアを独立変数にして、重回帰分析にぶち込んだ論文がいくつか存在していることを知った。 Yakubu, Idahor and Isopa(2009) https://www.researchgate.net/publication/50194274_Using_factor_scores_in_multiple_linear_regression_model_for_predict…

因子スコアを重回帰分析の独立変数にした論文

頻度主義で行うのと特に何かが変わるというわけではないものの、Mplusでベイズ推定で因子分析を行った。 コード DATA: FILE = "HS1939.dat"; VARIABLE: NAMES = x1 x2 x3 x4 x5 x6 x7 x8 x9; MISSING=.; ANALYSIS: ESTIMATOR = BAYES; PROCESS = 2; FBITER =…

ベイズ推定による因子分析[Mplus]

データとモデルの作成 library("lavaan") Data <- HolzingerSwineford1939[,c("x1","x2","x3","x4", "x5", "x6")] model <- ' f1 =~ x1 + x2 + x3 f2 =~ x4 + x5 + x6 ' fit <- sem(model, data = Data, std.lv = TRUE) lavaanパッケージに含まれているデー…

lavaanで潜在変数の相関係数を求める[R]

相関係数をピアソンではなく、ポリコリック相関係数を使う方法。例として簡単な確証的因子分析を行う。 サンプルデータの作成 library(lavaan) data.con <- HolzingerSwineford1939[,c("x1","x2","x3")] # visualのデータのみ使用 ## 3カテゴリにしてカテゴ…

lavaanで順序尺度のSEMを行う[R]

信頼性係数はα係数が使用されることが多いが、問題点が指摘されている。α係数は因子構造を無視して一次元性の検証をしているという点である。別の言い方をすると、各因子での真の得点が共通している(一元性)と仮定しているが、ほとんどのケースでは一元性は…

同族テストモデル・タウ等価モデル・平行テストモデル[R]

IPIP-NEOデータの呼び出しと格納 library("psych") data(bfi) # IPIP-NEOデータ d1 <- bfi[1:10] # 因子分析に使用するのは１～10列目。2つの因子のみ。 d1 <-na.omit(d1) # 欠損値のあるケースを削除 モデルと実行 library(lavaan) model <- ' Ag =~ A1 + A…

lavaanで順序カテゴリカル因子分析[R]

ω(オメガ)係数は内的整合性を確認するための指標である。近年はCronbachのαより良い方法と言われることが多い。Cronbachのαは因子構造を無視して、一次元性の検証していることが原因である。 データ psychパッケージに含まれるIPIP-NEOのデータbfiを使用する…

ω係数[R]

階層因子分析（hierarchical factor analysis）について松田淑美・狩野裕「高次因子分析モデルと階層因子分析モデルについて」から説明を引用しよう。 http://www.sigmath.es.osaka-u.ac.jp/~kano/research/meeting/20050827_BSJ33/matsudaABS.pdf 階層因子…

lavaanパッケージを用いて階層因子分析を推定する[R]

二次因子分析モデル(second order factor analysis)とは「斜交解による通常の因子分析の結果，因子聞に強い相闘が観察されたときには，それらの因子に共通して影響を与える上位の因子を想定することができます。このように，通常の因子分析における複数の因…

lavaanパッケージを用いて二次因子分析モデルを推定する[R]

複数の観測変数によって1つの構成概念が規定され、その構成概念が別の観測変数群に影響を与えているようなモデルをMIMIC(Multiple Indicator Multiple Cause）モデルと呼ぶ。 MIMICモデルのデモを見ていると、Mplusの例5.8を使ったものが多かったので、それ…

確証的因子分析におけるMIMICモデル[R][Mplus]

時々気になるサンプルサイズの話である。 pubmed.ncbi.nlm.nih.gov サンプルサイズのについての記述 Henrica et al.(2005)では、下記のように書かれてある。 サンプルサイズ 探索的因子分析も確認的因子分析も、信頼できる結果を得るためには相応のデータ量…

因子分析のサンプルサイズは項目数×7で求められるのか

前回のエントリ、２値のカテゴリカル因子分析と潜在クラス分析の結果の差異で因子分析と潜在クラス分析のグループ分けの違いを分析した。 今回は以前のエントリで使用した児童向けウェクスラー式知能検査（Wechsler Intelligence Scale for Children; WISC)…

２値カテゴリカル因子分析と潜在クラス分析のシミュレーション, その２

２値のカテゴリカル因子分析と２値の潜在クラス分析の結果どのように違うのかをシミュレーションしてみたい。因子分析も潜在クラス分析もグループ分けをする手技である。しかし、両者はグループの分け方が異なる。因子分析が変数の近さでグループを作るのに…

２値のカテゴリカル因子分析と潜在クラス分析の結果の差異

ノーマルな因子分析は連続変数を用いる。その拡張としてカテゴリカル変数でも因子分析は可能である。今回は、連続変数で推定した結果とカテゴリカル因子分析の推定はどのくらい異なるのかをシミュレーションしてみたい。 大前提としてカテゴリカル変数といっ…

因子分析とカテゴリカル因子分析の結果はとのくらい異なるか

追記： 019/09/18 以前、psychパッケージのfa.poly functionでカテゴリカル因子分析を書いていたが、非推奨の方法だったようだ。psychパッケージの仕様書には次のように書いてある。 fa.poly Deprecated Exploratory Factor analysis functions. Please use …

カテゴリカル因子分析

落穂拾い的なエントリ。 以前のエントリー(http://ides.hatenablog.com/entry/2019/04/08/171145)で取り上げた清水和秋「因子分析的研究におけるmisuseとartifact」の続き部分である。 kansai-u.repo.nii.ac.jp 平行分析もMAPも主成分を使用 Horn（1965）の…

因子分析に関するメモ

清水和秋(2018)「因子分析的研究におけるmisuseとartifact」『関西大学社会学部紀要』 49(2): 191-211. kansai-u.repo.nii.ac.jp 標本サイズ Cattell（1978）は、標本サイズと分析対象の変数の数との比として、3 対 1 を目安としている。この比についての基…

因子分析と標本サイズ

Mplusでの平行分析の実行の仕方は以前に書いたが、Rでも実行できる。 今回もbfiデータの1~25列目を使用する(参考: http://ides.hatenablog.com/entry/2019/03/19/093726)。 library("psych") data(bfi) d1<-bfi[1:25] 平行分析は次のように指定する。 fa.par…

Rで平行分析

因子分析をしていて混同しがちなのは因子負荷(因子パターン)と相関係数ではないかと思う。 誰も教えてくれなかった因子分析: 数式が絶対に出てこない因子分析入門作者: 松尾太加志,中村知靖出版社/メーカー: 北大路書房発売日: 2002/05/01メディア: 単行本（…

因子負荷と相関係数の違い

探索的因子分析において因子数を決定する基準として使われるMAPについて。 MAP:Minimum Average Partial correlationであり、日本語だと最小平均偏相関になる。 Velicerによつて開発された方法である。 link.springer.com 因子数を決める基準はどれが適切か…

MAP(最小平均偏相関)

Mplusで因子分析をしたことがなかったので、勉強がてらのメモ。 今回行うのは、下記のエントリで行った探索的因子分析の内容と同じである。 ides.hatenablog.com データはこちらから Mplusのコード TITLE: Exploratory Factor Analysis by Mplus, using NEO …

Mplusで因子分析 & 因子選択の基準

因子分析の参考書のブックガイドを書いておきたいと思う。 誰も教えてくれなかった因子分析: 数式が絶対に出てこない因子分析入門作者: 松尾太加志,中村知靖出版社/メーカー: 北大路書房発売日: 2002/05/01メディア: 単行本（ソフトカバー）購入: 1人 クリッ…

因子分析のブックガイド

因子分析の勉強メモ。 カイザー基準は、カイザー・ガットマン基準と呼ぶこともある。これはKaiser(1960）よりも先にGutttman(1954）が先に指摘があったため、連名での表現である。 固有値が１以上のものを因子とする方法である。それは、１つの因子にしか負…

カイザー基準についてのメモ

因子分析について少し勉強をしているところである。因子分析は心理学の人が使うことが多く、心理学徒ではない僕は因子分析をほとんど使ったことがなく、勉強もろくにしたことがないので自分の勉強を兼ねてのエントリーである。 今日、取り上げるのは下記の文…

柳井晴夫「因子分析法の利用をめぐる問題点を中心にして」

この原稿は統計学勉強会のためにかかれたもので、学校で習うあたりまでの簡単な解説である。 因子分析とは (探索的)因子分析とは、観測できる変数から観察できない因子（潜在変数）を見出す探索し発見する方法である。このように説明してもよくわからないの…

Rで因子分析　基礎