STK_ILinearReduct.cpp

Go to the documentation of this file.

/*--------------------------------------------------------------------*/
/*     Copyright (C) 2004-2007  Serge Iovleff

    This program is free software; you can redistribute it and/or modify
    it under the terms of the GNU Lesser General Public License as
    published by the Free Software Foundation; either version 2 of the
    License, or (at your option) any later version.

    This program is distributed in the hope that it will be useful,
    but WITHOUT ANY WARRANTY; without even the implied warranty of
    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
    GNU Lesser General Public License for more details.

    You should have received a copy of the GNU Lesser General Public
    License along with this program; if not, write to the
    Free Software Foundation, Inc.,
    59 Temple Place,
    Suite 330,
    Boston, MA 02111-1307
    USA

    Contact : Serge.Iovleff@stkpp.org
*/

/*
 * Project:  stkpp::AAModels
 * created on: 17 avr. 2010
 * Purpose:  Abstract class for the computation of the Index in the
 * AAM models.
 * Author:   iovleff, serge.iovleff@stkpp.org
 **/

#include "../include/STK_ILinearReduct.h"
#include "../../Algebra/include/STK_LinAlgebra3D.h"

namespace STK
{
/*
 * Constructor.
 * @param data the input data set
 */
ILinearReduct::ILinearReduct( Matrix const& data) : IReduct(data)
{}

/*
 * Destructor
 */
ILinearReduct::~ILinearReduct()
{}

/* Compute the Index.
 *  @param nbAxis number of Axis to compute
 */
bool ILinearReduct::run()
{
  try
  {
    // maximize the Index and compute the axis
    maximizeIndex();
    // project data
    projection();

  } catch (Exception error)
  {
    msg_error_ = error.error();
    return false;
  }
  return true;
}

/*
 * Compute the weighted index.
 * @param weights the weights to used
 * @param nbAxis number of Axis to compute
 */
bool ILinearReduct::run( Vector const& weights)
{
  try
  {
    // maximize the Index and compute the axis
    maximizeIndex(weights);
    // project data
    projection();

  } catch (Exception error)
  {
    msg_error_ = error.error();
    return false;
  }
  return true;
}

void ILinearReduct::projection()
{
  // free allocated memory
  if (!p_reduced_) p_reduced_ = new Matrix;
  // compute matrix multiplication
  mult( *p_y_, axis_, *p_reduced_);
}

} // namespace STK