Я работаю над сценарием R, в котором я хочу сравнить различные наборы инструментов для нейровизуализации. Я хочу вычислить стандартное отклонение по voxels (формат файла nifit).
Мы отсканировали 6 объектов 50 раз, поэтому у меня есть 3D-массив и столбец для 50 сканирований (4D-массив). Я подумал, что могу использовать функцию apply для вычисления SD по строкам. Я понятия не имею, почему это не работает.
Вот мой код:
VBMsd <- apply(VBMimg, 1, function(x) sd(as.vector(x)))
VBMimg - массив с размерами 121, 145, 121, 50.
Я хочу вычислить sd по четвертому измерению (время сканирования), чтобы увидеть, есть ли разница во времени (x1-x2-x3-x4-x5..., y1-y2-y3..., z1-z2-z3 etc...), поэтому в основном sd по x,y,z координатам. Есть идеи?