18

時系列の 1 つの列を含むデータ ファイルがあります。

'2012-02-01 17:42:44'
'2012-02-01 17:42:44'
'2012-02-01 17:42:44'

...時間の先頭にカウントがあるようにデータを分割したい。言う：

'2012-02-01 17:00:00'  20   
'2012-02-01 18:00:00'  30

「20」と「30」は、そのアウト期間の時系列エントリの数を表します。そして、時間とその「カウント」をグラフ化できるようにしたいと考えています。Rでこれを行うにはどうすればよいですか？

これが私の現在の折れ線グラフのプロットです。

library(ggplot2)

req <- read.table("times1.dat")
summary(req)

da <- req$V2
db <- req$V1

time <- as.POSIXct(db)

png('time_data_errs.png', width=800, height=600)
gg <- qplot(time, da) + geom_line()

print(gg)
dev.off()
4

1 に答える 1

31

cut1時間以内にいくつの値が発生するかを把握するために使用したいようです。

いくつかのサンプル データを提供していただけると、通常は役に立ちます。ここにいくつかあります：

set.seed(1) # So you can get the same numbers as I do
MyDates <- ISOdatetime(2012, 1, 1, 0, 0, 0, tz = "GMT") + sample(1:27000, 500)
head(MyDates)
# [1] "2012-01-01 01:59:29 GMT" "2012-01-01 02:47:27 GMT" "2012-01-01 04:17:46 GMT"
# [4] "2012-01-01 06:48:39 GMT" "2012-01-01 01:30:45 GMT" "2012-01-01 06:44:13 GMT"

tableand をcut(引数breaks="hour"(詳細については参照) とともに)使用?cut.Dateして、1 時間あたりの頻度を見つけることができます。

MyDatesTable <- table(cut(MyDates, breaks="hour"))
MyDatesTable
# 
# 2012-01-01 00:00:00 2012-01-01 01:00:00 2012-01-01 02:00:00 2012-01-01 03:00:00 
#                  59                  73                  74                  83 
# 2012-01-01 04:00:00 2012-01-01 05:00:00 2012-01-01 06:00:00 2012-01-01 07:00:00 
#                  52                  62                  64                  33 
# Or a data.frame if you prefer
data.frame(MyDatesTable)
#                  Var1 Freq
# 1 2012-01-01 00:00:00   59
# 2 2012-01-01 01:00:00   73
# 3 2012-01-01 02:00:00   74
# 4 2012-01-01 03:00:00   83
# 5 2012-01-01 04:00:00   52
# 6 2012-01-01 05:00:00   62
# 7 2012-01-01 06:00:00   64
# 8 2012-01-01 07:00:00   33

MyDatesTable最後に、オブジェクトのライン プロットを次に示します。

plot(MyDatesTable, type="l", xlab="Time", ylab="Freq")

ここに画像の説明を入力

cut時間間隔の範囲を処理できます。たとえば、30 分ごとに集計したい場合は、それbreaksを処理する引数を簡単に調整できます。

data.frame(table(cut(MyDates, breaks = "30 mins")))
#                   Var1 Freq
# 1  2012-01-01 00:00:00   22
# 2  2012-01-01 00:30:00   37
# 3  2012-01-01 01:00:00   38
# 4  2012-01-01 01:30:00   35
# 5  2012-01-01 02:00:00   32
# 6  2012-01-01 02:30:00   42
# 7  2012-01-01 03:00:00   39
# 8  2012-01-01 03:30:00   44
# 9  2012-01-01 04:00:00   25
# 10 2012-01-01 04:30:00   27
# 11 2012-01-01 05:00:00   33
# 12 2012-01-01 05:30:00   29
# 13 2012-01-01 06:00:00   29
# 14 2012-01-01 06:30:00   35
# 15 2012-01-01 07:00:00   33

アップデート

でプロットしようとしていたのでggplot2、ここに 1 つのアプローチがあります (必要なときにベース R のグラフィックスを通常使用するため、これが最適かどうかはわかりません)。

テーブルの を作成しdata.frame(上に示したように)、ダミーの「グループ」変数を追加して、次のようにプロットします。

MyDatesDF <- data.frame(MyDatesTable, grp = 1)
ggplot(MyDatesDF, aes(Var1, Freq)) + geom_line(aes(group = grp))

ここに画像の説明を入力

于 2012-11-30T17:30:56.800 に答える