Как смоделировать случайные числа Y из линейной модели с конкретными X и остатками?

Я хочу найти способ случайным образом сгенерировать 100 значений Y из линейной модели, где

Yi=  2−8Xi+ ei

Я хочу, чтобы остатки (ei) исходили из нормального распределения с указанным средним значением и дисперсией, а X был вектором значений от 1: 100.

Я знаю, как генерировать случайные переменные с помощью rnorm (), но я не уверен, как подойти к этому более сложному вопросу. Любые идеи о том, как я могу указать нужные мне параметры, будут приветствоваться.


r normal-distribution
person Community    schedule 29.11.2017    source источник

Ответы (1)


arrow_upward
2
arrow_downward

Это должно работать:

X <- 1:100
Y <-  2 − 8 * X + rnorm(100, mean = 0, sd = 2)

str(Y)
#num [1:100] -3.51 -12.03 -21.05 -31.38 -36.46 ...
person LyzandeR    schedule 29.11.2017
comment
Нет, остатки линейной модели всегда имеют нулевое среднее. - person Roland; 29.11.2017
comment
это отвечает на вопрос - person Andre Elrico; 29.11.2017
comment
На случай, если кому-то интересно, это математическое объяснение того, почему среднее всегда будет равно нулю. - person LyzandeR; 29.11.2017