0

私は次の機能を持っています。 CreateChronVectorそれが意味することを正確に行います。結果のベクトルは、デフォルトでは 1 時間間隔です。このRoundHour関数は、chron ベクトルを時間に切り上げます。

CreateChronVector  <-  function(chronFrom, chronTo, frequency = "hourly")  {
  library(chron)
  datesFrom <- dates(chronFrom)
  timesFrom <- (chronFrom - dates(chronFrom))
  datesTo <- dates(chronTo)
  timesTo <- (chronTo - dates(chronTo))
  if ((timesFrom != 0 || timesTo != 0) && frequency == "daily") {
      print("Error: The indicated dates have hour components while the given frequency is daily.")
  }
  else {
      if (timesTo == 0 && frequency == "hourly") {
          timesTo <- 23/24
      }
      if (frequency == "hourly") {
          chronFrom <- chron(dates = datesFrom, times = timesFrom, 
              format = c(dates = "m/d/y", times = "h:m:s"))
          chronTo <- chron(dates = datesTo, times = timesTo, 
              format = c(dates = "m/d/y", times = "h:m:s"))
          dateVector <- seq(chronFrom, chronTo, by = 1/24)
      }
      else if (frequency == "daily") {
          dateVector <- seq(datesFrom, datesTo)
      }
      return(dateVector)
  }
}

RoundHour  <- function(x)  {
  res <-  trunc(x,'hours', eps=1e-17)
  res <-  ifelse((x-res) > 0.5/24, res+1/24, res)
  return(as.chron(res))
}

私が直面している問題は、間隔が一貫していないことです。例として、以下のコードは 2 つの異なる間隔サイズを返します。

unique(diff(CreateChronVector(as.chron('2010-01-01'), as.chron('2010-01-01'))))

同様に、丸め関数を使用しても問題は解決しません。

unique(diff(RoundHour(CreateChronVector(as.chron('2010-01-01'), as.chron('2010-01-01')))))

この問題は丸め誤差に関係していると確信しています。trunc 関数とその eps パラメータを試してみましたが、うまくいきませんでした。

4

2 に答える 2

0

@Gからポイントを取ります。Grothendieckさん、これを試してみれば、彼が何について話しているかを見ることができます:

hours <- 1:23
dateVector <- sapply(hours , function(x){ chron( dates = "01/01/10" , times = paste0(x,":00:00") ) } )
head( dateVector )
[1] 14610.04166666666606034 14610.08333333333393966 14610.12500000000000000
[4] 14610.16666666666606034 14610.20833333333393966 14610.25000000000000000
unique(diff(dateVector))
[1] 0.04166666666787932626903 0.04166666666606033686548

これらの数値は浮動小数点で正確に表現できないため、実際にはできませんが、これが重要な理由はありますか?

于 2013-02-26T00:22:21.330 に答える