Standard deviation in 3D array from 1st to nth number array

2 views (last 30 days)
Hello I am trying to calculate residual noise from auditory brainstem responses (ABR).
The 3D array that I'm working in is in the size [2 x 493 x 6000] double which is a reflection of [Nchannels x Ntime x Nepochs].
I have 2 recording channels and each channel consists of 493 time data points and there are 6000 trials (or epochs).
Will call our 3D array ABR.
To calculate residual noise I have to:
  1. calculate the std of the data points (in rows) from the 1st data point of epoch no.1 (which is given by ABR(:,1,1)) to the last data point of epoch no. X (which is given by ABR(:, end, X)).
  2. The issue that I'm facing is how to specify calculation of std from 1st epoch to nth epoch.
  3. I've tried std(ABR,0,2,[ABR(:,1,1), ABR(:,end,1)]); but this is not the answer.
  4. I need to calcuate 6000 of these in the following format...
  5. std(ABR,0,2,[ABR(:,1,1), ABR(:,end,1)]);
  6. std(ABR,0,2,[ABR(:,1,1), ABR(:,end,2)]);
  7. std(ABR,0,2,[ABR(:,1,1), ABR(:,end,3)]);
  8. std(ABR,0,2,[ABR(:,1,1), ABR(:,end,4)]);
  9. std(ABR,0,2,[ABR(:,1,1), ABR(:,end,5)]); ... and so on until the last last number becomes 6000.
  10. In this case ABR is the data set, 0 is the weighting, 2 means calculation in rows and [] specifies what arrays (or epochs) to calculate std from.
This is what I am tryin to do graphically - using excel.
As you can see each std data point is the std result of all the epochs upto and including the position of the std.
What is the best way to do this?
  6 Comments
MinChul Park
MinChul Park on 22 Jun 2022
Thanks for that reply Adam! 😀
I've figured everything out and it works well.
Adam Danz
Adam Danz on 23 Jun 2022
Glad it worked out! I'll move my comment to the answers section so your quesitons is moved from the unanswered questions queue.

Sign in to comment.

Accepted Answer

Adam Danz
Adam Danz on 23 Jun 2022
To cumulatively compute the std across dimension 3,
data = rand(6, 493, 6000);
RN = zeros(size(data,[1,3]));
for i = 1 : size(data,3)
RN(:,i) = std(data(:, :, 1:i), 0, [2,3]);
end
size(RN)

More Answers (0)

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!