XTS

Question

Iは、2つの列のデータが設定されているを使用して、不規則な毎時データから毎日の平均値を計算します。時系列は不規則な毎時のデータです、私は毎日の平均データを計算する必要があります。

私の問題は、最終的なファイルではなく、全体Ozone列の一つの値だけを持っています。私はいくつかの異なるバージョンを試しましたが、私は問題を解決できません。ここで

dataset1 <- structure(list(Date = structure(1:3, .Label = c("7/11/2013 18:33", 
    "7/11/2013 18:34", "7/11/2013 18:35"), class = "factor"), ozone = c(40.7, 40.4, 
    40.9)), .Names = c("Date", "ozone"), row.names = c(NA, 3L), class = "data.frame") 
#convert dataframe to xts object 
library(xts) 
xt1 <- xts(dataset1[,-1], 
      order.by = as.POSIXct(dataset1$Date, format= "%m/%d/%Y %H:%M")) 
x_updated <- apply.daily(xt1, colMeans) 
#convert back the dataset to data.frame 
write.csv(as.data.frame(x_updated), file="daily") 

Answer 1

はdplyrソリューションです：

library(dplyr) 

df2 <- dataset1 %>% 
    mutate(date2 = as.Date(dataset1$Date)) %>% 
    group_by(date2) %>% 
    mutate(avg_ozone = mean(ozone)) %>% 
    select(Date, ozone, avg_ozone) 

Answer 2

# Get 'day' from 'datetime' 
# dataset1$day = as.POSIXct(dataset1$Date, format='%d/%m/%Y') # Apparently this caused problem 
dataset1$day = do.call(rbind, strsplit(as.character(dataset1$Date)," "))[,1] 

# Get Daily Averages 
daily_averages = aggregate(ozone~day, dataset1, mean) 

Answer 3

私はタイプミスを修正した後、あなたの出力に問題があるだろうと私に示唆してあなたの再現性のある例では何も（ありませんあなたのas.POSIXctコールへのformat引数 - 元は"%m\%d\%y %H:%M"でした）。つまり、あなたの例は1日のデータしか持っていなかったので、ファイルには1行の出力しかないはずです。ここで

は「毎日」ファイルに複数の観測があります示して再現性の例です。

# re-create example data 
set.seed(21) 
Data <- data.frame(Date=seq(ISOdate(2017, 1, 1), ISOdate(2017, 1, 3), by="hours")) 
Data$Ozone <- rnorm(NROW(Data)) 
Data$Date <- format(Data$Date, "%m/%d/%Y %H:%M") 
# convert to xts 
library(xts) 
xt1 <- xts(Data[,-1], as.POSIXct(Data$Date, format="%m/%d/%Y %H:%M")) 
x_updated <- apply.daily(xt1, colMeans) 
# convert to data.frame 
write.csv(as.data.frame(x_updated)) 
# "","V1" 
# "2017-01-01 23:00:00",0.038108046611944 
# "2017-01-02 23:00:00",-0.0380170075916872 
# "2017-01-03 12:00:00",0.612294079665861