bbolker wrote:
>
> krusty the klown 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 ...
>
> * I wrote it on the title! :-)
>
>> 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).
>
> * That's ok, thanks!
>
> Ben Bolker
>
> ______________________________________________
> R-help@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.
>
--
View this message in context:
http://r.789695.n4.nabble.com/dbetagen-function-in-mc2d-package-tp3535411p3538128.html
Sent from the R help mailing list archive at Nabble.com.
[[alternative HTML version deleted]]
______________________________________________
R-help@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.
