Пост
Даркин должен иметь list()
вместо c()
. Тогда я думаю, что это должно работать. Вот полный пример:
set.seed(12345)
data = expand.grid(trial=1:10, stim1=1:5, stim2=1:3)
data = data.frame(data, response=rnorm(nrow(data)))
with(data, tapply(response, list(stim1,stim2), mean))
Выход:
1 2 3
1 -0.13294415 0.27326245 -0.11120045
2 0.28597776 0.02338804 0.21280916
3 0.08338741 0.44086561 -0.08682628
4 0.72432003 0.84250712 0.28383378
5 -0.06290978 -0.02588252 -0.36364019
(Кстати, прошу прощения за новый ответ. У меня недостаточно репутации, чтобы комментировать, а изменение слишком маленькое, чтобы применить его непосредственно к предыдущему сообщению.)