使用lmer()的预测矩阵
我正试图用大量的预测值来拟合一个模型,这样在模型公式中枚举它们就很乏味了。这对于lm()来说很简单: 除了我想包含大量随机效应项之外。而不是必须键入使用lmer()的预测矩阵,r,lme4,R,Lme4,我正试图用大量的预测值来拟合一个模型,这样在模型公式中枚举它们就很乏味了。这对于lm()来说很简单: 除了我想包含大量随机效应项之外。而不是必须键入(1 | A)。。。(1 | ZZZ),我希望以类似于上文中用于olsModel2的矩阵方法的方式包含每个预测值。显然,以下措施不起作用: meModel2 <- lmer(someDV ~ (1 | as.matrix(indicatorMatrix))) meModel2我认为将公式构造为字符串,然后使用as.formula,类似于 re
(1 | A)。。。(1 | ZZZ)
,我希望以类似于上文中用于olsModel2
的矩阵方法的方式包含每个预测值。显然,以下措施不起作用:
meModel2 <- lmer(someDV ~ (1 | as.matrix(indicatorMatrix)))
meModel2我认为将公式构造为字符串,然后使用as.formula
,类似于
restring1 <- paste0("(1 | ",colnames(indicatorMatrix),")",collapse="+")
form <- as.formula(paste0("someDV ~",restring1))
meModel1 <- lmer(form, data = data.frame(someDV,indicatorMatrix))
restring1您也可以使用reformate(restring1,response='someDV')
这必须是最快满意答案的新记录。现在我已经学会了一种简便的公式构建新方法!非常感谢你们两位。
meModel2 <- lmer(someDV ~ (1 | as.matrix(indicatorMatrix)))
restring1 <- paste0("(1 | ",colnames(indicatorMatrix),")",collapse="+")
form <- as.formula(paste0("someDV ~",restring1))
meModel1 <- lmer(form, data = data.frame(someDV,indicatorMatrix))