为R中的每个ID生成日期和时间序列

时间:2019-03-22 23:51:48

标签: r time time-series posixct

我正在尝试找出以这种格式创建日期和时间序列的方式:2018-01-01 01:00至2018-03-30 01:00 每个患者,然后用随机数填充新的空值。

我的数据如下:

Patients       temperature     
Patient1           37          
Patient2           36
Patient3           35.4

我想让数据看起来像

Patients       temperature       Time  
Patient1           37           2018-01-01 01:00
Patient2           36           2018-01-01 01:00
Patient3           35.4         2018-01-01 01:00
Patient1           NA           2018-01-01 02:00
Patient2           NA           2018-01-01 02:00
Patient3           NA           2018-01-01 02:00
Patient1           NA           2018-01-01 03:00
Patient2           NA           2018-01-01 03:00
Patient3           NA           2018-01-01 03:00

因此,时间变量将一直持续到2018年3月30日01:00,温度可以为NA,然后我生成随机数,但不会重复每位患者的相同温度值。

我尝试了此命令,但是没有用,我也不知道如何为每个患者分配时间

Time <- seq (from=as.POSIXct("2018-1-1 01:00"), to=as.POSIXct("2018-3-30 01:00", tz="UTC"), by="hour")

我也尝试过此命令,但收到错误消息:

dt = data.table(ID = Sensor7$StationID,Time = seq (from=as.POSIXct("2018-01-01 02:00"), to=as.POSIXct("2018-03-30 01:00",format = "%Y-%m-%d %H:%M",by="hour")))

但是它给了我错误信息:

Error in seq.POSIXt(from = as.POSIXct("2018-01-01 00:00"), to = as.POSIXct("2018-03-30 23:00",  : 
  exactly two of 'to', 'by' and 'length.out' / 'along.with' must be specified

有人知道如何以我正在寻找的请求形式获取数据吗?

2 个答案:

答案 0 :(得分:0)

这是一种方法:

dat = data.frame(Patients=paste0("Patients", 1:3), temperature=c(37,36,35.4))

Time = seq(as.POSIXct("2018-01-01 01:00"), as.POSIXct("2018-03-30 01:00"), by="hour")

new.data = data.frame(
  Patient = rep(dat$Patients, each=length(Time)),
  Time = rep(Time, length(dat$Patients))
)

我不确定您要如何生成随机值,但这是一个通用方法:

new.data$Random.Temperature = rnorm(nrow(new.data), 35, 1)

答案 1 :(得分:0)

离您不太远。试试这个:

cp

这应该做您想要的:

# I reproduce your data:

library(data.table)
data = data.table::fread(input = 
"Patients,temperature     
Patient1,37          
Patient2,36
Patient3,35.4")

library(dplyr)

Time <- seq (from=as.POSIXct("2018-1-1 01:00"), to=as.POSIXct("2018-3-30 01:00", tz="UTC"), by="hour")