I am currently working a lot with R and Python. I am not able to access the C code the the R function lm_fit. I am wondering how is the linear regression optimize in R and python ?
I am pretty sure in R it is optimize with the Normal equation (since R is good with matrix multiplication) and I would think it is optimize with Gradient descent in Python ?
Am I right ?