3

編集：これまでに回答してくれた人に感謝します。私は R の初心者で、修士論文の大規模なプロジェクトに取り掛かったばかりなので、最初の処理に少し圧倒されています。私が使用しているデータは次のとおりです (WMO が公開している降雨データから)。


120 6272100 KHARTOUM 15.60 32.55 382 1899 1989 0.0
1899 0.03 0.03 0.03 0.03 0.03 1.03 13.03 12.03 9999 6.03 0.03 0.03
1900 0.03 0.03 0.03 0.03 0.03 23.03 80.03 47.03 23.03 8.03 0.03 0.03
1901 0.03 0.03 0.03 0.03 0.03 17.03 23.03 17.03 0.03 8.03 0.03 0.03
(...)
120 6272101 JEBEL AULIA 15.20 32.50 380 1920 1988 0.0
1920 0.03 0.03 0.03 0.00 0.03 6.90 20.00 108.80 47.30 1.00 0.01 0.03
1921 0.03 0.03 0.03 0.00 0.03 0.00 88.00 57.00 35.00 18.50 0.01 0.03
1922 0.03 0.03 0.03 0.00 0.03 0.00 87.50 102.30 10.40 15.20 0.01 0.03
(...)

There are ~100 observation stations that I'm interested in, each of which has a varying start and end date for rainfall measurements. They're formatted as above in a single data file, with stations separated by "120 (station number) (station name)".

I need first to separate this file by station, then to extract March, April, May and June for each year, then take a total of these months for each year. So far I'm messing around with loops (as below), but I understand this isn't the right way to go about it and would rather learn some better technique. Thanks again for the help!

(Original question:) I've got a large data set containing rainfall by season for ~100 years over 100+ locations. I'm trying to separate this data into more managable arrays, and in particular I want to retrieve the sum of the rainfall for March, April, May and June for each station for each year. The following is a simplified version of my code so far: 

a <- array(1,dim=c(10,12))
for (i in 1:5) {

  all data:
  assign(paste("station_",i,sep=""), a)

  #march - june data:
  assign(paste("station_",i,"_mamj",sep=""), a[,4:7])
}

So this gives me station_(i)__mamj_ which contains the data for the months I'm interested in for each station. Now I want to sum each row of this array and enter it in a new array called station_(i)_mamj_tot. Simple enough in theory, but I can't work out how to reference station_(i)_mamj so that it varies the value of i反復ごとに。どんな助けでも大歓迎です！

4

3 に答える 3

2

元の質問には、get() を使用します。

i <- 10
var <- paste("test", i, sep="_")
assign(10, var)
get(var)

Davidが言ったように、これはおそらく最善の道ではありませんが、時には役立つことがあります（IMOの割り当て/取得構造はeval（parse）よりもはるかに優れています）。

于 2012-05-14T17:44:36.317 に答える
1

、、などのassign変数を作成するために使用しているのはなぜですか？、、、などのリストにそれらを格納する方がはるかに簡単で直感的です。次に、それぞれのインデックスを使用してアクセスできます。station1station2station_3_mamjstations[[1]]stations[[2]]stations_mamj[[3]]

使用しているステーションごとのデータは同じサイズの行列であるように見えるため、3次元の行列として扱うこともできます。

ETA：ちなみに、本当にこの方法で問題を解決したい場合は、次のようにします。

eval(parse(text=paste("station", i, "mamj", sep="_")))

ただし、使用しないことevalはほとんどの場合悪い習慣であり、データに対して単純な操作を行うことさえ困難になります。

于 2012-05-14T17:21:40.913 に答える