r - 编写一个输出最大似然估计的 R 函数？

Question

我将如何编写一个 R 函数，它采用响应 vec y 和协变量矩阵 X 并输出系数 β 的最大似然估计向量，其中 μ = Xβ 和 E[Y] = μ，其中 Y 是指数分布的随机向量？

我可以使用 qr() 和相关函数，并且需要使用 QR 分解从第一原理对迭代重新加权最小二乘算法进行编程。