tree-tool/count.cpp at master · ncbi/tree-tool · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
// count.cpp

/*===========================================================================
*
*                            PUBLIC DOMAIN NOTICE
*               National Center for Biotechnology Information
*
*  This software/database is a "United States Government Work" under the
*  terms of the United States Copyright Act.  It was written as part of
*  the author's official duties as a United States Government employee and
*  thus cannot be copyrighted.  This software/database is freely available
*  to the public for use. The National Library of Medicine and the U.S.
*  Government have not placed any restriction on its use or reproduction.
*
*  Although all reasonable efforts have been taken to ensure the accuracy
*  and reliability of the software and data, the NLM and the U.S.
*  Government do not and cannot warrant the performance or results that
*  may be obtained by using this software or data. The NLM and the U.S.
*  Government disclaim all warranties, express or implied, including
*  warranties of performance, merchantability or fitness for any particular
*  purpose.
*
*  Please cite the author in any work or product based on this material.
*
* ===========================================================================
*
* Author: Vyacheslav Brover
*
* File Description:
*   Print statistics of a list of numbers
*
*/


#undef NDEBUG

#include "common.hpp"
#include "numeric.hpp"
using namespace Common_sp;
#include "version.inc"

#include "common.inc"


namespace
{


streamsize precision = 0;


void report (const string &attr,
             Real value)
{
  const ONumber on (cout, isInteger (value) ? 0 : precision, false);
  cout << attr << '\t' << value << '\n';
}


struct ThisApplication : Application
{
  ThisApplication ()
    : Application ("Statistics of a sequence of numbers from cin", false)
    {
      version = VERSION;
      addKey ("decimals", "Number of decimals", "2");
    }


	void body () const final
	{
	  precision = str2<streamsize> (getArg ("decimals"));


		MeanVar mv;
		size_t missing = 0;
		string s;
		while (cin >> s)
		  if (strNull (s))
		    missing++;
		  else
		  {
  			const Real x = str2real (s);
  			if (isNan (x))
  			  missing++;
  			else
  				mv << x;
  	  }

		report ("count",  mv. n);
		report ("mean",   mv. getMean ());
		report ("var",    mv. getVar ());
		report ("SD",     mv. getSD ());
	  report ("sum",    mv. s);
		report ("meanSD", mv. getSD () / sqrt (mv. n));
	  report ("min",    mv. v_min);
	  report ("max",    mv. v_max);

		report ("missing", (Real) missing);
	}
};


}  // namespace


int main (int argc,
          const char* argv[])
{
  ThisApplication app;
  return app. run (argc, argv);
}