diff options
| author | astrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675> | 2008-12-15 07:26:12 +0000 |
|---|---|---|
| committer | astrojhgu <astrojhgu@ed2142bd-67ad-457f-ba7c-d818d4011675> | 2008-12-15 07:26:12 +0000 |
| commit | 1f4a944064bc42284c33e6b755353d191cf288e8 (patch) | |
| tree | c8cb2253dea5f395e0f867aa6976433bd3eb00de /statistics/cstat.hpp | |
| download | opt-utilities-1f4a944064bc42284c33e6b755353d191cf288e8.tar.bz2 | |
git-svn-id: file:///home/svn/opt_utilities@1 ed2142bd-67ad-457f-ba7c-d818d4011675
Diffstat (limited to 'statistics/cstat.hpp')
| -rw-r--r-- | statistics/cstat.hpp | 66 |
1 files changed, 66 insertions, 0 deletions
diff --git a/statistics/cstat.hpp b/statistics/cstat.hpp new file mode 100644 index 0000000..b6b401d --- /dev/null +++ b/statistics/cstat.hpp @@ -0,0 +1,66 @@ +#ifndef CSTAT_HPP
+#define CSTAT_HPP
+#include <core/fitter.hpp>
+#include <iostream>
+
+using std::cout;using std::endl;
+namespace opt_utilities
+{
+ template<typename Ty,typename Tx,typename Tp,typename Ts,typename Tstr>
+ class cstat_poisson
+ :public statistic<Ty,Tx,Tp,Ts,Tstr>
+ {
+ private:
+ bool verb;
+ int n;
+
+ Ty lnfrac(Ty y)const
+ {
+ return y*log(y)-y;
+ }
+
+ public:
+ void verbose(bool v)
+ {
+ verb=v;
+ }
+ public:
+
+ statistic<Ty,Tx,Tp,Ts,Tstr>* do_clone()const
+ {
+ // return const_cast<statistic<Ty,Tx,Tp>*>(this);
+ return new cstat_poisson<Ty,Tx,Tp,Ts,Tstr>(*this);
+ }
+
+ Ts do_eval(const Tp& p)
+ {
+ Ts result(0);
+ for(int i=(this->datas()).size()-1;i>=0;--i)
+ {
+ Ty model_y=eval_model(this->datas().get_data(i).get_x(),p);
+ result+=model_y-this->datas().get_data(i).get_y()*log(model_y)+lnfrac(this->datas().get_data(i).get_y());
+ }
+
+ if(verb)
+ {
+ n++;
+ if(n%10==0)
+ {
+ cout<<result/this->p_fitter->get_dof()<<"\t";
+ for(int i=0;i<get_size(p);++i)
+ {
+ cout<<get_element(p,i)<<",";
+ }
+ cout<<endl;
+ }
+
+ }
+ return result;
+ }
+ };
+}
+
+#endif
+//EOF
+
+
|