У меня большой набор данных, и мне нужно сделать дизайн внутри предмета в пакете lmer с вложенной моделью.
Я полагаю, что это обозначение для дизайна предметов в R LM -
https://cran.r -project.org / web / packages / lme4 / vignettes / lmer.pdf страницы 5-7
lmer(DV ~ Condition:Country + (1 + Condition|subj) + (1 + Condition|Region), data = mydata)
На основании этой логики -
set.seed(5)
x <- rnorm(10)
x1<-x+rnorm(10)
x2<-x+1 + rnorm(10)
myDat <- data.frame(c(x1,x2), c(rep("x1", 10), rep("x2", 10)),
rep(paste("S", seq(1,10), sep=""), 2))
names(myDat) <- c("y", "x", "subj")
install.packages("lme4")
require(lme4)
mm1 <- lmer(y ~ x + (1|subj), data=myDat)
summary(mm1)
Я верю, что это сработает, теперь вот проблема. Мой DV находится в двух столбцах, это требует реструктуризации данных таким образом, что DV объединяется в 1 столбец и создается идентифицирующая переменная, указывающая условие. Другими словами, на данный момент у вас есть результат к DV1, результат к DV2. После реструктуризации у вас будет DV, а затем условие, указывающее, было ли это условием 1 или условием 2. Сложность этого заключается в том, что эта реструктуризация удвоит количество строк в наборе данных, что вызовет проблемы, когда мы захотим включить дополнительные столбцы в анализ, так как теперь в любом дополнительном числовом столбце будут дублированные значения.
Можно ли сгенерировать эту модель без реструктуризации набора данных?
Что-то вроде:
lmer(c(DV1,DV2) ~ Condition:Country + (1 + Condition|subj) + (1 + Condition|Region), data = mydata)