Re: [R] dbetagen function in mc2d package

From: Ben Bolker <bbolker_at_gmail.com>
Date: Thu, 19 May 2011 17:05:13 +0000

krusty the klown <alessandrobarbiero <at> libero.it> writes:

> I found this useful package for generalized beta, yet the function that
> calculates its density leaves me puzzled, especially when I plotted it:
>
> plot(function(y) dbetagen(y,4,1.2,min=0,max=40),xlim=c(0,40))
>
> The area between the density curve and the x axis does not seem to measure
> 1!!!

  It would help to tell us that this is in the "mc2d" package ...

> integrate(dbetagen,shape1=4,shape2=1.2,min=0,max=40,0,40)
40 with absolute error < 2e-06

  I would suggest that you report this as a bug to the package maintainers [maintainer("mc2d")]. Seems fairly easy to fix for yourself, though -- just scale by the size of the domain (max-min).

  Ben Bolker



R-help_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Thu 19 May 2011 - 17:08:16 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 20 May 2011 - 14:30:08 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive