Как вы выполняете ограниченные нелинейные наименьшие квадраты в R

Question

Я устанавливаю модель нелинейных наименьших квадратов в R. Я хочу свести к минимуму $ (Y - f (Xb))^2 $, где $ f $ является нелинейная монотонная дифференцируемая функция, $ X $ - набор признаков, а $ b $ - вектор параметра. Есть ли способ сделать это с ограничениями на $ b $? Я хочу ограничить $ b $ больше, чем 0, и я хочу, чтобы L1-стиль усадки некоторых элементов равнялся 0. Есть ли способ сделать это в R? nls() не допускает ограничений.

Answer 1

Вы можете преобразовать $ \ | \ boldsymbol {x} \ | _1 $ в свою цель, поставив ограничение на каждый элемент $ \ boldsymbol {x} $ в простую сумму и затем использовать quadprog для решения проблемы ,