<h3 dir="auto">Describe your issue.</h3>
<p dir="auto">Hi,</p>
<p dir="auto">Not sure what's going on here:</p>
<div class="highlight highlight-source-batchfile notranslate position-relative overflow-auto" dir="auto" data-snippet-clipboard-copy-content="&gt; python -c &quot;import scipy.stats; print('PASS')&quot;
PASS

&gt; python -c &quot;from scipy.stats import F_onewayResult; print('PASS')&quot;
Traceback (most recent call last):
  File &quot;&lt;string&gt;&quot;, line 1, in &lt;module&gt;
ImportError: cannot import name 'F_onewayResult' from 'scipy.stats' (C:\Users\runner\anaconda3\envs\...\lib\site-packages\scipy\stats\__init__.py)"><pre class="notranslate"><span class="pl-k">&gt;</span> python -c <span class="pl-s"><span class="pl-pds">"</span>import scipy.stats; print('PASS')<span class="pl-pds">"</span></span>
PASS

<span class="pl-k">&gt;</span> python -c <span class="pl-s"><span class="pl-pds">"</span>from scipy.stats import F_onewayResult; print('PASS')<span class="pl-pds">"</span></span>
Traceback (most recent <span class="pl-k">call</span> last):
  File <span class="pl-s"><span class="pl-pds">"</span>&lt;string&gt;<span class="pl-pds">"</span></span>, line 1, in <span class="pl-k">&lt;</span>module<span class="pl-k">&gt;</span>
ImportError: cannot import name 'F_onewayResult' from 'scipy.stats' (C:\Users\runner\anaconda3\envs\...\lib\site-packages\scipy\stats\__init__.py)</pre></div>
<p dir="auto">However, if I modify </p><div class="Box Box--condensed my-2">
  <div class="Box-header f6">
    <p class="mb-0 text-bold">
      <a href="https://github.com/scipy/scipy/blob/63472fd76eadf2c97b9aaf6f55ffc6ddce4b39b6/scipy/stats/__init__.py#L453">scipy/scipy/stats/__init__.py</a>
    </p>
    <p class="mb-0 color-fg-muted">
         Line 453
      in
      <a data-pjax="true" class="commit-tease-sha" href="/scipy/scipy/commit/63472fd76eadf2c97b9aaf6f55ffc6ddce4b39b6">63472fd</a>
    </p>
  </div>
  <div itemprop="text" class="Box-body p-0 blob-wrapper blob-wrapper-embedded data">
    <table class="highlight tab-size mb-0 js-file-line-container" data-tab-size="8" data-paste-markdown-skip="">

        <tbody><tr class="border-0">
          <td id="L453" class="blob-num border-0 px-3 py-0 color-bg-default" data-line-number="453"></td>
          <td id="LC453" class="blob-code border-0 px-3 py-0 color-bg-default blob-code-inner js-file-line"> <span class="pl-k">from</span> .<span class="pl-s1">_stats_py</span> <span class="pl-k">import</span> <span class="pl-c1">*</span> </td>
        </tr>
    </tbody></table>
  </div>
</div>
<br>
by adding the line:<p></p>
<div class="highlight highlight-source-python notranslate position-relative overflow-auto" dir="auto" data-snippet-clipboard-copy-content="from ._stats_py import F_onewayResult
from ._stats_py import *"><pre class="notranslate"><span class="pl-k">from</span> .<span class="pl-s1">_stats_py</span> <span class="pl-k">import</span> <span class="pl-v">F_onewayResult</span>
<span class="pl-k">from</span> .<span class="pl-s1">_stats_py</span> <span class="pl-k">import</span> <span class="pl-c1">*</span></pre></div>
<p dir="auto">I get a pass for both commands:</p>
<div class="highlight highlight-source-batchfile notranslate position-relative overflow-auto" dir="auto" data-snippet-clipboard-copy-content="&gt; python -c &quot;import scipy.stats; print('PASS')&quot;
PASS

&gt; python -c &quot;from scipy.stats import F_onewayResult; print('PASS')&quot;
PASS"><pre class="notranslate"><span class="pl-k">&gt;</span> python -c <span class="pl-s"><span class="pl-pds">"</span>import scipy.stats; print('PASS')<span class="pl-pds">"</span></span>
PASS

<span class="pl-k">&gt;</span> python -c <span class="pl-s"><span class="pl-pds">"</span>from scipy.stats import F_onewayResult; print('PASS')<span class="pl-pds">"</span></span>
PASS</pre></div>
<p dir="auto">Thanks in advance for your help!</p>
<h3 dir="auto">Reproducing Code Example</h3>
<div class="highlight highlight-source-python notranslate position-relative overflow-auto" dir="auto" data-snippet-clipboard-copy-content="from scipy.stats import F_onewayResult


print('PASS')"><pre class="notranslate"><span class="pl-k">from</span> <span class="pl-s1">scipy</span>.<span class="pl-s1">stats</span> <span class="pl-k">import</span> <span class="pl-v">F_onewayResult</span>


<span class="pl-en">print</span>(<span class="pl-s">'PASS'</span>)</pre></div>
<h3 dir="auto">Error message</h3>
<div class="highlight highlight-source-shell notranslate position-relative overflow-auto" dir="auto" data-snippet-clipboard-copy-content="Traceback (most recent call last):
  File &quot;&lt;string&gt;&quot;, line 1, in &lt;module&gt;
ImportError: cannot import name 'F_onewayResult' from 'scipy.stats' (C:\Users\runner\anaconda3\envs\...\lib\site-packages\scipy\stats\__init__.py)"><pre class="notranslate">Traceback (most recent call last):
  File <span class="pl-s"><span class="pl-pds">"</span>&lt;string&gt;<span class="pl-pds">"</span></span>, line 1, <span class="pl-k">in</span> <span class="pl-k">&lt;</span>module<span class="pl-k">&gt;</span>
ImportError: cannot import name <span class="pl-s"><span class="pl-pds">'</span>F_onewayResult<span class="pl-pds">'</span></span> from <span class="pl-s"><span class="pl-pds">'</span>scipy.stats<span class="pl-pds">'</span></span> (C:<span class="pl-cce">\U</span>sers<span class="pl-cce">\r</span>unner<span class="pl-cce">\a</span>naconda3<span class="pl-cce">\e</span>nvs<span class="pl-cce">\.</span>..<span class="pl-cce">\l</span>ib<span class="pl-cce">\s</span>ite-packages<span class="pl-cce">\s</span>cipy<span class="pl-cce">\s</span>tats<span class="pl-cce">\_</span>_init__.py)</pre></div>
<h3 dir="auto">SciPy/NumPy/Python version information</h3>
<p dir="auto">1.8.0 1.20.2 sys.version_info(major=3, minor=9, micro=10, releaselevel='final', serial=0)</p>
<p dir="auto">I have a fast Cython implementation of a Sobol low-discrepancy quasi-random number generator using Owen scrambling.  On top of that I have some transformations (inverse transform, Box-Muller) that make it easy to draw quasi-MC samples from, say, a multivariate normal distribution. This code has been used successfully in an industry setting.</p>
<p dir="auto">Before I go through the effort of creating a PR, I wanted to see whether this is something that people would like to see in scipy.</p>
<p dir="auto">The api could look something like the following:</p>
<div class="snippet-clipboard-content notranslate position-relative overflow-auto" data-snippet-clipboard-copy-content="from scipy.qmc import Sobol

sobol = Sobol(dimen=d, seed=1234)
samples = sobol.rvs(n)"><pre class="notranslate"><code class="notranslate">from scipy.qmc import Sobol

sobol = Sobol(dimen=d, seed=1234)
samples = sobol.rvs(n)
</code></pre></div>
<div class="snippet-clipboard-content notranslate position-relative overflow-auto" data-snippet-clipboard-copy-content="from scipy.qmc import MultivariateNormalQMC

mvn = MultivariateNormalQMC(mean=np.zeros(2), cov=np.array([[1, 0.1], [0.1, 1]]))
samples = mvn.rvs(n)"><pre class="notranslate"><code class="notranslate">from scipy.qmc import MultivariateNormalQMC

mvn = MultivariateNormalQMC(mean=np.zeros(2), cov=np.array([[1, 0.1], [0.1, 1]]))
samples = mvn.rvs(n)
</code></pre></div>