我正在运行一个混合模型来预测特定时间点的连续值（以及其他变量值，但这只是一个简单的示例）。假设我的模型只有一个截距和时间，而时间是一个随机效应。我仅显示模型以显示从中获得固定效果的位置，下面将其放置。那是我唯一需要的部分。
library(lme4)
library(lmerTest)
mixed1<-lmer(donor_egfr~test_yrs+(test_yrs|DONOR_ID),data=donor_cr)
来自fixef（mixed1）的模型Beta：
fixef<-c(65,0.2)
我想计算时间0、1、2，... 15的估算值。所以我建立了一个数据框：
pred<-data.frame('time'=c(0:15))
如果要计算在时间15的估算值，我将使用：
contest1D(mixed1,c(1,15))
但是我想在所有时间点计算这些值。我不能使用这个：
p<-contest1D(mixed1,c(1,pred$time))
由于长度不匹配。我可以使用：
p<-contest1D(mixed1,c(1,pred$time[16]))
获得与第一条语句相同的结果。
所以本质上，我想知道如何将此race1D函数分别应用于pred $ time的每个值。我应该使用某种函数还是套用？
谢谢