1

我有一个这种格式的时间序列数据集:

       Time Val1 Val2
     0 0.68 0.39
    30 0.08 0.14
    35 0.12 0.07
    40 0.17 0.28
    45 0.35 0.31
    50 0.14 0.45
   100 1.01 1.31
   105 0.40 1.20
   110 2.02 0.57
   115 1.51 0.58
   130 1.32 2.01

使用这个数据集,我想提取(而不是预测)FC1=1 和 FC2=1 的时间。这是我用要提取的注释点创建的图。

在此处输入图像描述

我正在寻找一种使用或函数来插值/截取以提取值的解决方案。例如,如果我在倍数变化 1 处画一条直线(比如在 y 轴上),我想提取 X 轴上直线相交的所有点。

期待您的建议并提前致谢!

4

1 回答 1

2

您可以使用approxfun进行插值并uniroot找到单根(线交叉的地方)。您需要多次运行 uniroot 才能找到所有交叉点,该rle功能可能有助于选择起点。

数据中的 FC 值永远不会接近 1,更不用说交叉了,因此您必须拥有比显示的更多的数据,或者意味着不同的值。

如果您可以提供更多详细信息(可能包括显示您想要的图表),那么我们可能会提供更详细的帮助。

编辑

好的,这里有一些 R 代码可以找到线交叉的位置:

con <- textConnection('           Time Val1 Val2
         0 0.68 0.39
        30 0.08 0.14
        35 0.12 0.07
        40 0.17 0.28
        45 0.35 0.31
        50 0.14 0.45
       100 1.01 1.31
       105 0.40 1.20
       110 2.02 0.57
       115 1.51 0.58
       130 1.32 2.01')

mydat <- read.table(con, header=TRUE)

with(mydat, {
    plot( Time, Val1, ylim=range(Val1,Val2), col='green', type='l' )
    lines(Time, Val2, col='blue')
})
abline(h=1, col='red')

afun1 <- approxfun( mydat$Time, mydat$Val1 - 1 )
afun2 <- approxfun( mydat$Time, mydat$Val2 - 1 )
points1 <- cumsum( rle(sign(mydat$Val1 - 1))$lengths )
points2 <- cumsum( rle(sign(mydat$Val2 - 1))$lengths )

xval1 <- numeric( length(points1) - 1 )
xval2 <- numeric( length(points2) - 1 )

for( i in seq_along(xval1) ) {
    tmp <- uniroot(afun1, mydat$Time[ points1[c(i, i+1)] ])
    	xval1[i] <- tmp$root
}

for( i in seq_along(xval2) ) {
    tmp <- uniroot(afun2, mydat$Time[ points2[c(i, i+1)] ])
    	xval2[i] <- tmp$root
}

abline( v=xval1, col='green' )
abline( v=xval2, col='blue')
于 2013-03-08T00:42:52.200 回答