[R]{psych} scoreItems() : 項目合成尺度のスコアの計算などを行う

説明

項目合成尺度のスコアおよびクロンバックα、グットマンλ6、項目-全体相関係数(item-whole-correlations)を計算する。

n個の項目に対するk個のキー(-1, 0, 1)を含む行列あるいはデータフレームと、N人✕m個の項目からなる行列あるいはデータフレームをもとに、各人および各スケールの合計および平均スコアを計算する。

上記に加え、クロンバックのα係数、グットマンのλ 6、平均r、尺度の相互相関、アイテム✕尺度の相関(生相関および項目のオーバーラップを補正した相関の両方)を報告する。

指定すれば欠損値を中央値あるいは平均値で置換(impute)することができる。逆転項目の調整も可能である。キー行列の生成のためにはmake.keys()という便利な関数も用意されている。インプットが正方行列の場合、インプットは分散共分散行列とみなされ、スコアの計算はなされない。しかし項目統計は報告される。(これはcluster.cor()と同様の機能である)。response.frequenciesは報告する、多値(polytomous)あるいは多選択肢である各応答カテゴリごとの項目裏書(endorsement)の頻度を報告する。

使用法

scoreItems(keys, items, totals = FALSE, ilabels = NULL,missing=TRUE, impute="median",
delete=TRUE, min = NULL, max = NULL, digits = 2)
score.items(keys, items, totals = FALSE, ilabels = NULL,missing=TRUE, impute="median",
delete=TRUE, min = NULL, max = NULL, digits = 2) 
response.frequencies(items,max=10,uniqueitems=NULL)

keys にキー行列というオブジェクトを、itemsに生データを渡すのが通常の使用方法である。(本関数の使用に先立ってはキー行列という特殊なオブジェクトの生成が必須である!)

実行例

data(bfi)
keys.list <- list(agree=c("-A1","A2","A3","A4","A5"),
                  conscientious=c("C1","C2","C3","-C4","-C5"),extraversion=c("-E1","-E2","E3","E4","E5"),
                  neuroticism=c("N1","N2","N3","N4","N5"), openness = c("O1","-O2","O3","O4","-O5")) 
keys <- make.keys(bfi,keys.list)
scores <- scoreItems(keys,bfi,min=1,max=6)
summary(scores)

上記実行すると以下の出力が得られる。

Call: scoreItems(keys = keys, items = bfi, min = 1, max = 6)

Scale intercorrelations corrected for attenuation 
 raw correlations below the diagonal, (unstandardized) alpha on the diagonal 
 corrected correlations above the diagonal:
              agree conscientious extraversion neuroticism openness
agree          0.70          0.36         0.63      -0.245     0.23
conscientious  0.26          0.72         0.35      -0.305     0.30
extraversion   0.46          0.26         0.76      -0.284     0.32
neuroticism   -0.18         -0.23        -0.22       0.812    -0.12
openness       0.15          0.19         0.22      -0.086     0.60

更に詳細な出力を得るにはprint()を使用する。

> print(scores)
Call: scoreItems(keys = keys, items = bfi, min = 1, max = 6)

(Unstandardized) Alpha:
      agree conscientious extraversion neuroticism openness
alpha   0.7          0.72         0.76        0.81      0.6

Standard errors of unstandardized Alpha:
      agree conscientious extraversion neuroticism openness
ASE   0.014         0.014        0.013       0.011    0.017

Average item correlation:
          agree conscientious extraversion neuroticism openness
average.r  0.32          0.34         0.39        0.46     0.23

 Guttman 6* reliability: 
         agree conscientious extraversion neuroticism openness
Lambda.6   0.7          0.72         0.76        0.81      0.6

Signal/Noise based upon av.r : 
             agree conscientious extraversion neuroticism openness
Signal/Noise   2.3           2.6          3.2         4.3      1.5

Scale intercorrelations corrected for attenuation 
 raw correlations below the diagonal, alpha on the diagonal 
 corrected correlations above the diagonal:
              agree conscientious extraversion neuroticism openness
agree          0.70          0.36         0.63      -0.245     0.23
conscientious  0.26          0.72         0.35      -0.305     0.30
extraversion   0.46          0.26         0.76      -0.284     0.32
neuroticism   -0.18         -0.23        -0.22       0.812    -0.12
openness       0.15          0.19         0.22      -0.086     0.60

 In order to see the item by scale loadings and frequency counts of the data
 print with the short option = FALSE

print()にshort=FALSEオプションを指定すると最も詳細な出力が得られる。

> print(scores,short=FALSE)
Call: scoreItems(keys = keys, items = bfi, min = 1, max = 6)

(Unstandardized) Alpha:
      agree conscientious extraversion neuroticism openness
alpha   0.7          0.72         0.76        0.81      0.6

Standard errors of unstandardized Alpha:
      agree conscientious extraversion neuroticism openness
ASE   0.014         0.014        0.013       0.011    0.017

Average item correlation:
          agree conscientious extraversion neuroticism openness
average.r  0.32          0.34         0.39        0.46     0.23

 Guttman 6* reliability: 
         agree conscientious extraversion neuroticism openness
Lambda.6   0.7          0.72         0.76        0.81      0.6

Signal/Noise based upon av.r : 
             agree conscientious extraversion neuroticism openness
Signal/Noise   2.3           2.6          3.2         4.3      1.5

Scale intercorrelations corrected for attenuation 
 raw correlations below the diagonal, alpha on the diagonal 
 corrected correlations above the diagonal:
              agree conscientious extraversion neuroticism openness
agree          0.70          0.36         0.63      -0.245     0.23
conscientious  0.26          0.72         0.35      -0.305     0.30
extraversion   0.46          0.26         0.76      -0.284     0.32
neuroticism   -0.18         -0.23        -0.22       0.812    -0.12
openness       0.15          0.19         0.22      -0.086     0.60

Item by scale correlations:
 corrected for item overlap and scale reliability
          agree conscientious extraversion neuroticism openness
A1        -0.40         -0.06        -0.11        0.14    -0.14
A2         0.67          0.23         0.40       -0.07     0.17
A3         0.70          0.22         0.48       -0.11     0.17
A4         0.49          0.29         0.30       -0.14     0.01
A5         0.62          0.23         0.55       -0.23     0.18
C1         0.13          0.53         0.19       -0.08     0.28
C2         0.21          0.61         0.17        0.00     0.20
C3         0.21          0.54         0.14       -0.09     0.08
C4        -0.24         -0.66        -0.23        0.31    -0.23
C5        -0.26         -0.59        -0.29        0.36    -0.10
E1        -0.30         -0.06        -0.59        0.11    -0.16
E2        -0.39         -0.25        -0.70        0.34    -0.15
E3         0.44          0.20         0.60       -0.10     0.37
E4         0.51          0.23         0.68       -0.22     0.04
E5         0.34          0.40         0.55       -0.10     0.31
N1        -0.22         -0.21        -0.11        0.76    -0.12
N2        -0.22         -0.19        -0.12        0.74    -0.06
N3        -0.14         -0.20        -0.14        0.74    -0.03
N4        -0.22         -0.30        -0.39        0.62    -0.02
N5        -0.04         -0.14        -0.19        0.55    -0.18
O1         0.16          0.20         0.31       -0.09     0.52
O2        -0.01         -0.18        -0.07        0.19    -0.45
O3         0.26          0.20         0.42       -0.07     0.61
O4         0.06         -0.02        -0.10        0.21     0.32
O5        -0.09         -0.14        -0.11        0.11    -0.53
gender     0.25          0.11         0.12        0.14    -0.07
education  0.06          0.03         0.01       -0.06     0.13
age        0.22          0.14         0.07       -0.13     0.10

 

コメント