statslib/logistic.cpp

Go to the documentation of this file.

 #include "statsLib.h"

dvariable plogis( const prevariable& x, const double& location, const double& scale )
{
  if( scale<=0 ) 
  {
    cerr<<"Standard deviation is less than or equal to zero in "
    "plogis( const dvariable& x, const double& location, const double& scale )\n";
    return 0;
  }
  return 1./(1.+mfexp((location-x)/scale));
}

dvariable plogis( const prevariable& x, const prevariable& location, const prevariable& scale )
{
  if( scale<=0 ) 
  {
    cerr<<"Standard deviation is less than or equal to zero in "
    "plogis( const dvariable& x, const dvariable& location, const dvariable& scale )\n";
    return 0;
  }
  return 1./(1.+mfexp((location-x)/scale));
}

dvar_vector plogis( const dvector& x, const prevariable& location, const prevariable& scale )
{
  if( scale<=0 ) 
  {
    cerr<<"Standard deviation is less than or equal to zero in "
    "plogis( const dvector& x, const dvariable& location, const dvariable& scale )\n";
    return 0;
  }
  return 1./(1.+mfexp((location-x)/scale));
}

dvector plogis( const dvector& x, const double& location, const double& scale )
{
  if( scale<=0 ) 
  {
    cerr<<"Standard deviation is less than or equal to zero in "
    "plogis( const dvector& x, const double& location, const double& scale )\n";
    return 0;
  }
  return 1./(1.+mfexp((location-x)/scale));
}

dvar_vector plogis( const dvar_vector& x, const prevariable& location, const prevariable& scale )
{
  if( scale<=0 ) 
  {
    cerr<<"Standard deviation is less than or equal to zero in "
    "plogis( const dvar_vector& x, const dvariable& location, const dvariable& scale )\n";
    return 0;
  }
  return 1./(1.+mfexp((location-x)/scale));
}

// template <typename T1, typename T2>
// T1 plogis( const T1& x, const T2& location, const T2& scale )
// {
//  if( scale<=0 ) 
//  {
//    cerr<<"Standard deviation is less than or equal to zero in "
//    "plogis( const dvar_vector& x, const dvariable& location, const dvariable& scale )\n";
//    return 0;
//  }
//  return 1./(1.+mfexp((location-x)/scale));
// }