Re: [R] dbetagen function in mc2d package

From: krusty the klown <alessandrobarbiero_at_libero.it>
Date: Fri, 20 May 2011 03:19:31 -0700 (PDT)

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_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.
>
--
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_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 Fri 20 May 2011 - 14:28:02 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:40: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