How to write better scientific code in Python?

Re:How to write better scientific code in Python? by ww_wpg on Saturday February 19, 2022 @04:14PM (#30400587)

> You cannot calculate the EV of a random variable X by taking the mean of random samples drawn from the distribution of X

My understanding of statistics is rudimentary so forgive me but doesn't the sample mean of a normally distributed variable tend towards the expected value for the population?

[ Parent | on HN ]

Re:How to write better scientific code in Python? by mturmon on Saturday February 19, 2022 @06:05PM (#30401384)

The sample mean of a Normal RV does tend towards the population expected value, sure.

The GP comment is talking about a Cauchy RV, which has heavy tails. So it has enough probability mass at large values that the expected value is infinite. Discarding constant scale factors, in this case:

  E[X] = Int{0..infty} x * p(x) dx 
       = Int{0..infty} x * (1/x^2) dx 
       = Int{0..infty} (1/x) dx 
       = +infty

So, the sample mean of Cauchy random variates will not converge to any real number.

[ Parent | on HN ]