Написание функции R, которая выводит оценки максимального правдоподобия?

Question

Как бы я написал функцию R, которая берет вектор отклика и ковариатную матрицу X и выводит вектор оценок максимального правдоподобия коэффициентов β, где μ = Xβ и E[Y] = µ, где Y - экспоненциально распределенный случайный вектор?

Я могу использовать qr() и связанные с ним функции, и мне нужно запрограммировать итерационно переоцененный алгоритм наименьших квадратов из первых принципов, используя QR-декомпозицию.