<div dir="ltr">It looks like that there is an error in stats.fatiguelife.rvs<br><br>Kolmogorov test fails<br>>>> stats.kstest('fatiguelife','',args=(5,),N=1000)<br>(0.093216666807115545, array(2.5853082230575808e-008))<br>
<br>Mean of sample<br><br>>>> stats.fatiguelife.stats(5,moments='m')<br>array(13.5)<br>>>> np.mean(stats.fatiguelife.rvs(5,size=1000))<br>26.683858360164475<br>>>> np.mean(stats.fatiguelife.rvs(5,size=10000))<br>
26.841525716395847<br>>>> np.mean(stats.fatiguelife.rvs(5,size=100000))<br>
26.730694604009678<br>>>> np.mean(stats.fatiguelife.rvs(5,size=100000)/2)<br>13.469823793800416<br><br>>>> stats.fatiguelife.stats(3,moments='m')<br>array(5.5)<br>>>> np.mean(stats.fatiguelife.rvs(3,size=100000))<br>
10.922712537094393<br>>>> np.mean(stats.fatiguelife.rvs(3,size=100000)/2)<br>5.5340854278553246<br><br>Variance of sample<br><br>>>> stats.fatiguelife.stats(3,moments='v')<br>array(110.25)<br>>>> np.var(stats.fatiguelife.rvs(3,size=1000000))<br>
440.1793356094052<br>>>> np.var(stats.fatiguelife.rvs(3,size=1000000)/2)<br>110.445022957997<br>>>> np.var(stats.fatiguelife.rvs(3,size=10000000)/2)<br>110.03364894832275<br><br>>>> stats.fatiguelife.stats(5,moments='v')<br>
array(806.25)<br>>>> np.var(stats.fatiguelife.rvs(5,size=1000000))<br>3222.4271388000293<br>>>> np.var(stats.fatiguelife.rvs(5,size=1000000)/2)<br>809.29193071702855<br><br>theoretical mean and cdf  look correct, according to <a href="http://www.itl.nist.gov/div898/handbook/eda/section3/eda366a.htm">http://www.itl.nist.gov/div898/handbook/eda/section3/eda366a.htm</a><br>
but random number generator, is wrong by approximately the scale of 1/2<br><br>Josef<br></div>