-
Notifications
You must be signed in to change notification settings - Fork 59
/
asa189.html
151 lines (126 loc) · 3.34 KB
/
asa189.html
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
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
<html>
<head>
<title>
ASA189 - Maximum Likelihood Parameter Estimation, Beta Binomial PDF
</title>
</head>
<body bgcolor="#EEEEEE" link="#CC0000" alink="#FF3300" vlink="#000055">
<h1 align = "center">
ASA189 <br> Maximum Likelihood Parameter Estimation <br> Beta Binomial PDF
</h1>
<hr>
<p>
<b>ASA189</b>
is a FORTRAN90 library which
does a maximum likelihood estimation of the parameters of
a Beta Binomial PDF,
by D Smith.
</p>
<p>
<b>ASA189</b> is Applied Statistics Algorithm 189. Source code for many
Applied Statistics Algorithms is available through
<a href = "http://lib.stat.cmu.edu/apstat">STATLIB</a>.
</p>
<p>
This algorithm analyzes data that was generated by a beta binomial
probability density function, and tries to estimate the values of the
parameters that determing that PDF.
</p>
<h3 align = "center">
Languages:
</h3>
<p>
<b>ASA189</b> is available in
<a href = "../../f77_src/asa189/asa189.html">a FORTRAN77 version</a> and
<a href = "../../f_src/asa189/asa189.html">a FORTRAN90 version</a>.
</p>
<h3 align = "center">
Author:
</h3>
<p>
Original FORTRAN77 version by D Smith;
FORTRAN90 version by John Burkardt.
</p>
<h3 align = "center">
Reference:
</h3>
<p>
<ol>
<li>
D Smith,<br>
Algorithm AS 189:
Maximum Likelihood Estimation of the Parameters of the Beta
Binomial Distribution,<br>
Applied Statistics,<br>
Volume 32, Number 2, 1983, pages 196-204.
</li>
</ol>
</p>
<h3 align = "center">
Source Code:
</h3>
<p>
<ul>
<li>
<a href = "asa189.f90">asa189.f90</a>, the source code;
</li>
<li>
<a href = "asa189.sh">asa189.sh</a>,
commands to compile the source code;
</li>
</ul>
</p>
<h3 align = "center">
Examples and Tests:
</h3>
<p>
<ul>
<li>
<a href = "asa189_prb.f90">asa189_prb.f90</a>, a sample problem;
</li>
<li>
<a href = "asa189_prb.sh">asa189_prb.sh</a>,
commands to compile, link and run the sample problem;
</li>
<li>
<a href = "asa189_prb_output.txt">asa189_prb_output.txt</a>,
the output file.
</li>
</ul>
</p>
<h3 align = "center">
List of Routines:
</h3>
<p>
<ul>
<li>
<b>BBL</b> calculates the beta binomial log likelihood.
</li>
<li>
<b>BBME</b> estimates MU and THETA of the beta binomial distribution.
</li>
<li>
<b>BBML</b> estimates the parameters of a beta binomial distribution
</li>
<li>
<b>GDER</b> computes derivatives of the log likelihood.
</li>
<li>
<b>SET</b> sets up the arrays RL and LM.
</li>
<li>
<b>TIMESTAMP</b> prints the current YMDHMS date as a time stamp.
</li>
</ul>
</p>
<p>
You can go up one level to <a href = "../f_src.html">
the FORTRAN90 source codes</a>.
</p>
<hr>
<i>
Last revised on 09 January 2008.
</i>
<!-- John Burkardt -->
</body>
</html>