我有两个data.frames,一个包含原始数据,另一个包含我从原始数据派生的建模系数.
更多细节:第一个data.frame"raw"包含"Time"(0s到900s)和"OD",用于许多Variants和4次运行.第二个data.frame"coef"每个Variant/run组合包含一行,该行中有各个系数("M","D.1"和"t0.1").
我已经绘制了每个Variant的原始数据分割并用runID着色,没问题.但现在我想根据runID覆盖模型曲线.
由于建模系数是在不同的data.frame(具有不同的维度,所以我不能只是cbind),stat_function对我不起作用.我可以一次显示曲线.
我尝试过for(){}循环,每次都添加一个stat_function图层:
p <- ggplot(temp, aes(Time, OD)) + geom_point(aes(colour = runID), size = 2) #works fine!
calc <- function(x){temp.n$M[ID] * (1 - exp(temp.n$D.1[ID] * temp.n$t0.1[ID] - x)))}
for(ID in 1:length(unique(temp.n$runID))) {
p <- p + stat_function(fun = calc)
}
print(p)
最后,所有"p"返回是原始数据和循环位的最终曲线的图.每当我尝试添加新的"stat_function"图层时,"p"似乎恢复到原始状态.
有任何想法吗？