**使用可重复数据编辑**
我有一个data.frame,其中包含50种实验性治疗随时间的增长曲线。我将它们绘制为多面5x10图网格。考虑到我的实验治疗方法,我还以合理的方式订购了它们。
我运行了回归函数以查找每种处理的增长率,并将斜率值保存在另一个数据框中。我已经绘制了数据,回归线和增长率值,但是我想根据该回归斜率值为各个多面图的背景着色,但是我不知道如何设置颜色以调用一个连续变量,尤其是来自不同df且行数不同的连续变量(原始df有300行,我要调用的df有50行-每种处理一个)。
我的代码如下:
Df:
df <- data.frame(matrix(ncol = 3,nrow=300))
colnames(df) <- c("Trt", "Day", "Size")
df$Trt <- rep(1:50, each=6)
df$Day <- rep_len(1:6, length.out=300)
df$Size <- rep_len(c(3,5,8,9,12,12,3,7,10,16,17,20),length.out = 300)
回归函数和输出数据框:
regression=function(df){
reg_fun<-lm(formula=df$Size~df$Day)
slope<-round(coef(reg_fun)[2],3)
intercept<-round(coef(reg_fun)[1],3)
R2<-round(as.numeric(summary(reg_fun)[8]),3)
R2.Adj<-round(as.numeric(summary(reg_fun)[9]),3)
c(slope,intercept,R2,R2.Adj)
}
library(plyr)
slopevalues<-ddply(df,"Trt",regression)
colnames(slopevalues)<-c ("Trt","slope","intercept","R2","R2.Adj")
图:
ggplot(data=df, aes(x=Day, y=Size))+
geom_line() +
geom_point() +
xlab("Day") + ylab("Size (μm)")+
geom_smooth(method="lm",size=.5,se=FALSE)+
geom_text(data=slopevalues,
inherit.aes=FALSE,
aes(x =1, y = 16,hjust=0,
label=paste(slope)))+
facet_wrap(~ Trt, nrow=5)
我要做的是根据渐变上的斜率值(slopevalues $ slope)为各个图表的背景着色。我的真实数据不只是重复的2个值,所以我想根据该值在渐变的颜色上进行此操作。
欢迎任何建议。
答案 0 :(得分:0)
您可以将geom_rect
与无限坐标结合使用:
ggplot(data=df, aes(x=Day, y=Size))+
## This is the only new bit
geom_rect(
aes(fill = slope, xmin = -Inf, xmax = Inf, ymin = -Inf, ymax = Inf),
slopevalues,
inherit.aes = FALSE
) +
## New bit ends here
geom_line() +
geom_point() +
xlab("Day") + ylab("Size (μm)")+
geom_smooth(method="lm",size=.5,se=FALSE)+
geom_text(data=slopevalues,
inherit.aes=FALSE,
aes(x =1, y = 16,hjust=0,
label=paste(slope)))+
facet_wrap(~ Trt, nrow=5)