【问题标题】:Converting xts to ts: Error in Round(frequency)将 xts 转换为 ts:回合误差(频率)
【发布时间】:2011-10-09 17:03:30
【问题描述】:

我有一些已导入的 csv 数据已转换为 xts 对象。如果我尝试将其转换为 ts 对象(最终目标是使用 acf 之类的函数),我会得到:

“回合错误(频率):数学的非数字参数 函数”

转换它的代码是:

library("zoo")
#Working With Milliseconds
op <- options(digits.secs=3)

#Rename Function
clean_perfmon = function(x, servername) {
  names(x)[names(x)=="X.PDH.CSV.4.0...Coordinated.Universal.Time..0."] <- "Time"
  x$Time = strptime(x$Time, "%m/%d/%Y %H:%M:%OS")
  return(x)
}

web02 = read.csv("/home/kbrandt/Desktop/Shared/web02_2011_07_20_1.csv")
web02 = clean_perfmon(web02, "NY.WEB02")
web02ts = xts(web02[,-1], web02[,"Time"])

时间大多是规律的,但 MS 有一些变化:

time(web02ts)[1:3]
[1] "2011-07-20 11:21:50.459 EDT" "2011-07-20 11:21:51.457 EDT" "2011-07-20 11:21:52.456 EDT"

部分数据有 NA 点:

> web02ts[1:3,1]
                        X..NY.WEB02.Process.Idle....Processor.Time
2011-07-20 11:21:50.459                                         NA
2011-07-20 11:21:51.457                                   1134.819
2011-07-20 11:21:52.456                                   1374.877

更新:
更改为每秒分辨率,非 na 子集无济于事:

> as.ts(web02ts[2:10,1])
Error in round(frequency) : Non-numeric argument to mathematical function
> web02ts[2:10,1]
                    X..NY.WEB02.Process.Idle....Processor.Time
2011-07-20 11:21:51                                   1134.819
2011-07-20 11:21:52                                   1374.877
2011-07-20 11:21:53                                   1060.842
2011-07-20 11:21:54                                   1067.092
2011-07-20 11:21:55                                   1195.205
2011-07-20 11:21:56                                   1223.328
2011-07-20 11:21:57                                   1121.774
2011-07-20 11:21:58                                   1187.393
2011-07-20 11:21:59                                   1378.001
>

另外,frequency(web02ts) 返回NULL

【问题讨论】:

    标签: r time-series xts


    【解决方案1】:

    strptime 创建一个 POSIXlt 类的对象。 as.ts 不支持它,并认为它是一个列表,因此抱怨非数字参数。改为POSIXct

    as.POSIXct(strptime(x$Time, "%m/%d/%Y %H:%M:%OS"))
    

    【讨论】:

      【解决方案2】:

      xts/zoo 对象必须是规则的,才能具有非 NULL 频率。

      您没有显示如何更改为每秒分辨率,但如果您通过options(digits.secs=0) 尝试,那将不起作用,因为它只会影响打印。你需要做这样的事情:

      # example data
      set.seed(21)
      web02ts <- xts(rnorm(10), Sys.time()+1:10+runif(10)/3)
      
      web02ts_reg <- align.time(web02ts,1)
      frequency(web02ts_reg)
      # [1] 1
      as.ts(web02ts_reg)
      # Time Series:
      # Start = 1 
      # End = 10 
      # Frequency = 1 
      #  [1]  0.793013171  0.522251264  1.746222241 -1.271336123  2.197389533
      #  [6]  0.433130777 -1.570199630 -0.934905667  0.063493345 -0.002393336
      

      【讨论】:

        猜你喜欢
        • 2017-12-17
        • 1970-01-01
        • 1970-01-01
        • 2017-07-04
        • 1970-01-01
        • 2018-06-20
        • 2010-12-01
        • 2012-08-27
        • 2020-11-06
        相关资源
        最近更新 更多