我有一个机器故障数据,其中一列定义了故障间隔时间(tbf)
structure(list(tbf = c(2441, 2934, 4285, 2285, 4027, 2419, 2437, 2519, 3294, 2858, 3023, 2567, 3112, 2283, 3068, 2215, 3915, 2354.290323, 2477, 2258, 2742.5, 5198, 2837, 3282, 2474, 2883, 3837, 5054, 4874, 3559.5, 2783, 4246, 2602)), .Names = "tbf", class = "data.frame", row.names = c(NA, -33L))
我想绘制累积发生图。我可以用
做到这一点library(ggplot2)
ggplot(mydf, aes(x = tbf)) + stat_ecdf()
所示的图
但是,我希望直线适合这个情节。我不想要不均匀的线条,而是需要一根直线。我试过了
library(dplyr)
# add cumulative time and failures
mydf <- mydf %>% mutate(cumm_time = cumsum(tbf), cumm_fmode = row_number())
# fit linear regression
fit <- lm(cumm_time ~ cumm_fmode, data = mydf)
# plot points
plot(mydf$cumm_time, mydf$cumm_time)
# plot straight line
abline(fit)
我的要求是获得类似::
的数字我在哪里弄错了?任何帮助将不胜感激。