1

我正在绘制残差与时间的关系,并在图上添加了一条黄土平滑曲线。当我将 x 比例更改为对数比例时,黄土平滑消失,之后我无法添加它。当我要求在 x 轴的对数转换后显示图形时,RStudio 中出现以下错误:

Error in seq.default(range[1], range[2], length = n) :'from' cannot be NA, NaN or infinite. 

我究竟做错了什么?如何在对数 x 尺度上添加黄土平滑?他是我的情节代码:

max.CWRES <- max(abs(data$CWRES),na.rm=T)
plotobj3 <- NULL
plotobj3 <-  ggplot(data[data$CONC!=0,])
plotobj3 <- plotobj3 + geom_point(aes(x=TIME, y=CWRES, colour=DOSE), shape=1, size=3)
plotobj3 <- plotobj3 + geom_abline(aes(x=TIME, y=CWRES),intercept=0, slope=0, colour="black")        

#Add loess smoothing line
plotobj3 <- plotobj3 + geom_smooth(aes(x=TIME, y=CWRES), method="loess", se=F, span=1.5,   colour="green")       
plotobj3 <- plotobj3+ scale_x_continuous(name="Time (hours)")
plotobj3 <- plotobj3+ scale_y_continuous(name="CWRES", limits=c(-max.CWRES ,max.CWRES))

# log scale for x-axis
plotobj3 <- plotobj3+ scale_x_log10(name="TIME (hours)")
plotobj3
4

1 回答 1

4

这里的问题似乎是您从TIME = 0. 当您log10-transform 您的x轴时,log-time 然后从-Inf. 您应该(至少对于 loes 部分)仅使用以下数据TIME > 0

plotobj3 + geom_smooth(aes(x = TIME, y=CWRES), 
                       data = data[data$TIME > 0 & data$CONC != 0, ], #this is the crucial part
                       method = "loess", 
                       se = FALSE, 
                       span = 1.5,   
                       colour = "green")        

请注意,这使用原始数据进行模型拟合。如果您想使用log10数据进行模型拟合,您可以aes(x=log10(TIME), y=CWRES)在任何地方使用,而不是aes(x=TIME, y=CWRES)

于 2014-11-07T08:57:21.120 回答