Error with reshape of number of elements

Hello,

I'm working with multi-class classification by using function fitcecoc.

Featureste is a 1063x22 matrix while Wte is a 1063x1 vector.

Running the code shown below, a message error is given, saying: 'To RESHAPE the number of elements must not change'

Error is given at the following line:

reshape(max(PosteriorRegion,[ ],2),size(x1Grid,1),size(x1Grid,2)))

How can I fix it? Thank you in advance

Mdl = fitcecoc(Featureste,Wte,'Learners',t,'FitPosterior',1,...
    'ClassNames',{'1','2','3'},...
    'Verbose',2)

%Cross-validate Mdl using 10-fold cross-validation.
CVMdl = crossval(Mdl);

%CVMdl is a ClassificationPartitionedECOC cross-validated ECOC classifier.
%Estimate the classification error.
loss = kfoldLoss(CVMdl); 

%Define a grid of values in the observed predictor space. 
%Predict the posterior probabilities for each instance in the grid.
xMax = max(Featureste);
xMin = min(Featureste);

x1Pts = linspace(xMin(1),xMax(1));
x2Pts = linspace(xMin(2),xMax(2));
[x1Grid,x2Grid] = meshgrid(x1Pts,x2Pts);

[~,~,~,PosteriorRegion] = predict(Mdl,[x1Grid(:),x2Grid(:)]);

%For each coordinate on the grid
%plot the maximum class posterior probability among all classes.
figure
contourf(x1Grid,x2Grid,...
        reshape(max(PosteriorRegion,[],2),size(x1Grid,1),size(x1Grid,2)))

hold on
gh = gscatter(Featureste(:,1),Featureste(:,2),Wte,'krk','*xd',8)
gh(2).LineWidth = 2;
gh(3).LineWidth = 2;
hold off

Answer 1

[Cris LaPierre]

Here's my guess. You create x1Pts and x2Pts without specifying the number of elements you want. When you do that, you get a vector with length 100. See this comment from the documentation page:

• y = linspace(x1,x2) returns a row vector of 100 evenly spaced points between x1 and x2.

Therefore, x1Grid and x2Grid and 100x100.

I have no idea what the size of PosteriorRange is, but max(PosteriorRegion,[],2) must NOT be a 10000 x 1 (size(x1Grid,1) = size(x1Grid,2) = 100).

[Fed]

Thank you for your answer.

max(PosteriorRegion,[],2) has 1063 elements, so I have understood that in order to reshape in the correct way, the product between size(x1Grid,1) and size(x2Grid,2) must give 1063 elements as well.

I tried to write

x1Pts = linspace(xMin(1),xMax(1),1);
x2Pts = linspace(xMin(2),xMax(2),1063);

but when I run the code a message error says that

Z must be at least a 2x2 matrix 

1063 is a prime number... is there another way to obtain 1063 elements without putting in

x1Pts = linspace(X,Y,n);

n = 1?

Thank you for your help

[Cris LaPierre]

A countour plot must be given a matrix for input Z. You can't create a contour with a vector. Without knowing anything about your data, I can't tell you how you should reshape your data to see what you want to see.

Consider adding a row of NaN to make your vector 1064 so you can reshape it to a matrix.

[Fed]

Thank you!