blob: 4c440db2833cd68eb8c7faeb99f60885a52d042a (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
|
function [skew, kurtosis] = pdf_stats(x,p)
% Define some useful bits
dx = x(2) - x(1);
% The expected value (From LibreTexts 4.1)
E = @(x) sum(x .* p) * dx;
% First moment (mean) (From LibreTexts 4.3)
mu = E(x)
% Second moment (variance) (From LibreTexts 4.3)
sigma = sqrt(E((x - mu).^2))
% Third moment (skewness) (From LibreTexts 4.4)
skew = E(((x - mu)/sigma).^3)
% Fourth moment (excess kurtosis) (From LibreTexts 4.4)
kurtosis = E(((x - mu)/sigma).^4) - 3
end
|
