Ignore:
Timestamp:
Dec 26, 2006, 10:56:26 AM (15 years ago)
Author:
Jari Häkkinen
Message:

Addresses #170.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/yat/statistics/AveragerPair.cc

    r703 r718  
    4747  }
    4848
     49  double AveragerPair::ccc(void) const
     50  {
     51    return ( (x_.variance() && y_.variance() && (x_.mean()-y_.mean()) ) ?
     52             ((2*covariance()) /
     53              ((x_.variance()+y_.variance()) +
     54               (x_.mean()-y_.mean())*(x_.mean()-y_.mean()))) : 0);
     55  }
     56
     57  double AveragerPair::correlation(void) const
     58  { return ((x_.std()>0 && y_.std()>0) ?
     59            (covariance() / (x_.std()*y_.std()) ) : 0);
     60  }
     61
     62  double AveragerPair::covariance(void) const
     63  {
     64    return (n()>1) ? (xy_ - x_.sum_x()*y_.mean()) / n(): 0;
     65  }
     66
     67  double AveragerPair::mean_xy(void) const
     68  {
     69    return xy_/n();
     70  }
     71
     72  double AveragerPair::msd(void) const
     73  {
     74    return ( x_averager().sum_xx()+y_averager().sum_xx()-2*sum_xy() )/n();
     75  }
     76
     77  unsigned long AveragerPair::n(void) const
     78  {
     79    return x_.n();
     80  }
     81
    4982  void AveragerPair::reset(void)
    5083  {
     
    5891  }
    5992
     93  double AveragerPair::sum_xy(void) const
     94  {
     95    return xy_;
     96  }
     97
     98  double AveragerPair::sum_xy_centered(void) const
     99  {
     100    return xy_-x_.sum_x()*y_.mean();
     101  }
     102
     103  const Averager& AveragerPair::x_averager(void) const
     104  {
     105    return x_;
     106  }
     107
     108  const Averager& AveragerPair::y_averager(void) const
     109  {
     110    return y_;
     111  }
     112
    60113  const AveragerPair& AveragerPair::operator+=(const AveragerPair& a)
    61114  {
Note: See TracChangeset for help on using the changeset viewer.