Я бы хотел использовать функцию nlrob () в пакете robustbase в R для регрессии LAD.
Оценщик LAD является M-оценщиком со следующей функцией Psi:
psi = sign(x)
По этой причине я хотел бы использовать следующий код для запуска регрессии:
result = nlrob(y ~ a + x ^b, data = testdata,method = "M" ,psi = myfunction)
Мой вопрос: как я могу определить myfunction для получения оценок LAD?