ugw/data/preprocess.m

   1 #!/usr/bin/octave -q
   2 ## -*-octave-*-
   3
   4 datafiles = argv();
   5 quantiles = [0 0.5 0.90 0.99];
   6
   7 mat = zeros(0, length(quantiles));
   8 names = {};
   9
  10 for i=1:length(argv),
  11   x = load(datafiles{i})*1e6;
  12   q = quantile (x, quantiles);
  13   mat = [ mat; q' ];
  14   [d, n, e, v] = fileparts(datafiles{i});
  15   names{i} = n;
  16 end
  17
  18 [mat, ind] = sortrows(mat, 2);
  19 result = [names(ind)', num2cell(mat)];
  20 cell2csv("result", result, ' ', '"');