2012-05-12 177 views
2

對不起,這個愚蠢的問題,我是新來R.我有一些這樣的甲酸,並保存在CSV:描述統計

%Y-%m-%d,st1,st2,st3,st4,st5,st6,st7,st8,st9,st10 
2005-09-20,38.75,48.625,48.5,23.667,45.5,48.75,18.75,33.25,43.455,76.042 
2005-09-21,39.482,49.3,49,23.9,46.15,50.281,18.975,34.125,44.465,78.232 
... 

我它導入[R

library(fPortfolio) 

Data <- readSeries(file = "data.csv", header = TRUE, sep = ",") 

我希望有一些描述性統計

library(psych) 
describe(Data) 

Error in x[!is.na(x[, i]), i] : 
    invalid or not-yet-implemented 'timeSeries' subsetting 

什麼建議嗎?

回答

1

你可能想讓它成爲一個時間系列第一對不對?

tS <- dummySeries() #make quick dummy time series 
describe(tS) # fails 

newtS<-as.ts(tS) 
describe(newtS) #works fine giving: 

     var n mean sd median trimmed mad min max range skew kurtosis se 
Series 1 1 12 0.49 0.25 0.44 0.48 0.29 0.13 0.89 0.76 0.24 -1.52 0.07 
Series 2 2 12 0.45 0.28 0.44 0.45 0.42 0.07 0.83 0.77 0.03 -1.74 0.08 
+0

謝謝!有用!在R中定義了幾個不同的時間序列結構,比如zoo(),as.ts(),timeSeries,它們有什麼不同? – onethird

+0

不太清楚差異,但有趣的是'as.zoo(tS):describe(tS)'也可以正常工作。如果你打印as.ts(tS)'和'as.zoo(tS)',兩者的結構看起來有點不同。 'as.ts'給出'2012年1月',但'動物園'給出'2012-01-01'。很高興它的作品 –