web-dev-qa-db-ja.com

Rの時間差を計算する

Start.timeとend.timeを2つの変数として持つ、300万を超えるレコードを持つデータがあります。最初の10個のobsは次のとおりです。

   start.date start.time   end.date end.time
1  2012-07-13   15:01:32 2012-07-13 15:02:42
2  2012-07-05   18:26:31 2012-07-05 18:27:19
3  2012-07-14   20:23:21 2012-07-14 20:24:11
4  2012-07-29   16:09:54 2012-07-29 16:10:48
5  2012-07-21   14:58:32 2012-07-21 15:00:17
6  2012-07-04   15:36:31 2012-07-04 15:37:11
7  2012-07-22   18:28:31 2012-07-22 18:28:50
8  2012-07-09   21:08:42 2012-07-09 21:09:02
9  2012-07-05   09:44:52 2012-07-05 09:45:05
10 2012-07-02   18:50:47 2012-07-02 18:51:38

Start.timeとend.timeの差を計算する必要があります。

次のコードを使用しました。

mbehave11$diff.time <- difftime(mbehave11$end.time, mbehave11$start.time, units="secs")

しかし、私はこのエラーを受け取っています：

Error in as.POSIXlt.character(x, tz, ...) : 
  character string is not in a standard unambiguous format
In addition: Warning messages:
1: In is.na.POSIXlt(strptime(xx, f <- "%Y-%m-%d %H:%M:%OS", tz = tz)) :
  Reached total allocation of 1535Mb: see help(memory.size)
rdatetime
22
user1702490

must日付/時刻演算を行う前に、文字列を日付オブジェクトに変換します。これを試して：

a）データの読み取り：

R> dat <- read.table(textConnection("start.date start.time end.date end.time
2012-07-13   15:01:32 2012-07-13 15:02:42
2012-07-05   18:26:31 2012-07-05 18:27:19 
2012-07-14   20:23:21 2012-07-14 20:24:11"), header=TRUE)

b）1つの観察に取り組む：

 R>  strptime( paste(dat[,1], dat[,2]), "%Y-%m-%d %H:%M:%S")
 [1] "2012-07-13 15:01:32" "2012-07-05 18:26:31" "2012-07-14 20:23:21"

c）セットでの作業、数値への変換：

 R> as.numeric(difftime(strptime(paste(dat[,1],dat[,2]),"%Y-%m-%d %H:%M:%S"),
                        strptime(paste(dat[,3],dat[,4]),"%Y-%m-%d %H:%M:%S"))) 
 [1] -70 -48 -50
 R>
40
Dirk Eddelbuettel

Lubridateパッケージを使用できると思います
ymd_hmsというメソッドがあります
これを使用して文字列から時間を取得できます。大規模なデータセットの場合ははるかに高速です

library(lubridate)
dat <- read.table(textConnection("start.date start.time end.date end.time
2012-07-13   15:01:32 2012-07-13 15:02:42
2012-07-05   18:26:31 2012-07-05 18:27:19 
2012-07-14   20:23:21 2012-07-14 20:24:11"), header=TRUE)
starttime = ymd_hms(paste(dat[,1], dat[,2]))
endtime = ymd_hms(paste(dat[,3], dat[,4]))
interval = difftime(endtime,starttime,units = "secs")

または、1行で実行できますが、大きなデータセットの場合は時間がかかります。

difftime(paste(dat[,3], dat[,4]),paste(dat[,1], dat[,2]),units = "secs")
0
cloudscomputes