r - 最尤推定値を出力する R 関数を作成しますか?

Question

応答 vec y と共変量行列 X を取り、係数 β の最尤推定値のベクトルを出力する R 関数をどのように記述すればよいでしょうか。ここで、μ = Xβ および E[Y] = μ、ここで、Y は指数分布のランダム ベクトルです。

qr() および関連する関数を使用できますが、QR 分解を使用して、第一原理から繰り返し再重み付けされた最小二乗アルゴリズムをプログラムする必要があります。