title	author	date	output
多项式混合效应模型绘图	sjy	2023-05-29	html_document

```{r} #加载所用r包 library(openxlsx)#导入xlsx数据 library(tidyverse)#数据处理 library(lme4) #混合效应模型 library(lmerTest) #输出p值 library(modelr)#结果提取 library(ggplot2) #绘图 ``` # 1.导入数据及数据处理 ```{r} origindata % mutate(pc=as.factor(pc)) #将pc转为因子 origindata=origindata %>% rename(childid=pc) #改名 str(origindata) ``` # 2.混合效应模型建模 ```{r} model1 % add_predictions(model1) origindata2 % add_predictions(model2) origindata3 % add_predictions(model3) origindata4 % add_predictions(model4) ggplot(origindata1,aes(monthage,pred))+geom_line(aes(group=childid),color="red") ggplot(origindata2,aes(monthage,pred))+geom_line(aes(group=childid),color="red") ggplot(origindata3,aes(monthage,pred))+geom_line(aes(group=childid),color="red") ggplot(origindata4,aes(monthage,pred))+geom_smooth(aes(group=childid),color="red",linewidth=0.2) ```