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