2

对于 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)

我有这个情节

情节 1

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

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

这就是我所拥有的：

情节2

我想要的是拟合一条曲线，就像我在图像 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()

但我没有从中得到任何东西。

有人可以帮我吗？

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

在此处输入图像描述

百万谢谢。

莉莉

4

0 回答 0