
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]]

