Question

我有一个与此类似（甚至更简单）的解决方案： How to interpolate data in R

但这对我不起作用。我已经使用dput拍摄了数据快照。

  a <- structure(list(
    time = structure(
       c(1566381574.097, 1566381542.104, 1566381510.109, 1566381390.134, 1566381330.118,
         1566381300.107, 1566381240.114, 1566381210.114, 1566381173.903, 1566381148.113
       ), class = c("POSIXct", "POSIXt"), tzone = ""),
    value = c(164.1, 162.3, 161.1, 158.6, 159.6, 157.6, 159, 157.8, 155.3, 155.3
  )), class = "data.frame", row.names = c(NA, 10L))

  b <- structure(list(time = structure(
    c(1566381120, 1566381180, 1566381240, 1566381300, 1566381360, 1566381420, 1566381480, 1566381540
    ), class = c("POSIXct", "POSIXt"), tzone = "")), row.names = c(NA, 8L), class = "data.frame")

  head(a)
  head(b)

  result <- b
  result$value <- approxfun(a$time,a$value)(b)

这就是我得到的：

>   head(a)
                 time value
1 2019-08-21 11:59:34 164.1
2 2019-08-21 11:59:02 162.3
3 2019-08-21 11:58:30 161.1
4 2019-08-21 11:56:30 158.6
5 2019-08-21 11:55:30 159.6
6 2019-08-21 11:55:00 157.6
>   head(b)
                 time
1 2019-08-21 11:52:00
2 2019-08-21 11:53:00
3 2019-08-21 11:54:00
4 2019-08-21 11:55:00
5 2019-08-21 11:56:00
6 2019-08-21 11:57:00
>   
>   result <- b
>   result$value <- approxfun(a$time,a$value)(b)
Error in .approxfun(x, y, v, method, yleft, yright, f) : 
  (list) object cannot be coerced to type 'double'

Answer 1

根据?approxfun

函数roxofun返回对给定数据点执行（线性或常数）插值的函数。对于给定的一组x值，此函数将返回相应的插值。   在这里，“ b”是一个data.frame，可能是我们需要的

x，y-数值向量，给出要插值的点的坐标。

基于解密，由于approxfun返回一个函数并且它基于一组x（numeric vector）进行插值，因此传递data.frame将不起作用。因此，提取特定的列作为插值向量

approxfun(a$time,a$value)(b$time)

R：roxofun抱怨“对象不能被强制键入'double'”

1 个答案: