假设一个数据集包含单个时间步长(以秒为单位)和一个距离(以m为单位)每个时间步长管理的跑步者:

time <- c(300,600,900) # time in sec
dist <- c(1050,1950,3100) # dist in m

这里有一个情节来说明它:

plot(x=time,y=dist,xlim=c(0,1800),ylim=c(0,5000))

基于零截距(假设在时间0处没有距离),线性回归模型可以评估距离与时间的关系,即可以计算跑步者的速度 .

lm1 <- lm(dist~0+time)
abline(lm1,col="red")

此外,还可以轻松预测和绘制该特定线性模型的置信区间:

new_time <- seq(0,2000)
prd <- predict(lm1,newdata=data.frame(time=new_time),interval = c("confidence"),level = 0.90,type="response")
lines(new_time,prd[,"lwr"],col="red",lty=2)
lines(new_time,prd[,"upr"],col="red",lty=2)

Plot:Dist_Vs_Time

但是,当跑步者达到一定距离时,我感兴趣吗?例如,我想知道跑步者何时超过5000米 . 当然,可以从回归系数中提取/计算平均时间:

5000/(coef(lm1))[["time"]]

但是,我如何得到相应的置信区间,即图中置信区间的线在哪一时间穿过dist = 5000 m线?