假设我们有这些时间序列的回归线,我想知道哪些更有可能首先达到某个值(例如4),或者达到某个值的时间概率分布。
date series1 series2
01-04 2 1.8
02-01 2.075 2.3
03-01 2.15 2.1
04-08 2.225 2.2
04-09 2.3 2.4
lr1=lm(series1~ date)
lr2=lm(series2~ date)
答案 0 :(得分:1)
由于date
列意味着我们将使用1:nrows为日期的问题没有解释。显示的结果是每个系列达到4的时间。
tt <- 1:nrow(DF)
co <- coef(lm(cbind(series1, series2) ~ tt, DF))
# solving 4 = a + b * t for t gives:
(4 - co[1,]) / co[2,]
## series1 series2
## 27.667 19.727
注意:我们将此用作输入DF
:
Lines <- "date series1 series2
01-04 2 1.8
02-01 2.075 2.3
03-01 2.15 2.1
04-08 2.225 2.2
04-09 2.3 2.4"
DF <- read.table(text = Lines, header = TRUE)