本日はR言語で仮想データを実際に解析しながら、混合モデル解析の手続きを紹介します。
第一弾の今回は、もっとも単純な一要因のデザインを考えます。
仮想データは、次のような実験デザインを想定しています。
調べたい要因(Condition) 条件a vs. 条件b
従属変数 反応時間(RT)
被験者(Subject) 12名
刺激セット数(Item) 12セット
要因配置(listnumber) ラテン方格デザイン(4リストを作成)
データは以下の通りです。また以下では解説を黒字、R内でのコマンドを赤字、出力結果を青字で示します。
Subject | Item | Condition | RT | listnumber |
1 | 1 | a | 363.5 | 1 |
1 | 2 | a | 362.75 | 1 |
1 | 3 | a | 360 | 1 |
1 | 4 | a | 355.25 | 1 |
1 | 5 | a | 357.5 | 1 |
1 | 6 | a | 358.75 | 1 |
1 | 7 | b | 412.25 | 1 |
1 | 8 | b | 420.5 | 1 |
1 | 9 | b | 421.75 | 1 |
1 | 10 | b | 423 | 1 |
1 | 11 | b | 434.75 | 1 |
1 | 12 | b | 424 | 1 |
2 | 1 | b | 410.5 | 2 |
2 | 2 | b | 400.25 | 2 |
2 | 3 | b | 407 | 2 |
2 | 4 | b | 408.25 | 2 |
2 | 5 | b | 420.5 | 2 |
2 | 6 | b | 429.25 | 2 |
2 | 7 | a | 386.25 | 2 |
2 | 8 | a | 378.5 | 2 |
2 | 9 | a | 382.75 | 2 |
2 | 10 | a | 389 | 2 |
2 | 11 | a | 376.25 | 2 |
2 | 12 | a | 391.5 | 2 |
3 | 1 | a | 374.5 | 3 |
3 | 2 | b | 412.25 | 3 |
3 | 3 | a | 364 | 3 |
3 | 4 | b | 435.25 | 3 |
3 | 5 | a | 366.5 | 3 |
3 | 6 | b | 423.25 | 3 |
3 | 7 | a | 376.25 | 3 |
3 | 8 | b | 448.5 | 3 |
3 | 9 | a | 394.75 | 3 |
3 | 10 | b | 441 | 3 |
3 | 11 | a | 392.25 | 3 |
3 | 12 | b | 444 | 3 |
4 | 1 | b | 422.5 | 4 |
4 | 2 | a | 367.75 | 4 |
4 | 3 | b | 419 | 4 |
4 | 4 | a | 385.25 | 4 |
4 | 5 | b | 446.5 | 4 |
4 | 6 | a | 378.75 | 4 |
4 | 7 | b | 445.25 | 4 |
4 | 8 | a | 383.5 | 4 |
4 | 9 | b | 443.75 | 4 |
4 | 10 | a | 400 | 4 |
4 | 11 | b | 444.75 | 4 |
4 | 12 | a | 392.5 | 4 |
5 | 1 | a | 382.5 | 1 |
5 | 2 | a | 386.75 | 1 |
5 | 3 | a | 381 | 1 |
5 | 4 | a | 379.25 | 1 |
5 | 5 | a | 391.5 | 1 |
5 | 6 | a | 380.75 | 1 |
5 | 7 | b | 436.25 | 1 |
5 | 8 | b | 443.5 | 1 |
5 | 9 | b | 435.75 | 1 |
5 | 10 | b | 435 | 1 |
5 | 11 | b | 440.75 | 1 |
5 | 12 | b | 445 | 1 |
6 | 1 | b | 443.5 | 2 |
6 | 2 | b | 428.25 | 2 |
6 | 3 | b | 434 | 2 |
6 | 4 | b | 442.25 | 2 |
6 | 5 | b | 444.5 | 2 |
6 | 6 | b | 444.25 | 2 |
6 | 7 | a | 403.25 | 2 |
6 | 8 | a | 397.5 | 2 |
6 | 9 | a | 411.75 | 2 |
6 | 10 | a | 405 | 2 |
6 | 11 | a | 413.25 | 2 |
6 | 12 | a | 410.5 | 2 |
7 | 1 | a | 396.5 | 3 |
7 | 2 | b | 449.25 | 3 |
7 | 3 | a | 405 | 3 |
7 | 4 | b | 450.25 | 3 |
7 | 5 | a | 391.5 | 3 |
7 | 6 | b | 443.25 | 3 |
7 | 7 | a | 412.25 | 3 |
7 | 8 | b | 455.5 | 3 |
7 | 9 | a | 409.75 | 3 |
7 | 10 | b | 470 | 3 |
7 | 11 | a | 414.25 | 3 |
7 | 12 | b | 462 | 3 |
8 | 1 | b | 454.5 | 4 |
8 | 2 | a | 400.75 | 4 |
8 | 3 | b | 436 | 4 |
8 | 4 | a | 398.25 | 4 |
8 | 5 | b | 448.5 | 4 |
8 | 6 | a | 409.75 | 4 |
8 | 7 | b | 448.25 | 4 |
8 | 8 | a | 420.5 | 4 |
8 | 9 | b | 466.75 | 4 |
8 | 10 | a | 415 | 4 |
8 | 11 | b | 475.75 | 4 |
8 | 12 | a | 422.5 | 4 |
9 | 1 | a | 396.5 | 1 |
9 | 2 | a | 409.75 | 1 |
9 | 3 | a | 401 | 1 |
9 | 4 | a | 415.25 | 1 |
9 | 5 | a | 401.5 | 1 |
9 | 6 | a | 414.75 | 1 |
9 | 7 | b | 455.25 | 1 |
9 | 8 | b | 470.5 | 1 |
9 | 9 | b | 458.75 | 1 |
9 | 10 | b | 465 | 1 |
9 | 11 | b | 463.75 | 1 |
9 | 12 | b | 481 | 1 |
10 | 1 | b | 460.5 | 2 |
10 | 2 | b | 464.25 | 2 |
10 | 3 | b | 458 | 2 |
10 | 4 | b | 465.25 | 2 |
10 | 5 | b | 479.5 | 2 |
10 | 6 | b | 461.25 | 2 |
10 | 7 | a | 414.25 | 2 |
10 | 8 | a | 432.5 | 2 |
10 | 9 | a | 434.75 | 2 |
10 | 10 | a | 425 | 2 |
10 | 11 | a | 428.25 | 2 |
10 | 12 | a | 425.5 | 2 |
11 | 1 | a | 418.5 | 3 |
11 | 2 | b | 465.25 | 3 |
11 | 3 | a | 421 | 3 |
11 | 4 | b | 469.25 | 3 |
11 | 5 | a | 416.5 | 3 |
11 | 6 | b | 478.25 | 3 |
11 | 7 | a | 436.25 | 3 |
11 | 8 | b | 474.5 | 3 |
11 | 9 | a | 436.75 | 3 |
11 | 10 | b | 477 | 3 |
11 | 11 | a | 436.25 | 3 |
11 | 12 | b | 495 | 3 |
12 | 1 | b | 479.5 | 4 |
12 | 2 | a | 424.75 | 4 |
12 | 3 | b | 459 | 4 |
12 | 4 | a | 421.25 | 4 |
12 | 5 | b | 481.5 | 4 |
12 | 6 | a | 417.75 | 4 |
12 | 7 | b | 482.25 | 4 |
12 | 8 | a | 435.5 | 4 |
12 | 9 | b | 480.75 | 4 |
12 | 10 | a | 439 | 4 |
12 | 11 | b | 491.75 | 4 |
12 | 12 | a | 441.5 | 4 |
記述が長くなりそうなので、いくつかに分割して掲載することにします。
(2)へ続く
0 件のコメント:
コメントを投稿