Issues with splitting a matrix into many smaller matrix

I am trying to group my origional matrix (size 11250 x 2) into many smaller matricies, grouped by the value in the first column.

I then want to take the mean and standard deviation of these values to report on the data that I have recieved.

PIN

Unrecognized function or variable 'PIN'.

PIN_SORT = sort(PIN.', 2).'

n = 1
for i = -4.5200:0.04: 0.7200
    M{n} = PIN_SORT(find(PIN_SORT(:,1)== i),:);
    %disp(M{'i'})
    %disp(i)
    n = n+1
end
avg{n} = mean(M{n})
stddev{n} = std(M{n})

The code does what I want it to do, for certain values, and will leave an empty spot in the array for the rest.

The two comments were attepmts at trouble shooting, but alas to no success.

Is there a glaringly obvious flaw in my approach or in my code, or is there an oversight on my behalf.

[Stephen23]

"Is there a glaringly obvious flaw in my approach or in my code, or is there an oversight on my behalf."

https://www.mathworks.com/matlabcentral/answers/57444-faq-why-is-0-3-0-2-0-1-not-equal-to-zero

https://www.mathworks.com/matlabcentral/answers/316889-why-two-equal-numbers-are-not-equal

http://www.mathworks.com/matlabcentral/answers/102419-how-do-i-determine-if-the-error-in-my-answer-is-the-result-of-round-off-error-or-a-bug

http://www.mathworks.com/matlabcentral/answers/140656-about-numerical-difficulty-of-equality-comparison-operation-between-two-double-real-numbers-in-matlab

This is worth reading as well:

https://docs.oracle.com/cd/E19957-01/806-3568/ncg_goldberg.html

Answer 1

[Walter Roberson]

for i = -4.5200:0.04: 0.7200

The i are floating point values.

    M{n} = PIN_SORT(find(PIN_SORT(:,1)== i),:);

You are comparing PIN_SORT(:,1) to a floating point value for bit-for-bit equality. Most values will not be found.

MATLAB does not calculate in decimal. The value -4.5200 cannot be exactly represented in binary double precision. The value 0.04 cannot be exactly represented in binary double precision. The not-exactly-0.04 will be repeatedly added to the original not-exiactly -4.5200 with the representation error for 0.04 accumulating each time.

You should not compare floating point numbers for equality -- except for cases where the values being compared are drawn from the same source. For example, it is valid to test

    PIN_SORT(:,1) == min(PIN_SORT(:,1))

because min() will pull out a bit-for-bit copy of one of the values, and you can be sure that bit-for-bit testing will work in that case (provided the array is not all NaN)

[Walter Roberson]

T1 = -4.5200:0.04: 0.7200;
T2 = (-452:4:72)/100;
mean(T1 == T2)

ans = 0.3030

Only about 30% of the entries in -4.5200:0.04:0.7200 will be bit-for-bit identical to what you would get if you calculated the corresponding integer divided by 100. The cumulative error on the 0.04 is a problem

You should use ismembertol or discretize

Answer 2

[KSSV]

You need to load the data from csv file into MATLAB first.

T = readtable('https://in.mathworks.com/matlabcentral/answers/uploaded_files/1455237/PIN.CSV');
PIN = T.(4)  % pick your column from table T 

PIN = 11252×1

NaN NaN -0.0021 -0.0020 -0.0020 -0.0020 -0.0020 -0.0020 -0.0020 -0.0020

[Lachlan Stephens]

Thanks KSSV, I have done that (just posts the issue in the online compiler), I have attached the fiel so you can see the issue.