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.