chart2/source/tools/MeanValueRegressionCurveCalculator.cxx

   1 /*************************************************************************
   2  *
   3  * DO NOT ALTER OR REMOVE COPYRIGHT NOTICES OR THIS FILE HEADER.
   4  *
   5  * Copyright 2008 by Sun Microsystems, Inc.
   6  *
   7  * OpenOffice.org - a multi-platform office productivity suite
   8  *
   9  * $RCSfile: MeanValueRegressionCurveCalculator.cxx,v $
  10  * $Revision: 1.8 $
  11  *
  12  * This file is part of OpenOffice.org.
  13  *
  14  * OpenOffice.org is free software: you can redistribute it and/or modify
  15  * it under the terms of the GNU Lesser General Public License version 3
  16  * only, as published by the Free Software Foundation.
  17  *
  18  * OpenOffice.org is distributed in the hope that it will be useful,
  19  * but WITHOUT ANY WARRANTY; without even the implied warranty of
  20  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  21  * GNU Lesser General Public License version 3 for more details
  22  * (a copy is included in the LICENSE file that accompanied this code).
  23  *
  24  * You should have received a copy of the GNU Lesser General Public License
  25  * version 3 along with OpenOffice.org.  If not, see
  26  * <http://www.openoffice.org/license.html>
  27  * for a copy of the LGPLv3 License.
  28  *
  29  ************************************************************************/
  30
  31 // MARKER(update_precomp.py): autogen include statement, do not remove
  32 #include "precompiled_chart2.hxx"
  33 #include "MeanValueRegressionCurveCalculator.hxx"
  34 #include "macros.hxx"
  35
  36 #include <rtl/math.hxx>
  37 #include <rtl/ustrbuf.hxx>
  38
  39 using namespace ::com::sun::star;
  40
  41 using ::rtl::OUString;
  42 using ::rtl::OUStringBuffer;
  43
  44 namespace chart
  45 {
  46
  47 MeanValueRegressionCurveCalculator::MeanValueRegressionCurveCalculator() :
  48         m_fMeanValue( 0.0 )
  49 {
  50     ::rtl::math::setNan( & m_fMeanValue );
  51 }
  52
  53 MeanValueRegressionCurveCalculator::~MeanValueRegressionCurveCalculator()
  54 {}
  55
  56 // ____ XRegressionCurveCalculator ____
  57 void SAL_CALL MeanValueRegressionCurveCalculator::recalculateRegression(
  58     const uno::Sequence< double >& /*aXValues*/,
  59     const uno::Sequence< double >& aYValues )
  60     throw (uno::RuntimeException)
  61 {
  62     const sal_Int32 nDataLength = aYValues.getLength();
  63     sal_Int32 nMax = nDataLength;
  64     double fSumY = 0.0;
  65     const double * pY = aYValues.getConstArray();
  66
  67     for( sal_Int32 i = 0; i < nDataLength; ++i )
  68     {
  69         if( ::rtl::math::isNan( pY[i] ) ||
  70             ::rtl::math::isInf( pY[i] ))
  71             --nMax;
  72         else
  73             fSumY += pY[i];
  74     }
  75
  76     m_fCorrelationCoeffitient = 0.0;
  77
  78     if( nMax == 0 )
  79     {
  80         ::rtl::math::setNan( & m_fMeanValue );
  81     }
  82     else
  83     {
  84         m_fMeanValue = fSumY / static_cast< double >( nMax );
  85
  86         // correlation coefficient: standard deviation
  87         if( nMax > 1 )
  88         {
  89             double fErrorSum = 0.0;
  90             for( sal_Int32 i = 0; i < nDataLength; ++i )
  91             {
  92                 if( !::rtl::math::isNan( pY[i] ) &&
  93                     !::rtl::math::isInf( pY[i] ))
  94                 {
  95                     double v = m_fMeanValue - pY[i];
  96                     fErrorSum += (v*v);
  97                 }
  98             }
  99             OSL_ASSERT( fErrorSum >= 0.0 );
 100             m_fCorrelationCoeffitient = sqrt( fErrorSum / (nMax - 1 ));
 101         }
 102     }
 103 }
 104
 105 double SAL_CALL MeanValueRegressionCurveCalculator::getCurveValue( double /*x*/ )
 106     throw (lang::IllegalArgumentException,
 107            uno::RuntimeException)
 108 {
 109     return m_fMeanValue;
 110 }
 111
 112
 113 uno::Sequence< geometry::RealPoint2D > SAL_CALL MeanValueRegressionCurveCalculator::getCurveValues(
 114     double min, double max, ::sal_Int32 nPointCount,
 115     const uno::Reference< chart2::XScaling >& xScalingX,
 116     const uno::Reference< chart2::XScaling >& xScalingY,
 117     ::sal_Bool bMaySkipPointsInCalculation )
 118     throw (lang::IllegalArgumentException,
 119            uno::RuntimeException)
 120 {
 121     if( bMaySkipPointsInCalculation )
 122     {
 123         // optimize result
 124         uno::Sequence< geometry::RealPoint2D > aResult( 2 );
 125         aResult[0].X = min;
 126         aResult[0].Y = m_fMeanValue;
 127         aResult[1].X = max;
 128         aResult[1].Y = m_fMeanValue;
 129
 130         return aResult;
 131     }
 132     return RegressionCurveCalculator::getCurveValues( min, max, nPointCount, xScalingX, xScalingY, bMaySkipPointsInCalculation );
 133 }
 134
 135 OUString MeanValueRegressionCurveCalculator::ImplGetRepresentation(
 136     const uno::Reference< util::XNumberFormatter >& xNumFormatter,
 137     ::sal_Int32 nNumberFormatKey ) const
 138 {
 139     OUStringBuffer aBuf( C2U( "f(x) = " ));
 140
 141     aBuf.append( getFormattedString( xNumFormatter, nNumberFormatKey, m_fMeanValue ));
 142
 143     return aBuf.makeStringAndClear();
 144 }
 145
 146 } //  namespace chart