diff options
author | Erich Schubert <erich@debian.org> | 2013-10-29 20:02:37 +0100 |
---|---|---|
committer | Andrej Shadura <andrewsh@debian.org> | 2019-03-09 22:30:37 +0000 |
commit | ec7f409f6e795bbcc6f3c005687954e9475c600c (patch) | |
tree | fbf36c0ab791c556198b487ca40ae56ae5ab1ee5 /src/de/lmu/ifi/dbs/elki/datasource/filter/transform/GlobalPrincipalComponentAnalysisTransform.java | |
parent | 974d4cf6d54cadc06258039f2cd0515cc34aeac6 (diff) | |
parent | 8300861dc4c62c5567a4e654976072f854217544 (diff) |
Import Debian changes 0.6.0~beta2-1
elki (0.6.0~beta2-1) unstable; urgency=low
* New upstream beta release.
* 3DPC extension is not yet included.
Diffstat (limited to 'src/de/lmu/ifi/dbs/elki/datasource/filter/transform/GlobalPrincipalComponentAnalysisTransform.java')
-rw-r--r-- | src/de/lmu/ifi/dbs/elki/datasource/filter/transform/GlobalPrincipalComponentAnalysisTransform.java | 50 |
1 files changed, 27 insertions, 23 deletions
diff --git a/src/de/lmu/ifi/dbs/elki/datasource/filter/transform/GlobalPrincipalComponentAnalysisTransform.java b/src/de/lmu/ifi/dbs/elki/datasource/filter/transform/GlobalPrincipalComponentAnalysisTransform.java index 18537a8d..3b4193ad 100644 --- a/src/de/lmu/ifi/dbs/elki/datasource/filter/transform/GlobalPrincipalComponentAnalysisTransform.java +++ b/src/de/lmu/ifi/dbs/elki/datasource/filter/transform/GlobalPrincipalComponentAnalysisTransform.java @@ -4,7 +4,7 @@ package de.lmu.ifi.dbs.elki.datasource.filter.transform; This file is part of ELKI: Environment for Developing KDD-Applications Supported by Index-Structures - Copyright (C) 2012 + Copyright (C) 2013 Ludwig-Maximilians-Universität München Lehr- und Forschungseinheit für Datenbanksysteme ELKI Development Team @@ -38,6 +38,7 @@ import de.lmu.ifi.dbs.elki.math.linearalgebra.VMath; import de.lmu.ifi.dbs.elki.math.linearalgebra.pca.EigenPairFilter; import de.lmu.ifi.dbs.elki.math.linearalgebra.pca.PCAResult; import de.lmu.ifi.dbs.elki.math.linearalgebra.pca.PCARunner; +import de.lmu.ifi.dbs.elki.utilities.Alias; import de.lmu.ifi.dbs.elki.utilities.exceptions.AbortException; import de.lmu.ifi.dbs.elki.utilities.optionhandling.AbstractParameterizer; import de.lmu.ifi.dbs.elki.utilities.optionhandling.OptionID; @@ -47,12 +48,16 @@ import de.lmu.ifi.dbs.elki.utilities.optionhandling.parameters.ObjectParameter; /** * Apply principal component analysis to the data set. * - * TODO: add dimensionality reduction + * This process is also known as "Whitening transformation". + * + * If you want to also reduce dimensionality, set the + * {@link Parameterizer#FILTER_ID} parameter! * * @author Erich Schubert * * @param <O> Vector type */ +@Alias({ "whiten", "whitening", "pca" }) public class GlobalPrincipalComponentAnalysisTransform<O extends NumberVector<?>> extends AbstractVectorConversionFilter<O, O> { /** * Class logger. @@ -101,7 +106,7 @@ public class GlobalPrincipalComponentAnalysisTransform<O extends NumberVector<?> @Override protected boolean prepareStart(SimpleTypeInformation<O> in) { - if(!(in instanceof VectorFieldTypeInformation)) { + if (!(in instanceof VectorFieldTypeInformation)) { throw new AbortException("PCA can only applied to fixed dimensionality vectors"); } dim = ((VectorFieldTypeInformation<?>) in).getDimensionality(); @@ -121,31 +126,30 @@ public class GlobalPrincipalComponentAnalysisTransform<O extends NumberVector<?> SortedEigenPairs eps = pcares.getEigenPairs(); covmat = null; - if(filter == null) { + if (filter == null) { proj = new double[dim][dim]; - for(int d = 0; d < dim; d++) { + for (int d = 0; d < dim; d++) { EigenPair ep = eps.getEigenPair(d); double[] ev = ep.getEigenvector().getArrayRef(); double eval = Math.sqrt(ep.getEigenvalue()); // Fill weighted and transposed: - for(int i = 0; i < dim; i++) { + for (int i = 0; i < dim; i++) { proj[d][i] = ev[i] / eval; } } - } - else { + } else { List<EigenPair> axes = filter.filter(eps).getStrongEigenPairs(); final int pdim = axes.size(); // Projection dimensionality if (LOG.isVerbose()) { - LOG.verbose("Reducing dimensionality from "+dim+" to "+pdim+" via PCA."); + LOG.verbose("Reducing dimensionality from " + dim + " to " + pdim + " via PCA."); } proj = new double[pdim][dim]; - for(int d = 0; d < pdim; d++) { + for (int d = 0; d < pdim; d++) { EigenPair ep = axes.get(d); double[] ev = ep.getEigenvector().getArrayRef(); double eval = Math.sqrt(ep.getEigenvalue()); // Fill weighted and transposed: - for(int i = 0; i < dim; i++) { + for (int i = 0; i < dim; i++) { proj[d][i] = ev[i] / eval; } } @@ -156,7 +160,7 @@ public class GlobalPrincipalComponentAnalysisTransform<O extends NumberVector<?> @Override protected O filterSingleObject(O obj) { // Shift by mean and copy - for(int i = 0; i < dim; i++) { + for (int i = 0; i < dim; i++) { buf[i] = obj.doubleValue(i) - mean[i]; } double[] p = VMath.times(proj, buf); @@ -171,12 +175,12 @@ public class GlobalPrincipalComponentAnalysisTransform<O extends NumberVector<?> @Override protected SimpleTypeInformation<? super O> convertedType(SimpleTypeInformation<O> in) { initializeOutputType(in); - if(proj.length == dim) { - return in; - } - else { - return new VectorFieldTypeInformation<O>(factory, proj.length); - } + return new VectorFieldTypeInformation<>(factory, proj.length); + } + + @Override + protected Logging getLogger() { + return LOG; } /** @@ -201,15 +205,15 @@ public class GlobalPrincipalComponentAnalysisTransform<O extends NumberVector<?> protected void makeOptions(Parameterization config) { super.makeOptions(config); - ObjectParameter<EigenPairFilter> filterP = new ObjectParameter<EigenPairFilter>(FILTER_ID, EigenPairFilter.class, true); - if(config.grab(filterP)) { + ObjectParameter<EigenPairFilter> filterP = new ObjectParameter<>(FILTER_ID, EigenPairFilter.class, true); + if (config.grab(filterP)) { filter = filterP.instantiateClass(config); } } @Override - protected Object makeInstance() { - return new GlobalPrincipalComponentAnalysisTransform<O>(filter); + protected GlobalPrincipalComponentAnalysisTransform<O> makeInstance() { + return new GlobalPrincipalComponentAnalysisTransform<>(filter); } } -}
\ No newline at end of file +} |