You just want the following. Given that I have a p probability of winning a single trial, what's the chance that, after n trials, I win more than n*p times (you run good if you win more than your fair share, which is n*p).
This can be calculated with the cumulative binomial distribution, which is implemented in many languages. In R for instance we can define the following function:
Code:
probRunGood<-function(n,p) {
Nover<-n*p
pbinom(Nover,n,p,lower.tail=FALSE)
}
and see how it varies depending on the number of trials:
Code:
n<-1:1000
y<-probRunGood(n,0.6)
plot(n,y,cex=.2,ty="l")
If you run the plot, you'll see that the probability keeps oscillating around 0.5 while approaching to it as n grows.