在意大利面条图中适合曲线

Question

对于R studio中的描述性情节，我想在我的意大利面条情节中拟合回归曲线。要创建我使用的意大利面条图：

library(lattice)

GCIP <- data_head$GCIP
time_since_on <- data_head$time_since_on
Patient <- data_head$Patient
Eye <-data_head$Eye

xyplot(GCIP~time_since_on, groups = Patient, type='b', data=data_head)

我有这个情节

然后我想拟合多项式曲线，所以我使用了这段代码：

plot.new<- plot(time_since_on,GCIP)
lines(lowess(GCIP ~ time_since_on))

这就是我所拥有的：

我想要的是拟合曲线，就像我在图像2中得到的曲线，但是在意大利面条图上（每个主题的纵向数据）。

我试过使用这段代码：

library(ggplot2)
library(reshape2)

GCIP <- data_head$GCIP
time_since_on <- data_head$time_since_on
Patient.ID <- data_head$Patient.ID
Eye <-data_head$Eye
Visit <-data_head$Visit
Patient<-data_head$Patient

ggplot(data = reprex, aes(x,y)) +
  geom_point(alpha=1, size=2) +
  aes(colour=Patient.ID) +
  geom_text(aes(label=label), size=2, colour='white') +
  geom_path(aes(group=Patient.ID))

ggplot(data= reprex, aes(x = time_since_on, y = GCIP)) + 
  geom_point(size = 2, alpha= 1,  aes(color = Patient.ID)) + #colour points by group
  geom_path(aes(group = Patient.ID)) + #spaghetti plot
  stat_smooth(method = "lm", formula = y ~ x, aes(group = Patient.ID, colour = group)) + #line of best fit by group
  ylab("GCIP (volume)") + xlab("time_since_on (months)") +
  theme_bw()

但我从中得不到任何结果。

COuld有人帮我吗？

这里有一个来自互联网的例子

百万谢谢。

利利