Changeset 404 for library/bdm/stat/exp_family.h

Timestamp:

07/02/09 22:16:05 (15 years ago)

Author:

smidl

Message:

Change in epdf: evallog returns -inf for points out of support. Merger is aware of it now.

Files:

• library/bdm/stat/exp_family.h (modified) (4 diffs)

library/bdm/stat/exp_family.h

@@ -51 +51 @@
                                 double tmp;
                                 tmp= evallog_nn ( val )-lognc();
-                                it_assert_debug ( std::isfinite ( tmp ),"Infinite value" ); 
+//                              it_assert_debug ( std::isfinite ( tmp ),"Infinite value" ); 
                                 return tmp;}
                         //!Evaluate normalized log-probability for many samples
@@ -280 +280 @@
                         double evallog_nn ( const vec &val ) const
                         {
-                                double tmp; tmp= ( beta-1 ) *log ( val );               it_assert_debug ( std::isfinite ( tmp ),"Infinite value" );
+                                double tmp; tmp= ( beta-1 ) *log ( val );
+//                              it_assert_debug ( std::isfinite ( tmp ),"Infinite value" );
                                 return tmp;
                         };
@@ -290 +291 @@
                                 for ( int i=0;i<beta.length();i++ ) {lgb+=lgamma ( beta ( i ) );}
                                 tmp= lgb-lgamma ( gam );
-                                it_assert_debug ( std::isfinite ( tmp ),"Infinite value" );
+//                              it_assert_debug ( std::isfinite ( tmp ),"Infinite value" );
                                 return tmp;
                         };
@@ -489 +490 @@
 
                         double eval ( const vec &val ) const  {return nk;}
-                        double evallog ( const vec &val ) const  {return lnk;}
+                        double evallog ( const vec &val ) const  {
+                                if (any(val<low) && any(val>high)) {return inf;}
+                                else return lnk;
+                        }
                         vec sample() const
                         {