Здравствуйте,
помогите, пожалуйста, с выполнением задания. Я очень плох в R и не могу понять, что я делаю не так.
Задача такая:
"Изучите влияние ложного ограничения на оценку наименьших квадратов с помощью моделирования Монте-Карло (10000 итераций). Смоделируйте для каждой итерации 50 наблюдений для модели:
yᵢ = b₀ + b₁x₁ᵢ + b₂x₂ᵢ + uᵢ ,
где b₀ = -3, b₁ = 1, b₂ = 1. Ошибки uᵢ независимы и одинаково распределены как N(0,1) (стандартное нормальное распределение), x₁ᵢ распределён как N(-1, 10) для всех наблюдений и итераций, x₂ᵢ распределён как Bin(1, 0.4) для всех наблюдений и итераций, причём i = 1, ... , 50 (50 наблюдений в каждой итерации).
Необходимо оценить коэффициенты b₀, b₁, b₂, имея в виду следующие ограничения:
(1) без ограничений
(2) b₀ + b₁ + b₂ = -5/4
(3) b₁ + b₂ = 0
(4) b₁ = 1/2.
Сравните оценки методом наименьших квадратов, посмотрев на среднее арифметическое, эмпирическую дисперсию и среднеквадратичную ошибку".
Зачем нужно 10000 итераций и как их потом впихнуть в модель (или должно быть 10000 моделей, но зачем?), я не понял, но основная проблема в другом. Допустим, я смог построить модель на 50 наблюдениях:
b0 = -3
b1 = 1
b2 = 1
x1 = rnorm(50, mean = -1, sd = 10)
x2 = rbinom(50, 1, 0.4)
u = rnorm(50, 0, 1)
y = b0 + b1*x1 + b2*x2 + u
model50 = lm(y~x1+x2)
summary(model50),
вроде всё правильно, но не получается задать ограничения. Помогите, пожалуйста. Какую функцию надо использовать? Если пытаюсь использовать linearHypothesis, то всё время выдаёт ошибку, что я не так задаю ограничения. К тому же потом для сравнения нужно "вытащить" из результатов среднее, дисперсию и среднеквадратичную ошибку... В общем, помогите, пожалуйста, буду очень благодарен!
Задать ограничения для коэффициентов МНК в R
-
- Новоприбывший
- Сообщения: 1
- Зарегистрирован: 03 май 2021, 20:01
- Репутация: 0
- Откуда: Россия
-
- Гуру
- Сообщения: 4064
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1061
- Ваше звание: программист
- Откуда: Казань
-
- Гуру
- Сообщения: 4064
- Зарегистрирован: 15 окт 2010, 08:33
- Репутация: 1061
- Ваше звание: программист
- Откуда: Казань
Re: Задать ограничения для коэффициентов МНК в R
да, linearHypothesis - это про другое. Строим 100500 моделей, коеффициенты записываем в матрицу, протом по этой матрице считаем требуемые для зачета статистики
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 5 гостей