我想在以下data.frame df中为每个 id 拟合二次方(Time,SkinTemp)。每个ID都有不同数量的Time,SkinTemp条目,因此我坚持使用'预测'
df<-data.frame(Time=seq(65),
SkinTemp=rnorm(65,37,0.5),
id=rep(1:10,c(5,4,10,6,7,8,9,8,4,4)))
到目前为止,我有:
#Fit the model y=x^2+x+C
fitted_models = df %>% group_by(id) %>% do(model = lm(SkinTemp ~ Time+I(Time^2), data = .))
到目前为止一切顺利。这就是我被困的地方。如何将原始时间数据传递到下面的预测函数中?
#Predict data points for each quadratic
predQ<-sapply(unique(df$id), function(x) predict(fitted_models$model[[x]]))
答案 0 :(得分:2)
使用fitted
:
lapply(fitted_models$model, fitted)