An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110519/b96bb695/attachment.pl>
dbetagen function in mc2d package
3 messages · krusty the klown, Ben Bolker
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
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110520/26219b0c/attachment.pl>