convolution code snippets in r