On 8/27/2008 11:14 AM, [EMAIL PROTECTED] wrote:
hello,
gamma(a) works fine for a<173, but gamma(173) returns
[1] Inf
Warning message:
value out of range in 'gammafn'
What I really need is the ratio
gamma(nu/2) / gamma((nu+1)/2)
is anyone aware of a way to re-express this in a form that pushes the
threshold past which a warning message is generated?
Use logs:
exp( lgamma(nu/2) - lgamma((nu+1)/2) )
Duncan Murdoch
______________________________________________
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.
