我想用以下数据绘制密度时间序列:
密度矢量(4,2,5,8,4,6,4) 采样周期矢量(单位:秒)(2,2,2,2,3,2,2)
如您所见,采样周期不是恒定的。我只知道开始日期和时间。
我不知何故需要将开始时间日期分配给第一次测量,然后计算以下测量的以下日期和时间,但我不知道如何编码它。
答案 0 :(得分:1)
尝试首先转换ts中的所需向量,提供初始启动时间和句点的cumsum。 我假设你采样一个连续的过程(没有跨越/死亡时间)
require (lubridate)
require (tidyr)
require (ggplot2)
require (ggfortify)
require (timetk)
density <- c (4,2,5,8,4,6,4)
seconds <- c (2,2,2,2,3,2,2)
starttime <- 0
time <- 0 + cumsum (seconds)
df <- as.data.frame (cbind (time, seconds, density))
df$time <- as_datetime(df$time)
df$ts <- tk_ts (df, select = density)
autoplot (df$ts, ts.geom = 'bar', fill = 'blue')
答案 1 :(得分:0)
根据添加到开始的秒数的累积总和绘制密度。
dens <- c(4,2,5,8,4,6,4)
secs <- c(2,2,2,2,3,2,2)
st <- as.POSIXct("2000-01-01 00:00:00")
plot(st + cumsum(secs), dens, xlab = "", type = "l")