2017-11-27 24 views
0

unsolved problemへの解決策を探していたら、私はacsパッケージに出くわしました。私は、choroplethパッケージには、[city, state]の形式のデータから郡情報を取得する方法がないと仮定します。そのため、acsを使った前処理が必要です。Rパッケージ "acs":郡名、FIPSを取得しますか?

私は都市の郡の情報を取得するには、次のコードを試してみました:

library(acs) 
geo.lookup(state="CA", place="San Francisco") 

> geo.lookup(state="CA", place="San Francisco") 
    state state.name   county.name place    place.name 
1  6 California     <NA> NA      <NA> 
2  6 California San Francisco County 67000  San Francisco city 
3  6 California  San Mateo County 73262 South San Francisco city 

私たちが知っているように、都市は異なる郡の一部とすることができます。おそらく、私は第二に行くだろう

> geo.lookup(state="CA", place="San Francisco")[2,] 

    state state.name   county.name place   place.name 
2  6 California San Francisco County 67000 San Francisco city 

デフォルトでは、

私の質問: 状態の略語、郡名、郡FIPSも取得できますか?私はdocumentationで答えを見つけることができませんでした。

また、さらなる処理(choroplethrとのマッチング)のために、county.nameの最後の "County"とplace.nameの "city"を削除する必要があります。

答えて

1

州の省略名、郡名、郡FIPSを例に追加する方法は次のとおりです。 Rには、状態名と状態略語の変数が組み込まれています。 FIPSコードについては、国勢調査局のウェブサイトからcsvファイルを読みました。

library(acs) 
library(tidyverse) 

states <- cbind(state.name, state.abb) %>% tbl_df() 

fips <- 
    read_csv(
    "https://www2.census.gov/geo/docs/reference/codes/files/national_county.txt", 
    col_names = c("state.abb", "statefp", "countyfp", "county.name", "classfp") 
) 

query <- geo.lookup(state = "CA", place = "San Francisco")[2, ] %>% 
    tbl_df() %>% 
    left_join(states, by = "state.name") %>% 
    left_join(fips, by = c("county.name", "state.abb")) 

query 

# # A tibble: 1 x 9 
# state state.name   county.name place   place.name state.abb statefp countyfp classfp 
# <chr>  <chr>    <chr> <int>    <chr>  <chr> <chr> <chr> <chr> 
# 1  6 California San Francisco County 67000 San Francisco city  CA  06  075  H6 

あなたがあなたの質問の最後に注意点として、あなたはそれがchoroplethrに合うようにするためにもう少し、このデータをクリーンアップする必要があるかもしれません。

+0

これはかなりうまくいきます。私はstatefp + countyfpをマージし、choroplethrに一致させるために最初の "0"を削除する必要があるように見えます。最後に、これを使ってMSAのfipsコードも入手できますか?改善のアイデア:Census csv: 'install.packages( 'tigris') ライブラリ(チグリス) データ(fips_codes)' – Christopher

+0

(実際には、MSA/CBSAシェイプを視覚化できるかどうかは現在わかりません。 choroplethrZipの 'zip.regions.'で' cbsa'を実行しましたが、例は見つかりませんでした。) – Christopher

関連する問題