Я создаю набор псевдоданных из 'mvrnorm', и я хотел бы создать новые наборы данных для каждого моделирования, чтобы сравнить мощность различных регрессионных моделей.В настоящее время я надеюсь сравнить три регрессионные модели: 1) линейная регрессия, 2) смешанная модель и 3) GEE.
Скажите, что мой набор данных называется "df".Для каждой итерации симуляции у меня есть
for(i in 1:nSim){
plM = summary(lm(rowMeans(df ~ pred))$coef[2,4]
pLMM = ...
pGEE = ...}
, так что мой цикл симуляции возвращает разные p-значения (в настоящее время я получаю одинаковые p-значения для каждой регрессионной модели).
Как мне создать новый набор данных в цикле моделирования?
Спасибо, ваша помощь будет высоко оценена.