="op">]</span><span class="op">)</span>     <span class="co"># verify orthogonality</span></span></code></pre></div>
<pre><code><span><span class="co">##           [,1]</span></span>
<span><span class="co">## [1,] 7.276e-12</span></span></code></pre>
<p>Continue in the same way, subtracting the projections of
<code>X[,3]</code> on the previous columns, and so forth</p>
<div class="sourceCode" id="cb8"><pre class="downlit sourceCode r">
<code class="sourceCode R"><span><span class="va">Z</span><span class="o