Strange results for Beta Distribution

Dear All,
I got these results from the example in the function "dbeta":

x <- seq(0, 1, length=21)
     dbeta(x, 1, 1)

 [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1
Any Idea?

The density function for Beta(p,q) is (x^(p-1))*((1-x)^(q-1))/B(p,q).
In your example, p = q = 1, so the density function is 1 everywhere.

Best wishes,
Ted.


--------------------------------------------------------------------
E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk>
Fax-to-email: +44 (0)870 094 0861  [NB: New number!]
Date: 09-Nov-04                                       Time: 10:18:41
------------------------------ XFMail ------------------------------