我一直在尝试这样做,但不知道该怎么做。 所以这是我可重复的数据:
library(ggplot2)
library(nlme)
mydata <- data.frame(SID=sample(1:150,400,replace=TRUE),age=sample(50:70,400,replace=TRUE), sex=sample(c("Male","Female"),200, replace=TRUE),time= sample(1:3,400, replace=TRUE), Vol =rnorm(400),HCD =rnorm(400))
mydata$time <- as.numeric(mydata$time)
model = lme(Vol ~ age+sex*time+time* HCD, random=~1|SID, na.action="na.omit",data=mydata);summary(model)
gp <- ggplot(mydata, aes(x=age, y=Vol, group=SID))+geom_line(aes(group=time, col=sex))
print(gp)
我试图将时间点限制为两点,但仍然看起来不太好。 关于如何改进这个的任何建议?
非常感谢!