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

1 次查看（过去 30 天）
MinChul Park
MinChul Park 2022-6-21
回答： Adam Danz 2022-6-23
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 个评论
显示 4更早的评论
MinChul Park
MinChul Park 2022-6-22
Thanks for that reply Adam! 😀
I've figured everything out and it works well.
Adam Danz
Adam Danz 2022-6-23
Glad it worked out! I'll move my comment to the answers section so your quesitons is moved from the unanswered questions queue.

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Adam Danz
Adam Danz 2022-6-23
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)

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Antennas, Microphones, and Sonar Transducers 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by