一直使用的是visreg(data)多元回归分析直接作图,没有修改aes信息,下述内容实现两个目的:
①回归曲线是否显著?若不显著,将回归曲线标注为灰色
②回归曲线的R2标注在图上
首先实现判断回归曲线是否显著,使用fit_rs$coefficient[2,4],得到p值
然后对符合条件的p值图像进行颜色修改
if(p<0.05){
visreg(fit_rs,t_scale[j],line=list(col="grey"))
}else{
visreg(fit_rs,t_scale[j],line=list(col="blue"))
}
}
最后,对所有的多元回归图像进行R2标注
(暂时没有实现,或许可以通过分区进行?)
目前出现了一个新的问题,是图片的清晰度,远没有直接Export的高,需要调整图片的像素。
①使用png而不是jpeg,图像清晰度有一定的提升,但是仍然较模糊
②
③EPS格式的图可以用Adobe illustrator打开、编辑,在后续可以进行更复杂美观的图像编辑
这是与plot visreg相关的代码
# S3 method for visreg
plot(x, overlay=FALSE, print.cond=FALSE,
whitespace=0.2, partial=identical(x$meta$trans, I), band=TRUE,
rug=ifelse(partial, 0, 2), strip.names=is.numeric(x$fit[,x$meta$by]),
legend=TRUE, top=c('line', 'points'), gg=FALSE, line.par=NULL,
fill.par=NULL, points.par=NULL, ...)