%PRTEST test routine for PRTools commands
%
% This functions runs countless basic Prtools operations on basic
% Prtools objects. Set variable 'item' for specific choices.
% Possible 'item' values are: dataset, labtype, mapping, fixed, classifiers
if ~exist('item'), item = 'prdataset'; end
echo on
switch item
case 'prdataset'
disp('dataset definition')
r = rand(10,2);
lab = genlab([3,7])
a = prdataset(r,lab)
lab = genlab([3,7],['apple1';'apple2'])
a = prdataset(r,lab)
lab = genlab([3,7],{'pear';'apple'})
a = prdataset(r,lab)
struct(a)
getdata(a)
a.data
get(a,'data')
getlabels(a)
a.labels
get(a,'labels')
getnlab(a)
a.nlab
get(a,'nlab')
getfeatlab(a)
a.featlab
get(a,'featlab')
a = set(a,'featlab',{'feat1';'f2'})
get(a,'featlab')
getprior(a)
a.prior
get(a,'prior')
getlablist(a)
a.lablist
get(a,'lablist')
getlabtype(a)
a.labtype
get(a,'labtype')
size(a)
getobjsize(a)
a.objsize
get(a,'objsize')
getfeatsize(a)
a.featsize
get(a,'featsize')
getident(a)
a.ident
get(a,'ident')
getversion(a)
a.version
a.version{1}
a.version{2}
get(a,'version')
a = setname(a,'test of dataset definition')
getname(a)
a.name
get(a,'name')
a = setuser(a,'some rubbish')
getuser(a)
a.user
get(a,'user')
case 'labtype'
a = gendath(5)
a = setlabels(a,genlab([3 2 3 2],[4 3 1 2]'))
getlablist(a)'
getlabels(a)'
gettargets(a)'
b = setlabtype(a,'soft')
getlablist(b)'
getlabels(b)'
gettargets(b)'
c = setlabtype(a,'targets')
getlablist(c)'
getlabels(c)'
gettargets(c)'
a = setlabtype(b,'crisp')
getlablist(a)'
getlabels(a)'
gettargets(a)
a = setlabtype(b,'targets')
getlablist(a)'
getlabels(a)'
gettargets(a)
a = setlabtype(c,'crisp')
getlablist(a)'
getlabels(a)'
gettargets(a)
a = setlabtype(c,'soft')
getlablist(a)'
getlabels(a)'
gettargets(a)
case 'prmapping'
a = prdataset(rand(10,2),genlab([5,5]))
w = affine(rand(2,2),[1 0])
struct(w)
getmapping_file(w)
w.mapping_file
get(w,'mapping_file')
getmapping_type(w)
w.mapping_type
get(w,'mapping_type')
getdata(w)
w.data
get(w,'data')
getlabels(w)
w.labels
get(w,'labels')
getsize_in(w)
w.size_in
get(w,'size_in')
getsize_out(w)
w.size_out
get(w,'size_out')
getsize(w)
w.size
get(w,'size')
getscale(w)
w.scale
get(w,'scale')
getout_conv(w)
w.out_conv
get(w,'out_conv')
getname(w)
w.name
get(w,'name')
getuser(w)
w.user
get(w,'user')
v = [w w w]
v = [w; w; w]
case 'fixed'
a = prdataset(rand(5,2),[1 1 2 2 2]');
disp('Testing fixed mappings')
disp('Dataset')
+a
disp('feature selection')
b = [a a a];
+b
w = cmapm(size(b,2),[1 3 4 6])
+(b*w)
disp('sigmoid')
w = sigm([],2)
+(a*w)
disp('exp')
w = cmapm(2,'exp')
+(a*w)
disp('nexp')
w = cmapm(2,'nexp')
+(a*w)
disp('log')
w = cmapm(2,'log')
+(a*w)
disp('randrot')
w = cmapm(2,'randrot')
+(a*w)
disp('rot')
w = cmapm([1 1;1 -1],'rot')
figure
scatterd(a);
figure
scatterd(a*w)
disp('shift')
w = cmapm([100 50],'shift');
disp([mean(a);mean(a*w)])
w = cmapm([100,0.01],'scale');
disp([std(a);std(a*w)])
case 'classifiers'
W = {fisherc, quadrc, qdc, ldc, udc, nmc, nmsc, knnc, knnc([],1),parzenc, ...
klldc, pcldc, svc, loglc, treec, neurc, lmnc, bpxnc}
echo off
for j=1:length(W)
w = W{j};
prtestc
end
case 'mappings'
W = {cmapm, scalem, klm, pca, klms, reducm, fisherm, nlfisherm, kernelm, ...
sigm, invsigm, distm, proxm, lmnm, normm, svm}
case 'combining'
echo off
W = {maxc, minc, meanc, prodc, medianc, votec}
a = gendath;
[b,c] = gendat(a,0.2);
w = [fisherc, nmc, nmsc, qdc, knnc([],3)];
for j=1:length(W)
v1 = b*w*W{j};
v2 = b*(w*W{j});
v3 = b*w*classc*W{j};
v4 = b*(w*classc*W{j});
e1 = c*{v1 v2 v3 v4}*testc;
e2 = c*(b*w*W{j})*testc;
e3 = c*(b*w*classc*W{j})*testc;
e4 = c*(b*w*W{j}*classc)*testc;
disp(getname(W{j}))
disp([e1{:} e2 e3 e4])
end
end
echo off