3

Rセッションの例でおそらく最もよく説明されている、かなり奇妙な問題があります。以下のコメントで求められているように、これを再現可能にしようとしました。

meto <- structure(c(30, 25, 25, 25, 20, 20, 20, 20, 20, 20), index = structure(c(12796, 
            12796.0416666667, 12796.0833333333, 12796.125, 12796.1666666667, 
            12796.2083333333, 12796.25, 12796.2916666667, 12796.3333333333, 
            12796.375), format = structure(c("d/m/y", "h:m:s"), .Names = c("dates", 
            "times")), origin = structure(c(1, 1, 1970), .Names = c("month", 
            "day", "year")), class = c("chron", "dates", "times")), class = "zoo")

サンプル データセットは次のようになります。

> meto
(13/01/05 00:00:00) (13/01/05 01:00:00) (13/01/05 02:00:00) (13/01/05 03:00:00) (13/01/05 04:00:00) 
                 30                  25                  25                  25                  20 
(13/01/05 05:00:00) (13/01/05 06:00:00) (13/01/05 07:00:00) (13/01/05 08:00:00) (13/01/05 09:00:00) 
                 20                  20                  20                  20                  20 
> str(meto)
‘zoo’ series from (13/01/05 00:00:00) to (13/01/05 09:00:00)
  Data: num [1:10] 30 25 25 25 20 20 20 20 20 20
  Index: Classes 'chron', 'dates', 'times'  atomic [1:10] 12796 12796 12796 12796 12796 ...
  ..- attr(*, "format")= Named chr [1:2] "d/m/y" "h:m:s"
  .. ..- attr(*, "names")= chr [1:2] "dates" "times"
  ..- attr(*, "origin")= Named num [1:3] 1 1 1970
  .. ..- attr(*, "names")= chr [1:3] "month" "day" "year"

XTS に変換すると、次のようになります。

m <- as.xts(meto)

これにより、次の出力が得られます。

> str(m)
An ‘xts’ object from NA to NA containing:
  Data: num [1:10, 1] 30 25 25 25 20 20 20 20 20 20
  Indexed by objects of class: [chron,dates,times] TZ: 
  xts Attributes:  
 NULL
> summary(m)
     Index          m       
 Min.   :NA   Min.   :20.0  
 1st Qu.:NA   1st Qu.:20.0  
 Median :NA   Median :20.0  
 Mean   :NA   Mean   :22.5  
 3rd Qu.:NA   3rd Qu.:25.0  
 Max.   :NA   Max.   :30.0  
 NA's   :10                 
Warning message:
In data.row.names(row.names, rowsi, i) :
  some row.names duplicated: 2,3,4,5,6,7,8,9,10 --> row.names NOT used

ご覧のとおり、zoo の時系列には、chron オブジェクトによってインデックス付けされた大量のデータが含まれています。xtsただし、 を使用して時系列に変換すると、最初は問題as.xtsないように見えますが、コマンドは NA を表示し、 の概要をとstr比較すると、インデックスで 36,000 を超える NA が作成されていることがわかります!metom

なぜこれが起こっているのか、またはそれを解決するために何ができるのか、誰にも分かりますか?

4

1 に答える 1

4

問題は、インデックスがクラスであるということですchron。についてはほとんど知りませんchronが、AFAIKでは、RでPOSIX日時オブジェクトを使用することをお勧めしPOSIXctますPOSIXlt。

zooからクラスへxtsの変換のどこかでchron情報が破棄されます。

インデックスをクラスに変換POSIXctすると、問題が解決します。

index(meto) <- as.POSIXct(index(meto)) 
as.xts(meto)

                    [,1]
2005-01-13 00:00:00   30
2005-01-13 01:00:00   25
2005-01-13 01:59:59   25
2005-01-13 03:00:00   25
2005-01-13 04:00:00   20
2005-01-13 04:59:59   20
2005-01-13 06:00:00   20
2005-01-13 07:00:00   20
2005-01-13 07:59:59   20
2005-01-13 09:00:00   20

Rの日付と時刻のクラスの操作の詳細については、、または-を参照?DateTimeClassesし?POSIXctて?strptimeください。これらはすべて同じヘルプページにつながります。

編集

からインポートするときにオブジェクトxtsを処理することになっている場合は、関数にバグが見つかっている可能性があります。chronzooxts::xts

この行で問題が発生します。

if (inherits(order.by, "dates")) 
    index <- as.numeric(as.POSIXct(strptime(as.character(order.by), 
        "(%m/%d/%y %H:%M:%S)")))

chronしかし、あなたのオブジェクトのフォーマットが("d/m/y", "h:m:s")-私はあなたからこれを知っていることに注意してくださいstr(meto)。注意深く見てください-日と月の間にずれがあります。

これはロケールの問題である可能性があります。パッケージの作者は、標準フォーマットがm / d / yである米国に住んでいると思いますが、他の多くの場所では、標準フォーマットはd / m/yです。

したがって、どういうわけかzooとxtsの間の変換では、変換コードはユーザーのロケールに合わせて調整する必要があります。

この情報をパッケージの作成者に連絡することをお勧めします。

于 2011-12-09T14:21:25.800 に答える