 2] - beta[k + 2]) / 6<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Nneighs[k] &lt;- 6<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;}<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;betamean[K - 1] &lt;- (2 * beta[K] + 4 * beta[K - 2] - beta[K - 3]) / 5<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Nneighs[K - 1] &lt;- 5<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;betamean[K] &lt;- 2 * beta[K - 1] - beta[K - 2] <br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;Nneighs[K] &lt;- 1<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;for (k in 1 : 