Merge pull request #13 from jhamrick/master

Define X in confidence interval example
This commit is contained in:
Nicolas P. Rougier 2016-07-14 15:42:36 -05:00 committed by GitHub
commit bab502a047
2 changed files with 2 additions and 0 deletions

View File

@ -956,6 +956,7 @@ the rows which only contain integers and which sum to n. (★★★)</p>
<pre class="code python literal-block">
<span class="comment single"># Author: Jessica B. Hamrick</span>
<span class="name">X</span> <span class="operator">=</span> <span class="name">np</span><span class="operator">.</span><span class="name">random</span><span class="operator">.</span><span class="name">randn</span><span class="punctuation">(</span><span class="literal number integer">100</span><span class="punctuation">)</span> <span class="comment single"># random 1D array</span>
<span class="name">N</span> <span class="operator">=</span> <span class="literal number integer">1000</span> <span class="comment single"># number of bootstrap samples</span>
<span class="name">idx</span> <span class="operator">=</span> <span class="name">np</span><span class="operator">.</span><span class="name">random</span><span class="operator">.</span><span class="name">randint</span><span class="punctuation">(</span><span class="literal number integer">0</span><span class="punctuation">,</span> <span class="name">X</span><span class="operator">.</span><span class="name">size</span><span class="punctuation">,</span> <span class="punctuation">(</span><span class="name">N</span><span class="punctuation">,</span> <span class="name">X</span><span class="operator">.</span><span class="name">size</span><span class="punctuation">))</span>
<span class="name">means</span> <span class="operator">=</span> <span class="name">X</span><span class="punctuation">[</span><span class="name">idx</span><span class="punctuation">]</span><span class="operator">.</span><span class="name">mean</span><span class="punctuation">(</span><span class="name">axis</span><span class="operator">=</span><span class="literal number integer">1</span><span class="punctuation">)</span>

View File

@ -1113,6 +1113,7 @@ Thanks to Michiaki Ariga, there is now a
# Author: Jessica B. Hamrick
X = np.random.randn(100) # random 1D array
N = 1000 # number of bootstrap samples
idx = np.random.randint(0, X.size, (N, X.size))
means = X[idx].mean(axis=1)