# How can I calculate Kurtosis of Sound data?

1 view (last 30 days)
Serhat Sayim on 3 Jan 2021
Commented: dpb on 3 Jan 2021
I have sound data of people who have asthma.I made separate plots windows thanks to the start and end indexes of wheezing and non-wheezing sounds. How can I calculate the kurtosis value of each wheeze or nonwheeze windows.
wheezing starting indexes
2239
19576
30338
46537
56438
69869
97067
107264
123725
134174
wheezing ending indexes
3239
24791
31829
49212
57662
72755
99823
108708
126971
136379

Star Strider on 3 Jan 2021
Not certain what you want.
Try this:
startidx = [2239
19576
30338
46537
56438
69869
97067
107264
123725
134174];
endidx = [3239
24791
31829
49212
57662
72755
99823
108708
126971
136379];
framelen = endidx - startidx;
fl_mean = mean(framelen);
fl_varv = var(framelen);
fl_kurt = kurtosis(framelen);
.
dpb on 3 Jan 2021
Kin=arrayfun(@(2239,3239)kurtosis(X(2239:3239)),iStart,iEnd);
That's not the code I gave you...just plug your variables in where I told you:
Kin=arrayfun(@(i1,i2)kurtosis(X(i1:i2)),iStart,iEnd);
Use whatever is your array/variable name for the data in place of "X". You plotted the data so you have the variable already, use it again.

dpb on 3 Jan 2021
Kin=arrayfun(@(i1,i2)kurtosis(X(i1:i2)),iStart,iEnd);
where X are your data and the two indexing arrays above for the "in group" sections.
The "out group" indices can be derived from the above in group values by adjustment of point numbers by one and adding the first, last points of the X vector. The same logic/functional form will then work for those sections.
Will leave as "exercise for Student" to consider the details...
##### 2 CommentsShowHide 1 older comment
dpb on 3 Jan 2021
Edited: dpb on 3 Jan 2021
That's exactly the code I gave you for in-between the two indices -- you have to do a minimal amount of adjustment of those two arrays to get the start-stop for the other segments from those values.

### Categories

Find more on Signal Attributes and Indexing in Help Center and File Exchange

### Community Treasure Hunt

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

Start Hunting!

Translated by