2012-04-10 37 views
1

我有一個名稱數組,我想使用這些名稱作爲數據框的列名,但我收到了一些錯誤。我不確定究竟該如何做到這一點,但這是我目前爲止所做的。來自數組值的數據幀列名稱

windspeeds = data.frame() 
cities <- c("albuquerque_nm", "boston_ma", "charlotte_nc", "dallas_tx", "denver_co", "helena_mt", "louisville_ky", "pittsburgh_pa", "salt_lake_city_ut", "seattle_wa") 
for(i in 1:10){ 
    fastest <- read.delim(paste("http://www.itl.nist.gov/div898/winds/data/nondirectional/datasets/", cities[i], ".prn", sep=""), col.names=c("NULL", "fastest", "NULL", "NULL"), skip=4, header=F, sep="")$fastest 
    windspeeds$cities[i] = fastest 
} 

我收到此錯誤:

Error in `$<-.data.frame`(`*tmp*`, "cities", value = 59L) : 
replacement has 1 rows, data has 0 
In addition: Warning message: 
In windspeeds$cities[i] = fastest : 
number of items to replace is not a multiple of replacement length 

我必須將數組轉換爲某種類型的字符串或恆定的?

+0

循環中的對象data是什麼?嗯,你已經更新了... – Ben 2012-04-10 03:24:43

回答

3

您的問題之一是您的查詢不會爲每個城市返回相同數量的記錄(免責聲明,我對您的數據或它應該看起來不瞭解)。無論如何,這裏有一個方法來讀取數據到一個列表對象,它可能是一個更「R-ISH」的方式來做事:

x <- lapply(cities, function(x) 
    read.delim(paste("http://www.itl.nist.gov/div898/winds/data/nondirectional/datasets/", x, ".prn", sep=""), 
      col.names=c("NULL", "fastest", "NULL", "NULL"), skip=4, header=FALSE, sep="")$fastest 
      ) 

X現在看起來像:

> str(x) 
List of 10 
$ : int [1:46] 59 49 51 52 57 52 45 54 49 64 ... 
$ : int [1:42] 50 55 79 56 53 41 51 51 65 62 ... 
$ : int [1:29] 33 42 40 52 42 48 52 51 46 51 ... 
$ : int [1:32] 48 45 46 45 43 53 43 58 46 43 ... 
$ : int [1:33] 42 51 49 44 47 47 50 44 44 54 ... 
$ : int [1:48] 58 58 58 58 70 55 56 62 59 70 ... 
$ : int [1:39] 40 39 50 53 50 51 54 54 51 50 ... 
$ : int [1:18] 47 56 60 44 54 50 42 52 47 47 ... 
$ : int [1:46] 53 49 40 53 55 40 49 46 61 41 ... 
$ : int [1:10] 38 44 35 46 42 45 41 45 42 43 

而且具有描述性統計:

> do.call(rbind, lapply(x, summary)) 
     Min. 1st Qu. Median Mean 3rd Qu. Max. 
[1,] 45 49.50 53.0 55.02 57.00 85 
[2,] 41 49.25 54.5 56.26 60.75 85 
[3,] 33 39.00 42.0 44.86 51.00 65 
[4,] 39 45.75 48.0 49.16 51.50 67 
[5,] 42 44.00 48.0 48.67 51.00 61 
[6,] 42 49.00 55.0 54.04 58.00 71 
[7,] 38 43.50 49.0 48.74 52.50 66 
[8,] 39 45.00 47.0 48.44 53.50 60 
[9,] 40 45.25 49.0 50.41 54.00 69 
[10,] 35 41.25 42.5 42.10 44.75 46 

無論你是否應該有相同數量的每個城市的記錄是未知的,但希望這將讓你正確的道路。

+0

感謝您的提示,我沒有編程太多的R和不知道任何公約。 – trev9065 2012-04-10 12:49:50

+0

@ trev9065 - 別擔心!每個人都必須在某個時候學習。 R與許多您可能熟悉的編程語言絕對不同......一旦您開始思考R,您會驚喜於能用相對少量的代碼完成多少工作! – Chase 2012-04-10 16:55:48

+0

有沒有一種方法來命名行並將它們與城市相匹配? – trev9065 2012-04-11 14:01:39

相關問題