2

R スタジオの記述プロットの場合、スパゲッティ プロットに回帰曲線を当てはめたいと思います。私が使用したスパゲッティプロットを作成するには：

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