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

2 ビュー (過去 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 件のコメント
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 ExchangeSpectral Estimation についてさらに検索

Community Treasure Hunt

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

Start Hunting!

Translated by