how matlab determine skewness ?

16 ビュー (過去 30 日間)
crixus
crixus 2015 年 3 月 24 日
コメント済み: crixus 2015 年 3 月 25 日
hi all, i tabulated out my data residual and i do a boxplot on the residual, what i find is the median line is closer to the 75th percentile so by right it should be positively skewed and i use matlab function skewness(data) and got a positive value which validate the residual boxplot as positively skewed. However, if its positively skewed it should mean there are more observation on the positive side, but of all the residual that i have(124 of them) only 30+ of them is positive. I wondering if it still considered as positively skewed ? or is there something i missed out ? (and yes, my mean is smaller than median)
  1 件のコメント
crixus
crixus 2015 年 3 月 24 日

サインインしてコメントする。

採用された回答

Image Analyst
Image Analyst 2015 年 3 月 24 日
Please post a screenshot, and attach your data file and m-code to read it in and plot it, if you can.
Be aware that skewness is determined not only by how many data points are to the right and left of the mode but also how far away they are . So more points that are on the left but close to the mode, may not overwhelm a few points that are on the right but much farther away, giving an overall positive skewness even though more points are on the left. Here is the formula for an image:
% Get the skew.
skewness = sum((GLs - meanGL) .^ 3 .* pixelCounts) / ((numberOfPixels - 1) * stdDev^3);
  5 件のコメント
Image Analyst
Image Analyst 2015 年 3 月 24 日
I did it manually - not using skewness() function. I get -1.22.
s=load('matlab.mat')
array = s.err_sptb
referenceValue = 0;
numElements = length(array)
stdDev = std(array)
skewness = sum((array - referenceValue) .^ 3) / ((numElements - 1) * stdDev^3)
numElements =
123
stdDev =
345.66713081467
skewness =
-1.22297782162426
crixus
crixus 2015 年 3 月 25 日
thanks for your help ! really appreciate it.

サインインしてコメントする。

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeLogical についてさらに検索

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by