コメント

1
尺度化について
狩野 裕
大阪大学人間科学部
岡さんのまとめから
■尺度作成における項目選択 (本研究の場合⋯)
・従来の尺度の項目数: 30項目 (各因子15項目ずつ)
・SEFAによる項目選択: 14項目 (各因子 7項目ずつ)
1つの因子に含まれる項目数は,どの程度が妥当であるのか?
■信頼性に関して (α係数の算出による内的整合性の検討)
・従来の尺度のα係数: .89〜.90 (15項目)
・新しい尺度のα係数: .77〜.80 ( 7項目)
α係数がどの程度の値であると,最適だと言えるのか?
2
3
安部さんのまとめから
• 十分な信頼性を確保することは、尺度作成
において非常に重要である
• 被験者の回答への負担を考慮することは
研究者としての責務である!
• 項目数を少なく、しかも十分な信頼性と
妥当性を有した尺度を作成するためには?
新たな基準の必要性
4
いくつかの概念
• 信頼性
– 項目の中に含まれる真の値の割合
– 真の値のばらつき/項目のばらつき
– V(X)=V(x)+V(e)
• α係数
– 信頼性の簡便推定値
• 尺度の(内容的)妥当性
– 測定すべき側面を全て聞いているか
– バラエティのある項目になっているか
• 適合度
– 因子分析モデルがデータと適合しているかどうか
5
尺度の使い方
• 構成概念の「ものさし」を作成し,合計得点
で被験者を「ものさし」上に布置する
– 信頼性と内容的妥当性
• 共分散構造分析(SEM)の指標として使う
6
補足:SEMの指標
• 指標の数
– 4以上が薦められる
– 2,3でも可能
• 信頼性
– 高い方がよい
– 低くても「誤差を懐柔して」分析可能
7
1つの因子に含まれる項目数-1
• 理論的には項目数は問題ではない
– 信頼性・内容的妥当性を達するよう項目選択
– 項目数はその結果
– 実行上は少ない方がよい
• 誤差(回答の揺れ)を取り除くことが主目的
の場合は項目数=2で十分
(ただしSEMを使うことが前提)
8
1つの因子に含まれる項目数-2
• 信頼性と内容的妥当性の折衷
– 項目が少ない=>両妥当性とも低い
– 信頼性を重視して項目を増やす
=>内容的妥当性が低くなる
– 内容的妥当性を重視して項目を増やす
=>信頼性が低くなる
• 内容的妥当性を保ったままで信頼性が
最も高くなるように項目選択する
– 内容的妥当性の数理的測度はない
9
冗長性の検討 r>0.8?
• データの再現性はどのくらい?
• 反復測定(一週間おいて2回測定)
2回の相関
男子 女子
身長
0.997 0.992
タレントの好感度
0.84 0.84
減税†
0.85 0.62
別姓††
0.85 0.73
項目
† 財政改革を止め減税すべきである 賛成 どちらかといえば賛成 どちらでもない どちらかといえば反対 反対
†† 男女別姓を認めるべきである
賛成 どちらかといえば賛成 どちらでもない どちらかといえば反対 反対
10
冗長性の検討 r>0.8?
• 態度・意見・印象などを問う質問紙では,
同一項目間の相関係数は0.8程度である
ことが多い
• 0.8以上の項目があれば1つの項目を外す
というのはひとつの見識
– 異なった項目間の相関が0.8以上であったとし
ても,同一項目とみなしてもよいかの吟味は
必要
11
α係数の値-1
• 信頼性は(他の条件が同じならば)
高い方がよい
– 信頼性が高い(すぎる?)尺度では
• 内容的妥当性がとぼしい
• 同じことを聞いている項目を削ることが可能
ということが多い
– このことから「適度な大きさのα係数」という
基準がでてくるのだが...
12
0.70
α係数の値-2
• 信頼性が低いと
尺度間の相関
や回帰係数など
に偏りが生じる
• どの程度の偏り
を認めるか?
構成概念
A
構成概念
B
0.63
尺度A
α=0.9
尺度B
α=0.9
0.56
尺度A
α=0.8
尺度A
α=0.5
0.35
尺度B
α=0.8
尺度B
α=0.5
13
α係数の値-3
• 信頼性とα係数は異なるが
– 一般には良い近似である
– 上記が正しくあるためには,因子モデルの
適合度が高くないといけない
• 適合度の吟味はα係数を使うための
必要条件
• 適合度が悪いときは
– α係数の値が信頼できない
– 尺度の1次元性が確立できていない
14
「恩恵の知覚」尺度:n=653
除いた変数
カイ2乗値
df
GFI
AGFI
CFI
RMSEA
α係数
真の信頼性係数
悪い適合度による
15
-
460.11
90
0.914
0.874
0.892
0.080
0.90
0.86
分析変数の数
12
11
2 6 14
2 9 13 14
172.96
107.37
54
44
0.955
0.969
0.935
0.954
0.945
0.967
0.058
0.047
0.86
0.86
0.85
0.85
7
たくさん
14.27
14
0.994
0.987
1.000
0.006
0.80
0.80
15
真の信頼性係数について
• 適合度が低い
– 1因子のモデルで説明しきれていない
• 独自(誤差)因子間に相関を設定
– LM検定を実行し,Cov(ei,ej) がゼロでない組を
検出する
16
誤差相関
適合度:χ2=250.375; GFI=0.950 AGFI=0.928;
CFI=0.952; RMSEA=0.056
参考:変数選択で落とした変数:2,4,13,14
17
Summary and recommendation
• 15変数のモデルは1因子モデルの適合度が
悪いため,α係数の値(0.90)は信頼できない
– 実際,適合度の良い「誤差相関」1因子モデルによる
と信頼性は0.86
• 誤差相関があること(適合度が悪い)は尺度の
1次元性に反するのでいくつかの項目を落とす
• 不適切な3 or 4 項目を落とすことで,適合度の
よいモデルが構成できる
– 1次元性の確保
– 信頼性=0.86...15変数での信頼性を確保
18
まとめ
• α係数の大きさは,その後の解析で生じる
バイアスをどの程度許容するかに対応する
• α係数が信頼できるためには1因子モデル
の適合度が高いことが必要
• SEFAは適合度によって変数を選択する
プログラム
http://koko15.hus.osaka-u.ac.jp/~harada/factor/
• 内容的妥当性をにらみながら,適合度と信頼
性が高くなるよう変数選択する
19
補足:適合度の評価方法
• 中小標本のとき
– カイ2乗値(有意確率)をみて,モデルが棄却
されないことを確認
• 大標本のとき
– 適合度指標をみる
– GFI,AGFI,CFI...「0.90以上」は最低条件.
0.95前後以上が望ましい
– RMSEA...0.05以下が望ましい
20
補足:尺度とSEMの指標
信頼性
妥当性
適合度
X1+X2+X3+X4
合計得点 SEMの指標
高
そこそこ
高
高
そこそこ
高
21
多母集団の同時分析
狩野 裕
大阪大学人間科学部
22
検証的因子分析のメリット
•
•
•
•
•
多母集団の同時分析
因子平均の統計的比較
因子に関する種々の仮説の統計的検討
因子回転が不必要
探索的因子分析不可能でも,検証的因子
分析が可能な場合がある
23
多母集団の同時分析のメリット
• 母集団間の違いを統計的に検討できる
– 因子構造の違い
– 因子平均の違い
• ある変数に欠測が集中したデータの分析
• 個別分析できないほど小さなデータの分析
• その他
24
「山田・豊本」分析のながれ
• 二つの母集団の因子構造を「統計的」に
比較したい
– 探索的因子分析結果の比較は記述的
• 検証的因子分析
– 個別分析
• 母集団ごとに分析
– 同時分析
• 因子負荷の相等性の検定
因子不変性
• 配置不変
– 全母集団に,同じパスを引く
モデルが適合
• 測定不変(弱因子不変)
– 全母集団に,因子負荷の値が
同一のモデルが適合
X1
X2
X3
X4
X5
X6
因子
相関
母集団A
F1 F2
*
*
*
*
*
*
1
*
*
1
母集団B
F1 F2
*
*
*
*
*
*
1
*
*
1
X1
X2
X3
X4
X5
X6
因子
相関
母集団A
F1 F2
0.5
0.6
0.7
0.6
0.7
0.8
1
*
*
1
母集団B
F1 F2
0.5
0.6
0.7
0.6
0.7
0.8
1
*
*
1
X1
X2
X3
X4
X5
X6
因子
相関
母集団A
F1 F2
0.5
0.6
0.7
0.6
0.7
0.8
1 0.4
0.4
1
母集団B
F1 F2
0.5
0.6
0.7
0.6
0.7
0.8
1 0.4
0.4
1
• 強因子不変
– 測定不変かつ,因子相関行列が
同一のモデルが適合
25
モデル修正
• 初期モデルでは適合
が十分でない
• V8、F3のパスを追加
する
どのパスを加えると
最も適合がよくなる
のか
震
初期モデル
V8,F3のパスを追加
度 不安 うつ 混乱 不安 うつ 混乱
7
F1
F2
F3
F1
F2
F3
V1
.96
.95
V8
.91
.46
.64
V9
.98
.97
V10 1.22
1.23
V11
.81
.82
V18 -.41
-.43
V19 1.17
1.17
V20 1.05
1.05
V21 1.03
1.04
V3
.71
.72
V5
.53
.53
V7
.80
.80
V12
.63
.63
V14
.52
.52
V15
.43
.44
V2
.76
.76
V4
.36
.36
V6
.78
.78
V13
.89
.90
V16
.77
.77
V17
.72
.71
736.1
χ^2値
792.1
0.926
GFI
0.919
0.892
CFI
0.881
0.056
RMSEA
0.059
26
27
スタートのモデル
震度7のデータ
の個別分析(1)
適合度
χ^2値
792.1
df
186
P値
<0.001
GFI
0.919
CFI
0.881
RMSEA
0.059
LM検定結果
パス
χ^2値
V8,F3 55.776
V19,F3 47.520
V16,F1 46.504
V22,F2 33.012
V1,F2 32.650
不安
うつ
混乱
V V V V V V V V V
1 8 9 10 11 19 20 21 22
V V V V V V
3 5 7 12 14 15
V V V V V V
2 4 6 13 16 17
28

初期モデルの
 2値
ー
792.1
ー
2
値について
2
V8,F3の  値
(LM検定による
予測値)
736.324
55.776
≒
V8,F3を加えた
≒
736.1
モデルの  2 値
LM 検定はいつ終了させるの
か?
• 内容的吟味
• 統計的基準
– χ12(0.05)=3.841
– χ12(0.005)=7.882…..default of LISREL
– χ12(0.05/検定の総数)….検定の多重性を考慮
• χ12( 0.05/42)=**
[註:42=21×3-21]
– AIC でみるときは 2
– CAICでみるときは 1+log n
• 1+log956=7.86
– スクリープロットの要領で
• 飛びぬけて大きな統計量がなくなるまで続ける
29
30
個別分析の
まとめ
• 配置不変は成立し
ない
• 「配置不変」とは
– パス(因子負荷の
位置が一致すること
震度4
不安 うつ 混乱
.87
.50
.46
.84
1.26
-.39
.89
-.72
.52
.97
.79
1.80 -.97
.71
.56
.83
.49
.67
.25
.74
-1.31 1.82
-1.08 1.84
.91
.93
.73
V1
V8
V9
V10
V11
V18
V19
V20
V21
V3
V5
V7
V12
V14
V15
V2
V4
V6
V13
V16
V17
震度7
不安 うつ 混乱
.68 .37
.53
.58
.97
1.22
.83
-.78
.54
1.17
1.04
1.39 -.44
.72
.54
.79
.62
.51
.42
.79
-.84 1.20
.81
.90
.45
.41
.70
31
同時分析
それぞれの因子負荷
が母集団間で等しい
(測定不変)モデルか
らスタート
V1
V8
V9
V10
V11
V18
V19
V20
V21
V3
V5
V7
V12
V14
V15
V2
V4
V6
V13
V16
V17
震度4
不安 うつ 混乱
.79 .16
.53
.50
.91
1.24
-.20
.85
-.75
.53
1.07
.91
1.52 -.64
.72
.55
.82
.55
.59
.33
.76
-1.10 1.53
-.64 1.41
.89
.35
.58
.71
震度7
不安 うつ 混乱
.79 .16
.53
.50
.91
1.24
-.20
.85
-.75
.53
1.07
.91
1.52 -.64
.72
.55
.82
.55
.59
.33
.76
-1.10 1.53
-.64 1.41
.89
.35
.58
.71
32
最終モデル(標準
解)
しんぱいで いらいらして おちつかない
ねむれなかったり とちゅうで めがさめる
じしんのゆめや こわいゆめを みる
ふいに じしんを おもいだす
テレビや しんぶんを みると こわい
こわくない じぶんは へいきだ(逆転項目)
ちいさな おとに びっくりする
いつも ゆれているような きがする
また じしんがくるのが しんぱいだ
わけもなく かなしくて なにもしたくない
あじや いたさを かんじなくなった
ひとりぼっちに なったような きがする
じぶんの せいだと おもう
だれとも はなしたくない
じしんあそびを する
Q01
Q08
Q09
Q10
Q11
Q19
Q20
Q21
Q22
Q03
Q05
Q07
Q12
Q14
Q15
むしゃくしゃして らんぼうになり すぐかっとするように なった Q02
テレビゲームや どくしょなどに ねっちゅうする
Q04
べんきょうが つまらなくなった
Q06
すぐわすれたり おもいだせない
Q13
あたまいたや おなかいたなど からだのぐあいが わるい Q16
かんたんなことが できなくなった
Q17
m
χ^2値
CON6wald 1113.9
不安
.58
.32
.62
.74
.62
-.48
.64
.54
.92
震度4
うつ 混乱
.30
-.20
.34
-.39
.63
.58
.65
.53
.72
.22
-.65
-.38
.24
.54
.91
.88
.58
.44
.56
df
P値
GFI
381 <0.001 0.944
震度7
不安 うつ 混乱
.58
.15
.32
.30
.62
.74
.62
-.48
.34
.64
.69
.92 -.39
.63
.58
.65
.53
.54
.37
.54
-.65
.91
-.38
.88
.58
.24
.30
.56
CFI RMSEA
0.927 0.032
33
Wald検定
震度4
推定値
不安 うつ 混乱 不安
.79
14.73
.53
.50 10.14
.90
27.70
1.23
-.33 23.96
.85
27.86
-.74
.52 -13.33
1.07
29.13
.81
19.21
1.52 -.64
19.97
.72
.55
.82
.55
.68
.25
.76
-1.08 1.51
-.61 1.39
.89
.36
.67
7.59
.71
震度7
t-値
うつ
混乱
9.35
-5.30
9.12
-8.28
27.75
25.06
28.54
22.54
22.48
6.72
-8.41
-5.61
22.77
11.86
12.86
24.71
11.18
24.22
V1
V8
V9
V10
V11
V18
V19
V20
V21
V3
V5
V7
V12
V14
V15
V2
V4
V6
V13
V16
V17
推定値
不安 うつ 混乱 不安
.79 .25
14.73
.53
.50 10.14
.90
27.70
1.23
23.96
.85
27.86
-.74
.52 -13.33
1.07
29.13
1.02
21.14
1.52 -.64
19.97
.72
.55
.82
.55
.51
.42
.76
-1.08 1.51
-.61 1.39
.89
.36
.46
7.59
.71
t-値
うつ
4.10
混乱
9.35
9.12
-8.28
27.75
25.06
28.54
22.54
17.89
9.99
-8.41
-5.61
22.77
11.86
12.86
24.71
7.67
24.22
34
(統計的)解釈
• 構造や値が異なる箇所が明示的にわかる
ので解釈しやすくなる
• 配置不変は成立しない.しかし
– 異なる箇所の推定値は小さい
– 異なる箇所のt-値は大きくない
• 非常に近い構造と考えられる