https://researchdatapod.com/how-to-do-matrix-multiplication-in-r/