4

毎日のストリームフローを含む csv ファイルがあります。毎日の値を毎月に組み合わせる必要があります。「hydroTSM」パッケージの「daily2monthly」機能を利用しようとしています。BRPT2.csv からのサンプル データ:

_date,_time,_value,_flag
 10/2/1959,0:00:00,0,2
 10/3/1959,0:00:00,0,2
 10/4/1959,0:00:00,1540,2
 10/5/1959,0:00:00,16100,2
 10/6/1959,0:00:00,6680,2
 10/7/1959,0:00:00,3100,2
 10/8/1959,0:00:00,2060,2

次のコマンドを使用しました。

qme<- read.csv(file = "BRPT2.csv",header = T,sep = ",") #read in csv file
date<- as.Date(qme$X_date,format("%Y-%m-%d")) #convert date column to date format from factor
flow<- qme[,3]
flow_2<-replace(flow,flow==-999,0) #replace the missing values (-999) with 0
df<- data.frame()
df<- rbind(df,data.frame(date,flow_2,stringsAsFactors = FALSE))
daily2monthly(df,FUN=sum,dates=1)

そして、次のエラーメッセージが表示されます。

rownames<-(のエラー*tmp*、値 = c("Oct-1959"、"Nov-1959"、"Dec-1959"、: 次元のないオブジェクトに 'rownames' を設定しようとしました

誰かがこれに関して私を助けてくれませんか。前もって感謝します。

4

1 に答える 1

5

df最初にあなたをzooオブジェクトに変換してみてください：

z <- zoo(df[, -1], df[, 1])
daily2monthly(z, FUN=sum, dates=1)
# 1959-10-01 1959-11-01 1959-12-01 1960-01-01 1960-02-01 1960-03-01 1960-04-01 
#    31440.0      411.9     1199.3     4373.0     1466.0     1904.0      741.0

のdata.frame方法はdaily2monthly、データの列が複数ある場合にのみ機能します。あなたの場合、1 つしかないため、関数はそのエラー メッセージを返します。もう 1 列のデータがあるとします。

df$station_3 <- sample(1:1000, nrow(df))
daily2monthly(df, FUN=sum, dates=1)
#           flow_2 station_3
# Oct-1959 31440.0    132423
# Nov-1959   411.9    149305
# Dec-1959  1199.3    157622
# Jan-1960  4373.0    176413
# Feb-1960  1466.0    143373
# Mar-1960  1904.0    166102
# Apr-1960   741.0    141861
于 2016-09-21T06:19:27.487 に答える