Я знаю, что если я хочу смоделировать данные из нормального распределения, например, Y ~ N (x, y 2 ), я могу использовать
Y <- rnorm(n = 30,mean = x, sd = y)
и аналогично для других распределений.
Но возможно моделировать значения из рандомизированного полного уравнения проектирования блоков, например
Y xy = µ + T x + B y + E xy
где T x - лечебный эффект x, B y - блочный эффект y и E ху это остатки?