
<span class="sd">    &gt;&gt;&gt; y, c, yhat = simulate_y_c_yhat(0.3, 0.2, 0.2, n=3, random_state=42)</span>
<span class="sd">    &gt;&gt;&gt; print(y, c, yhat)</span>
<span class="sd">    [ 0.36959213 -1.16147869  2.34335464] [-0.86093366 -0.04858751  0.80259507] [-0.02156406 -0.43175879  0.61910925]</span>

<span class="sd">    &quot;&quot;&quot;</span>
    <span class="n">rng</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">random</span><span class="o">.</