時系列のデータを折れ線グラフにし、グラフの平らな部分を検出したい
時系列のデータを折れ線グラフにするスクリプトは、下記のように書きました。
a <- read.table("Sequence2.csv",sep=",", header=T)
names(a) <- c("NO","Time","SequenceNO")
plot( a[,c("TIME")],a[,c("SequenceNO")],type="l")
グラフはこのようになったのですが、グラフの傾きが平らな部分を検出するにはどうすればいいでしょうか。
元のデータ:testfile.csv
TIME SequenceNO
1 0.000023 7
2 0.003774 1455
3 0.004244 2903
4 0.004645 4351
5 0.004667 5799
6 0.0047 7247
7 0.00472 8695
8 0.004818 13039
9 0.004839 18831
10 0.004936 20279
11 0.011538 27519
12 0.029395 34759
13 0.029433 49239
14 0.035552 60823
15 0.039416 60823
16 0.040157 60823
17 0.042644 60823
18 0.049411 60825
19 0.050914 60826
20 0.05136 76751
21 0.059166 101367
22 0.060155 120191
23 0.072393 141911
24 0.095405 152047
25 0.095636 191143
26 0.096089 211415
27 0.11215 262095
28 0.11589 276575