Loop function lapply for culculating thousands of correlation coefficients

Question

I am using a function called icc9, to calculate intraclass correlation coefficients for features over 3 time points. I would like to run the function icc9(y=response,subje