summaryrefslogtreecommitdiff
path: root/src/de/lmu/ifi/dbs/elki/datasource/filter/transform/GlobalPrincipalComponentAnalysisTransform.java
diff options
context:
space:
mode:
authorErich Schubert <erich@debian.org>2013-10-29 20:02:37 +0100
committerAndrej Shadura <andrewsh@debian.org>2019-03-09 22:30:37 +0000
commitec7f409f6e795bbcc6f3c005687954e9475c600c (patch)
treefbf36c0ab791c556198b487ca40ae56ae5ab1ee5 /src/de/lmu/ifi/dbs/elki/datasource/filter/transform/GlobalPrincipalComponentAnalysisTransform.java
parent974d4cf6d54cadc06258039f2cd0515cc34aeac6 (diff)
parent8300861dc4c62c5567a4e654976072f854217544 (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.java50
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
+}