PCA: Quel est le problème avec cet algorithme?

Question

Quelqu'un peut-il vous plaît confirmer ou corriger cet algorithme de Wikipedia pour calculer le premier directeur composant? Je veux une implémentation simple de PCA dans D, qui n'a pas toutes les bibliothèques existantes pour PCA AFAIK. J'ai essayé de mettre en œuvre cela, et il ne semble pas que mes résultats sur des exemples simples correspondent à des choses que je reçois de R ou Octave. Je suis sûr que ce n'est pas un bug de mise en œuvre parce que je l'ai vérifié sur ma mise en œuvre à plusieurs reprises et même réécrit à partir de zéro, et il est un algorithme joli simple.

Answer 1

Votre matrice de données est zéro DIM Un, tel que requis par les conditions préalables?

L'extrait de code est clair sur ce point, mais la façon dont la matrice de données est à l'origine décrit, je pense que les moyens de colonne doivent être zéro pour la méthode de mise au travail.