[R] Collapsing across trials

Pedro Alcocer pealco at gmail.com
Sun Dec 17 05:20:19 CET 2006


Hello,

My ultimate goal is a repeated measures (mixed model) ANOVA, however,
my present question is about how to reorganize my data into the format
that the ANOVA commands expect. In particular, how to collapse across
trials. (I am using the tutorial at
[http://personality-project.org/r/r.anova.html] for the mixed model
ANOVA)

The data I am using looks like this. A subject sees 10 trials per
condition and there are 2 conditions. I want the average of all the
reaction times (RTs) from a subject looking a one condition. I also
want to retain List in the final output.

Subj	List	Condition	RT
2	1	C		338
2	1	C		227
2	1	C		430
2	1	C		621
2	1	C		255
2	1	C		348
2	1	C		280
2	1	C		356
2	1	C		272
2	1	C		346
3	2	C		489
3	2	C		426
3	2	C		352
3	2	C		351
3	2	C		349
3	2	C		403
3	2	C		336
3	2	C		278
3	2	C		365
3	2	C		271
2	1	D		360
2	1	D		374
2	1	D		326
2	1	D		363
2	1	D		290
2	1	D		458
2	1	D		295
2	1	D		362
2	1	D		285
2	1	D		277
3	2	D		354
3	2	D		352
3	2	D		362
3	2	D		360
3	2	D		334
3	2	D		365
3	2	D		335
3	2	D		391
3	2	D		272
3	2	D		618

The result should look like this:
2	1	C		(AVG)
3	2	C		(AVG)
2	1	D		(AVG)
3	2	D		(AVG)

Where (AVG) is the average of the 10 trials.

The above is a simplified case. How can I do this with multiple RT
measurements per subject? In other words, the above, but with more
than one RT column per subject.

Resulting in:

Subj	List	Condition	RT1		RT2		RT3		RT4		RT5
2	1	C		(AVG1)	(AVG2)	(AVG2)	(AVG2)	(AVG2)
3	2	C		(AVG1)	(AVG2)	(AVG2)	(AVG2)	(AVG2)
2	1	D		(AVG1)	(AVG2)	(AVG2)	(AVG2)	(AVG2)
3	2	D		(AVG1)	(AVG2)	(AVG2)	(AVG2)	(AVG2)

I've come across the apply and aggregate functions in online
documentation, and I have the suspicion that they may be called for
here, but their application isn't clear to me. I am fairly new to R.

(I am using R 2.4.0 on Mac OS X.)

I'd appreciate any insights.

Pedro Alcocer
University of Florida



More information about the R-help mailing list