Dear R Development Core Team,
I would like to submit a bug in pbeta and, consequently, in pf.
R version 2.9.2 (2009-08-24)
Piotr Pokarowski, Ph. D.
Institute of Informatics,
University of Warsaw
Thank you for the report (I see you did report this a year ago already);
I had started looking at it back then, and am looking at it further now.
It needs a slight change inside the toms708 algo. -- Martin Mächler
resolved (for upcoming R 2.11.0) with svn revision 51327
I think I had a similar error with higher parameter values:
In pbeta(q, shape1, shape2, lower.tail, log.p) : production de NaN
and a discontinuous curve
R version 2.12.1 (2010-12-16)
Platform: x86_64-unknown-linux-gnu (64-bit)