<div><font face="comic sans ms,sans-serif">Dear scipy users, </font></div><div><font face="comic sans ms,sans-serif"></font> </div><div><font face="comic sans ms,sans-serif">The time for a diffusion process to reach a single evidence threshold a is often modeled as an inverse Gaussian distribution with mean (a/<span style>ó) and variance (<span style>a*ó<font><sup>2</sup>/ì<sup>3</sup></font></span>
</span>

), where  ì represents the mean drift rate and  ó<font><sup>2</sup>  the variance of the accumlulation process. How could I reparametrise the scipy.stats.<tt class="descname">invgauss</tt><em class="property"> </em> function to manipulate those parameters? </font></font></div>
<div><font face="Comic Sans MS"></font> </div><div><font face="Comic Sans MS">Cheers,</font></div><div><font face="Comic Sans MS">Mathieu</font></div><div><font size="3" face="Calibri"></font> </div><div><font size="3" face="Calibri"></font> </div>