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
|
/* BLURB lgpl
Coda File System
Release 5
Copyright (c) 1987-1999 Carnegie Mellon University
Additional copyrights listed below
This code is distributed "AS IS" without warranty of any kind under
the terms of the GNU Library General Public Licence Version 2, as
shown in the file LICENSE. The technical and financial contributors to
Coda are listed in the file CREDITS.
Additional copyrights
#*/
/*
IBM COPYRIGHT NOTICE
Copyright (C) 1986
International Business Machines Corporation
All Rights Reserved
This file contains some code identical to or derived from the 1986
version of the Andrew File System ("AFS"), which is owned by the IBM
Corporation. This code is provided "AS IS" and IBM does not warrant
that it is free of infringement of any intellectual rights of any
third party. IBM disclaims liability of any kind for any damages
whatsoever resulting directly or indirectly from use of this software
or of any derivative work. Carnegie Mellon University has obtained
permission to modify, distribute and sublicense this code, which is
based on Version 2 of AFS and does not contain the features and
enhancements that are part of Version 3 of AFS. Version 3 of AFS is
commercially available and supported by Transarc Corporation,
Pittsburgh, PA.
*/
#include <stdio.h>
#include <math.h>
double lowCI[20] = {-1.0, 3.078, 1.886, 1.638, 1.533, 1.476, 1.440,1.415, 1.397, 1.383, 1.372,
1.363, 1.356, 1.35, 1.345, 1.341, 1.337, 1.333, 1.33, 1.328};
/* Magic numbers used by CIFactor(); can't make it local */
double CIFactor(dFreedom)
int dFreedom; /* no of degrees of freedom; == (NoOfSamples - 1) */
{
/* Approximate; errs conservatively; see Law & Kelton, pg 386 */
if (dFreedom >= 100) return(1.29);
if (dFreedom >= 50) return(1.3);
if (dFreedom >= 30) return(1.31);
if (dFreedom >= 20) return (1.325);
return(lowCI[dFreedom]);
}
#if 0
void CumStats(count, sum, sum2)
int count;
double sum, sum2;
{
double mean, stddev, c90;
if (count < 2) return;
mean = sum/count;
stddev = sqrt((count*sum2 - sum*sum)/(count*(count-1)));
c90 = stddev*CIFactor(count - 1);
printf("SAMPLES = %d MEAN = %7.2f STD DEV = %7.2f 90%% CI = %7.2f\n",
count, mean, stddev, c90);
}
#endif
|
