iag(1, 1 + N) --&gt;</span>
<span class="co"># &lt;!-- L1[cbind(1+(1:N), 1:N)] &lt;- -alpha --&gt;</span>
<span class="co"># &lt;!-- K1 &lt;- L1 %*% t(L1) --&gt;</span>
<span class="co"># &lt;!-- V1 &lt;- solve(K1) --&gt;</span>
<span class="co"># &lt;!-- sparsify(K1) --&gt;</span>
<span class="co"># &lt;!-- sparsify(V1) --&gt;</span>
<span class="co"># </span>
<span class="co"># &lt;!-- ## Dynamic linear models --&gt;</span>
<span class="co"># &lt;!-- N &lt;- 3 --&gt;</span>
<span class="co"># &lt;!-- L2 &