我有一些时间序列数据,我在ggplot2中拟合黄土曲线,如附图所示.数据采用"S"曲线的形状.我真正需要知道的是数据开始趋于平稳的日期,它看起来正好在'550'或'600'时间
是否有某种定量方式可以在图中标出？
指向数据集的链接:http://dl.dropbox.com/u/75403/stover_data.txt
A dput()数据集:
structure(list(date = c(211L, 213L, 215L, 217L, 218L, 221L, 222L,
223L, 224L, 225L, 226L, 228L, 229L, 230L, 231L, 232L, 233L, 234L,
235L, 236L, 237L, 238L, 239L, 240L, 241L, 242L, 244L, 246L, 247L,
248L, 249L, 250L, 251L, 253L, 254L, 255L, 256L, 258L, 259L, 260L,
261L, 262L, 263L, 264L, 265L, 266L, 267L, 268L, 269L, 270L, 271L,
272L, 273L, 274L, 275L, 276L, 277L, 278L, 279L, 281L, 282L, 283L,
285L, 286L, 287L, 288L, 290L, 291L, 292L, …我使用ggplot2制作了一张包含四条生长曲线的图表.
如果有人想尝试,希望下面的代码应该生成图表.
我想找到每条线上最大斜率的值,比如4个时间点.
任何人都可以提出任何想法如何去做？
library(ggplot2)
dat <- structure(list(TIME = c(0L, 2L, 4L, 6L, 8L, 10L, 12L, 14L, 16L,
18L, 20L, 22L, 24L, 26L, 28L, 30L, 0L, 2L, 4L, 6L, 8L, 10L, 12L,
14L, 16L, 18L, 20L, 22L, 24L, 26L, 28L, 30L, 0L, 2L, 4L, 6L,
8L, 10L, 12L, 14L, 16L, 18L, 20L, 22L, 24L, 26L, 28L, 30L, 0L,
2L, 4L, 6L, 8L, 10L, 12L, 14L, 16L, 18L, 20L, 22L, 24L, 26L,
28L, 30L), OD600 = c(0.2202, 0.2177, 0.2199, 0.2471, 0.2834, …