Changeset 78


Ignore:
Timestamp:
May 4, 2004, 1:26:48 PM (19 years ago)
Author:
Peter
Message:

modified get_p_approx

Location:
trunk/src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/src/ROC.cc

    r75 r78  
    6262  {
    6363    double x = area - 0.5;
    64     double sigma = (pow((value_.size()-nof_pos_)* nof_pos_ *
    65                         (value_.size()+1)/12, 0.5) /
     64    // Not integrating from the middle of the bin, but from the inner edge.
     65    if (x>0)
     66      x -= 0.5/nof_pos_/(value_.size()-nof_pos_);
     67    else if(x<0)
     68      x += 0.5/nof_pos_/(value_.size()-nof_pos_);
     69
     70    double sigma = (std::sqrt((value_.size()-nof_pos_)* nof_pos_ *
     71                        (value_.size()+1)/12) /
    6672                    (value_.size() - nof_pos_ ) / nof_pos_);
    6773    double p = gsl_cdf_gaussian_Q(x, sigma);
  • trunk/src/ROC.h

    r77 r78  
    6161    ///
    6262    ///
    63     /// Implemented as in MatLab 12.1
     63    /// Implemented as in MatLab 13.1
    6464    /// @return the p-value
    6565    ///
Note: See TracChangeset for help on using the changeset viewer.