Changeset 151 for prdatasets/flowcyto.m

Timestamp:

01/11/20 00:14:42 (5 years ago)

Author:

bduin

Message:

 

File:

• prdatasets/flowcyto.m (modified) (1 diff)

prdatasets/flowcyto.m

@@ -25 +25 @@
 
 if nargin < 1 || isempty(n)
-  n = 1;
+  n = 0;
 end
 
-if any(n<1) || any(n>4)
+if any(n<0) || any(n>4)
   error('Dataset number should be larger than 0 and smaller than 5')
 end
 
-d = pr_getdata('http://prtools.tudelft.nl/prdatasets/FlowCyto_sort.mat',5);
-
-a = cell(1,numel(n));
-for j=1:numel(n)
-  x = prdataset(d{n(j)});
-  x = setname(x(:,2:255),['Flow Cytometry ' num2str(n(j))]);
-  desc = ['These datasets are based on 612 FL3-A DNA flowcytometer histograms' ...
-  'from breast cancer tissues in 256 resolution. The initial data have been' ...
-  'acquired by M. Nap and N. van Rodijnen of the Atrium Medical Center in' ...
-  'Heerlen, The Netherlands, during 2000-2004, using tubes 3-6 of a DACO' ...
-  'Galaxy flowcytometer (N = 1-4). Histograms are labeled in 3 classes:' ...
-  'aneuploid (335 patients), diploid (131) and tetraploid (146). The first' ...
-  'and the last bin (1 and 256) are set to 0, as they contsin noise. After' ...
-  'that histograms are normalized (sum to one) resulting in a dataset with' ...
-  '254 features.'];
-  ref = {['']};
-  x = setprior(x,getprior(x,0));
-  x = setuser(x,desc,'desc');
-  x = setuser(x,ref,'ref');
-  a{j} = x;
+a = pr_loadmatfile(['flowcyto_' num2str(n)]);
+if isempty(a)
+  d = pr_getdata('http://prtools.tudelft.nl/prdatasets/FlowCyto_sort.mat',5);
+  matfile = fullfile(fullfile(fileparts(which(mfilename)),'data'),'flowcyto');
+  b = [];
+  for i=1:4
+    a = setname(prdataset(d{i}),['Flow Cytometry Tube_' num2str(i+2)]);
+    save([matfile '_' num2str(i)],'a');
+    b = [b a];
+  end
+  a = setname(b,'Flow Cytometry');
+  save([matfile '_0'],'a');
+  a = pr_loadmatfile(['flowcyto_' num2str(n)]);
 end
-
-if numel(a) == 1
-  a = a{1};
-end
-
-