我的研究是一个为期12周的研究，涉及两种饮食（饮食）。在时间（时间）= 0,6和12周的间隔进行结果测量，例如体重，腰围等。 （代码）代表主题。
我尝试下面的代码创建一个名为weight0的新变量，以仅在时间= 0时选择结果值。
weight0 <- dat2$weight[dat2$time==0]
我创建了以下模型。
p1 <- lme4::lmer(weight ~ diet * time + weight0 + weight0 * time + (time | code), REML = FALSE, data = dat2)
Error in model.frame.default(data = dat2, drop.unused.levels = TRUE, formula = weight ~ :
variable lengths differ (found for 'weight0')
我想我了解发生了什么事。权重0的值与结果权重的值长度不同。它是否正确？我该如何纠正这个问题，并使用基线作为固定效果？
感谢您抽出宝贵的时间阅读本文。
答案 0 :(得分：0)
您可以执行以下操作来定义data_female <- data[data$Gender == "Female",]
F <- cor(data_female[,-1], use = "pairwise.complete.obs")
F <- round(F, 1)
corrplot(F,
method="color",
col=col(200),
type="upper",
order="hclust",
addCoef.col = "black",
tl.col="black",
number.cex = 0.7,
tl.cex = 0.6,
tl.srt=45,
p.mat =p.mat,
sig.level = 0.5,
insig = "label_sig")
：
weight0然后拟合模型。