Question

我定期进行一系列测量，但有些测量遗失：

      measurement_date value
1  2011-01-17 13:00:00     5
2  2011-01-17 13:04:00     5
3  2011-01-17 13:08:00     7
4  2011-01-17 13:12:00     8
5  2011-01-17 13:16:00     4
6  2011-01-17 13:24:00     6
7  2011-01-17 13:28:00     5
8  2011-01-17 13:32:00     6
9  2011-01-17 13:36:00     9
10 2011-01-17 13:40:00     8
11 2011-01-17 13:44:00     6
12 2011-01-17 13:48:00     6
13 2011-01-17 13:52:00     4
14 2011-01-17 13:56:00     6

我有一个函数可以处理值并且可以处理缺失的值，但是行必须在那里，所以我生成一个每分钟有一行的数组，如下所示： / p>

times <- timeSequence(from=.., length=60, by="min")

现在我每小时都有一行，但我需要合并数据。我尝试过这样的事情，但不能说得对：

lapply(times, function(time) {
    n <- as.numeric(time)
    v <- Position(function(candidate) {
        y <- as.numeric(candiated)
        n == y
    }

    .. insert the value into the row here ..
}

但我只是得到错误和警告。我是否以正确的方式解决问题？我真的想要一个完整的＆＃34;具有每分钟值的数组，因为将有许多不同的函数将运行读数，如果它们可以假设它们全部存在，那么它就更容易实现它们。

Answer 1

DF <- data.frame(measurement_date = seq(as.POSIXct("2011-01-17 13:00:00"),
                                        as.POSIXct("2011-01-17 13:56:00"),
                                        by = "mins")[seq(1, 57, by = 4)][-6],
                 value = c(5,5,7,8,4,6,5,6,9,8,6,6,4,6))
full <- data.frame(measurement_date = seq(as.POSIXct("2011-01-17 13:00:00"),
                                          by = "mins", length = 60),
                   value = rep(NA, 60))

可以使用两种方法，第一种方法是merge：

> v1 <- merge(full, DF, by.x = 1, by.y = 1, all = TRUE)[, c(1,3)]
> names(v1)[2] <- "value" ## I only reset this to pass all.equal later
> head(v1)
     measurement_date value
1 2011-01-17 13:00:00     5
2 2011-01-17 13:01:00    NA
3 2011-01-17 13:02:00    NA
4 2011-01-17 13:03:00    NA
5 2011-01-17 13:04:00     5
6 2011-01-17 13:05:00    NA

第二个是通过使用%in%得出的指标变量：

> want <- full$measurement_date %in% DF$measurement_date
> full[want, "value"] <- DF[, "value"]
> head(full)
     measurement_date value
1 2011-01-17 13:00:00     5
2 2011-01-17 13:01:00    NA
3 2011-01-17 13:02:00    NA
4 2011-01-17 13:03:00    NA
5 2011-01-17 13:04:00     5
6 2011-01-17 13:05:00    NA
> all.equal(v1, full)
[1] TRUE

合并版本强烈首选，但需要一些工作。 %in%解决方案仅适用于此处，因为数据在DF和full中都按时间顺序排列，因此我之前的“首选”。然而，很容易按时间顺序获得/确保这两个对象，因此这两种方法都需要一些精细的工作。我们可以修改%in%方法以按顺序获取两个变量（重新开始full）：

full2 <- data.frame(measurement_date = seq(as.POSIXct("2011-01-17 13:00:00"),
                                           by = "mins", length = 60),
                    value = rep(NA, 60))
full2 <- full2[order(full2[,1]), ] ## get full2 in order
DF2 <- DF[order(DF[,1]), ]         ## get DF in order
want <- full$measurement_date %in% DF$measurement_date
full2[want, "value"] <- DF2[, "value"]

>     all.equal(full, full2)
[1] TRUE
>     all.equal(full2, v1)
[1] TRUE
>

Answer 2

在你的函数中，as.numeric（candiated）应该是as.numeric（候选者）。还有一个支架丢失。我不知道你在函数中想要实现的是什么，但它对我来说看起来非常复杂。

尝试

merge(Data,times,by.x=1,by.y=1,all.y=T)

这应该给你一些工作。

如何在R中组合两个不同长度的向量

2 个答案: