It turns out that the sum of normally distributed and exponentially distributed random numbers gives rise to Carpenter’s distribution:
>> t=.2+.02*randn(1,10000) + exprnd(.05,[1 10000]);
>> hist(t,30)
>> hist(1./t,30)
>> t_reaccion2histog_carpenter(t,20,’r.-‘);
Leave a Reply