We're running into trouble with taking draws from the beta distribution on 2c.
Our function for 2b gives us nice looking curves, but when alpha & beta < 1,
the curve goes to infinity at the extremes. So when we take draws from the
beta distribution, how are we to weigh the probability of those infinite
extremes? rbeta() seems to give a result consistent with the shape of the
curves - i.e.
for (alpha & beta = really small) --> returns 0 and 1s
for (alpha & beta = really big) --> returns values very close to 0.5
Or, basically, is there a forumla for the probability mass function for a beta
distribution?
Thanks,
Phillip
-------------------------------------------------
Phillip Y. Lipscy
Perkins Hall Room #129
35 Oxford Street
Cambridge, MA 02138
(617)493-4893 DORM
(617)851-8220 CELL
lipscy(a)fas.harvard.edu
http://www.people.fas.harvard.edu/~lipscy/
First Year Student, Ph.D. Program
Harvard University, FAS, Department of Government
-------------------------------------------------