/span><span class="op">)</span><span class="op">)</span>     <span class="co"># check orthogonality</span></span></code></pre></div>
<pre><code><span><span class="co">##       [,1] [,2] [,3] [,4]</span></span>
<span><span class="co">## [1,] 57888    0    0    0</span></span>
<span><span class="co">## [2,]     0 3249    0    0</span></span>
<span><span class="co">## [3,]     0    0    7    0</span></span>
<span><span class="co">## [4,]     0    0    0    2</span></span></code></pre>
<p>We make standardize co