Pearson correlation between number of siblings and not-being-a-firstborn for a population where family size is uniformly distributed on the integers from 1 to n [...] ~0.57 for n=1

I'm confused: does this make sense for n=1? (Your code suggests that that should be n=2, maybe?)