MATLAB Central Newsreader - classify in 3D

classify in 3D

basalganglia Sch?necker — Sat, 31 Oct 2009 13:35:04 -0400

Hallo,

I have data similar to that of the fisheriris data in the classification demo on http://www.mathworks.com/products/statistics/demos.html?file=/products/demos/shipping/stats/classdemo.html
The demo is exactly what I want to do with my data. However, I use 3 variables instead of the 2 shown in the demo (3D position in space). The classification works fine but now the problem is how to visualize the classification "line" in 3D when applying the quadratic type of classification.

In the demo in with 2 variables, the curve is described by the following formula 'f':
K = coeff(1,2).const;
L = coeff(1,2).linear;
Q = coeff(1,2).quadratic;
f = sprintf('0 = %g+%g*x+%g*y+%g*x^2+%g*x.*y+%g*y.^2',...
K,L,Q(1,1),Q(1,2)+Q(2,1),Q(2,2));
h2 = ezplot(f,[4.5 8 2 4]);

How would I extend the formula for 3 variables ? It seems that the order of the coefficients provided by classify is not further documented in detail.

Much thanks, Thomas

Re: classify in 3D

Tom Lane — Mon, 02 Nov 2009 23:02:10 -0500

> I have data similar to that of the fisheriris data in the classification
> demo
. . .
> How would I extend the formula for 3 variables ? It seems that the order
> of the coefficients provided by classify is not further documented in
> detail.

Thomas, the example picks apart the coefficient arrays into separate scalars
so they can be conveniently used in a two-dimensional example. The intent is
that the coefficients are documented using the following matrix-oriented
expression given in the help:

0 < K + x*L + x*Q*x'

To plot this, the isosurface function may be useful. I don't know of an
easier way. Try this:

% Get 3 measurements from 2 iris species
load fisheriris
rows = 51:150;
X = meas(rows,1:3);
x = X(:,1); y = X(:,2); z = X(:,3);
s = species(rows);
a = 1:50; b = 51:100;

% Do quadratic classification and get coefficients of boundary
[c,err,post,logl,str] = classify(X,X,s,'quadratic');
K = str(1,2).const;
L = str(1,2).linear;
Q = str(1,2).quadratic;

% Plot the data and curve K + [x,y,z]*L + [x,y,z]*Q*[x,y,z]' = 0:
xv = linspace(4.5,8,10); % vectors to cover the range of each column
yv = linspace(2,4,10);
zv = linspace(2.5,7,10);
[xx,yy,zz] = meshgrid(xv,yv,zv);
f = @(x,y,z) K + [x y z]*L + sum([x y z] .* ([x y z]*Q), 2);
v = f(xx(:),yy(:),zz(:));
v = reshape(v,size(xx));
plot3(x(a),y(a),z(a),'rv', x(b),y(b),z(b),'b^');
hold on
isosurface(xx,yy,zz,v,0);
hold off

-- Tom

Re: classify in 3D

Thomas schoenecker — Tue, 03 Nov 2009 00:37:02 -0500

Hallo Tom,

Your help solved my problem. Thank You very much.

Thomas
"Tom Lane" <tlane@mathworks.com> wrote in message <hcnodi$aob$1@fred.mathworks.com>...
> > I have data similar to that of the fisheriris data in the classification
> > demo
> . . .
> > How would I extend the formula for 3 variables ? It seems that the order
> > of the coefficients provided by classify is not further documented in
> > detail.
>
> Thomas, the example picks apart the coefficient arrays into separate scalars
> so they can be conveniently used in a two-dimensional example. The intent is
> that the coefficients are documented using the following matrix-oriented
> expression given in the help:
>
> 0 < K + x*L + x*Q*x'
>
> To plot this, the isosurface function may be useful. I don't know of an
> easier way. Try this:
>
> % Get 3 measurements from 2 iris species
> load fisheriris
> rows = 51:150;
> X = meas(rows,1:3);
> x = X(:,1); y = X(:,2); z = X(:,3);
> s = species(rows);
> a = 1:50; b = 51:100;
>
> % Do quadratic classification and get coefficients of boundary
> [c,err,post,logl,str] = classify(X,X,s,'quadratic');
> K = str(1,2).const;
> L = str(1,2).linear;
> Q = str(1,2).quadratic;
>
> % Plot the data and curve K + [x,y,z]*L + [x,y,z]*Q*[x,y,z]' = 0:
> xv = linspace(4.5,8,10); % vectors to cover the range of each column
> yv = linspace(2,4,10);
> zv = linspace(2.5,7,10);
> [xx,yy,zz] = meshgrid(xv,yv,zv);
> f = @(x,y,z) K + [x y z]*L + sum([x y z] .* ([x y z]*Q), 2);
> v = f(xx(:),yy(:),zz(:));
> v = reshape(v,size(xx));
> plot3(x(a),y(a),z(a),'rv', x(b),y(b),z(b),'b^');
> hold on
> isosurface(xx,yy,zz,v,0);
> hold off
>
> -- Tom
>