2012年12月17日月曜日

R言語を用いた混合モデル解析の紹介その1:1要因デザインの場合(1)


本日はR言語で仮想データを実際に解析しながら、混合モデル解析の手続きを紹介します。
第一弾の今回は、もっとも単純な一要因のデザインを考えます。

仮想データは、次のような実験デザインを想定しています。

調べたい要因(Condition) 条件a vs. 条件b
従属変数 反応時間(RT)
被験者(Subject) 12名
刺激セット数(Item) 12セット
要因配置(listnumber) ラテン方格デザイン(4リストを作成)

データは以下の通りです。また以下では解説を黒字、R内でのコマンドを赤字、出力結果を青字で示します。

SubjectItemConditionRTlistnumber
11a363.51
12a362.751
13a3601
14a355.251
15a357.51
16a358.751
17b412.251
18b420.51
19b421.751
110b4231
111b434.751
112b4241
21b410.52
22b400.252
23b4072
24b408.252
25b420.52
26b429.252
27a386.252
28a378.52
29a382.752
210a3892
211a376.252
212a391.52
31a374.53
32b412.253
33a3643
34b435.253
35a366.53
36b423.253
37a376.253
38b448.53
39a394.753
310b4413
311a392.253
312b4443
41b422.54
42a367.754
43b4194
44a385.254
45b446.54
46a378.754
47b445.254
48a383.54
49b443.754
410a4004
411b444.754
412a392.54
51a382.51
52a386.751
53a3811
54a379.251
55a391.51
56a380.751
57b436.251
58b443.51
59b435.751
510b4351
511b440.751
512b4451
61b443.52
62b428.252
63b4342
64b442.252
65b444.52
66b444.252
67a403.252
68a397.52
69a411.752
610a4052
611a413.252
612a410.52
71a396.53
72b449.253
73a4053
74b450.253
75a391.53
76b443.253
77a412.253
78b455.53
79a409.753
710b4703
711a414.253
712b4623
81b454.54
82a400.754
83b4364
84a398.254
85b448.54
86a409.754
87b448.254
88a420.54
89b466.754
810a4154
811b475.754
812a422.54
91a396.51
92a409.751
93a4011
94a415.251
95a401.51
96a414.751
97b455.251
98b470.51
99b458.751
910b4651
911b463.751
912b4811
101b460.52
102b464.252
103b4582
104b465.252
105b479.52
106b461.252
107a414.252
108a432.52
109a434.752
1010a4252
1011a428.252
1012a425.52
111a418.53
112b465.253
113a4213
114b469.253
115a416.53
116b478.253
117a436.253
118b474.53
119a436.753
1110b4773
1111a436.253
1112b4953
121b479.54
122a424.754
123b4594
124a421.254
125b481.54
126a417.754
127b482.254
128a435.54
129b480.754
1210a4394
1211b491.754
1212a441.54



記述が長くなりそうなので、いくつかに分割して掲載することにします。
(2)へ続く

0 件のコメント:

コメントを投稿