R программирование - Линейная регрессия (два набора фиктивных переменных) с ограничениями равенства

Question

Я в настоящее время работает на линейной регрессии с двумя набором фиктивных переменных -

model1 < - лм (у ~ factor1 + factor2, данные = mydata)

где factor1 содержит 30 категорных переменных и factor2 содержит 5 категориальных переменных.

У меня есть два ограничения - сумма (коэффициент каждого factor1 * некоторые постоянный каждый член factor1) = 0 и сумму (коэффициент каждого factor2 * некоторые постоянный каждый член factor2) = 0.

я понимаю что решение каким-то образом связано с добавлением «контрастов» к функции lm, но я не знаю, как правильно его закодировать. Кто-нибудь, пожалуйста, просветите меня. Или, пожалуйста, не стесняйтесь предлагать альтернативные методы для работы с такими проблемами

Большое спасибо.

Answer 1

Можете ли вы использовать пакет glmc? Но я не знаю, как установить ограничения.