将xts转换为ts:Round(频率)出错

时间:2011-07-20 19:59:05

标签: r time-series xts

我有一些导入的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

2 个答案:

答案 0 :(得分:3)

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

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

答案 1 :(得分:1)

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