I am trying to do parallel computing using the apply family of functions with the following lines of code. The objective is to fit each column of the silum_ matrix to my specifi