[R] birthday problem (factorial limit)

[Uwe Ligges]

Jörg Groß wrote:
> Hi,
> 
> I tried to calculate the formula for the birthday problem
> (the probability that at least two people out of a group of n people 
> share the same birthday)
> 
> But the factorial-function allows me only to calculate factorials up to 
> 170.
> 
> 
> So is there a way to push that limit?
> 
> to solve this formula:
> 
> (factorial(365) / factorial((365-23))) / (365^23)

Obviously you can easily rewrite this formula to:

prod(343:365) / (365^23)

or

factorial(23) * choose(365, 23) / (365^23)

Uwe Ligges


> (n=23)
> 
> ______________________________________________
> 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.