Changeset 1073

Timestamp:

May 15, 2009, 2:11:37 PM (14 years ago)

Author:

Jari Häkkinen

Message:

Addresses #118. q=max(10,min(100,target_size/10))

File:

• plugins/base2/net.sf.basedb.normalizers/trunk/src/c++/bin/qQN.cc (modified) (1 diff)

plugins/base2/net.sf.basedb.normalizers/trunk/src/c++/bin/qQN.cc

@@ -127 +127 @@
   std::transform(data_iterator(m.begin()), data_iterator(m.end()),
                  data_iterator(m.begin()), theplu::yat::utility::Log<double>());
-  qQuantileNormalizer qqn(target.begin(), target.end(), 100);
+  // q = min(100,target_size/10) but no smaller than 10
+  unsigned int q=target.size()/10;
+  q = q>100 ? 100 : ( q<10 ? 10 : q );
+  qQuantileNormalizer qqn(target.begin(), target.end(), q);
   ColumnNormalizer<qQuantileNormalizer> cn(qqn);
   MatrixWeighted result(m.rows(),m.columns());