R 随机模型公式对象

R 随机模型公式对象,r,formula,lme4,R,Formula,Lme4,我想把公式放在随机模型中,但我认为下面的错误是由于错误的公式对象(?),但无法修复它 set.seed(1234) mydata <- data.frame (A = rep(1:3, each = 20), B = rep(1:2, each = 30), C = rnorm(60, 10, 5)) mydata$A <- as.factor(mydata$A) mydata$B <- as.factor(mydata$B) myf

我想把公式放在随机模型中,但我认为下面的错误是由于错误的公式对象(?),但无法修复它

set.seed(1234)
mydata <- data.frame (A = rep(1:3, each = 20), B = rep(1:2, each = 30),
                      C = rnorm(60, 10, 5))
mydata$A <- as.factor(mydata$A)
mydata$B <- as.factor(mydata$B)
myfunction <- function (mydata, yvars, genovar, replication) { 
              require("lme4")
              formula = paste ("yvars" ~ 1|"genovar" + 1|"replication")
              model1 <- lmer(formula, data = dataframe, REML = TRUE)
              return(ranef(model2))
}

myfunction(mydata=dataf, yvars = "C", genovar = "A", replication = "B")
set.seed(1234)

我的数据这里有一些不可靠的东西,但我认为这是你想要的

set.seed(1234)
mydata <- data.frame (A = factor(rep(1:3, each = 20)),
                      B = factor(rep(1:2, each = 30)),
                      C = rnorm(60, 10, 5))
require("lme4")
myfunction <- function (mydata, yvars, genovar, replication) { 
  formula <- paste (yvars,"~ (1|",genovar,") + (1|",replication,")")
  model1 <- lmer(as.formula(formula), data = mydata, REML = TRUE)
  return(ranef(model1))
}
myfunction(mydata=mydata, yvars = "C", genovar = "A", replication = "B")
set.seed(1234)

我的数据这里有一些不可靠的东西,但我认为这是你想要的

set.seed(1234)
mydata <- data.frame (A = factor(rep(1:3, each = 20)),
                      B = factor(rep(1:2, each = 30)),
                      C = rnorm(60, 10, 5))
require("lme4")
myfunction <- function (mydata, yvars, genovar, replication) { 
  formula <- paste (yvars,"~ (1|",genovar,") + (1|",replication,")")
  model1 <- lmer(as.formula(formula), data = mydata, REML = TRUE)
  return(ranef(model1))
}
myfunction(mydata=mydata, yvars = "C", genovar = "A", replication = "B")
set.seed(1234)

mydata感谢Ben解决了关于复制大小的问题。提供的示例只是为了显示我得到的错误,但实际数据有大量的复制…再次感谢您的建议感谢Ben解决了这个问题以及关于复制大小的问题。提供的示例只是为了显示我得到的错误,但实际数据有大量复制…再次感谢您的建议