How to Use QR Factorization with qr.solve() Function in R for Linear Regression Lines