[R] Households per Census block

[Zack Almquist]

Hi Anthony and Keith Weintraub,

Here is a way to do what you are asking using the UScensus2010 packages:

## latest version of the package, not yet on CRAN
install.packages("UScensus2010", repos="http://R-Forge.R-project.org")
library(UScensus2010)
install.blk()
library(UScensus2010blk)
### You will want the H0010001 variable (see help(alabama.blk10))
### Other variables are also available
### You can use the new api function in UScensus2010 to get arbitrary
variables from SF1 and acs

data(states.names)
head(states.names)
state.blk.housing<-vector("list",length(states.names))
## notice this could be greatly spead up using the library(parallel)
## with mclapply
## This will be somewhat slow b/c of so much spatial data
for(i in 1:length(states.names)){
data(list=paste(states.names[i],"blk10",sep="."))
temp<-get(paste(states.names[i],"blk10",sep="."))
#unique b/c more shapefiles than fips
state.blk.housing[[i]]<-unique(temp at data[,c("fips","H0010001")])
print(i)
rm(paste(states.names,"blk10",sep="."))
}

###########
# alternatively Using the US Census API function in the new UScensus2010
package
###########

## Get all states fips code
data(countyfips)
state.fips<-unique(substr(countyfips$fips,1,2))
head(state.fips)
length(state.fips) ## will be 51=50 (states)+ 1(DC)
## You will need a census key
key<-"YOUR KEY HERE"
housing<-CensusAPI2010(c("H0010001"), state.fips=state.fips, level =
c("block"), key, summaryfile = c("sf1"))

Best,

-- Zack
---------------------------------------------------------
Zack W.  Almquist
Assistant Professor
Department of Sociology and School of Statistics
Affiliate, Minnesota Population Center
University of Minnesota


On Mon, Aug 3, 2015 at 12:43 PM, Anthony Damico <ajdamico at gmail.com> wrote:

> hi, ccing the package maintainer.  one alternative is to pull the HU100
> variable directly from the census bureau's summary files: that variable
> starts at position 328 and ends at 336.  just modify this loop and you'll
> get a table with one-record-per-census-block in every state.
>
>
> https://github.com/davidbrae/swmap/blob/master/how%20to%20map%20the%20consumer%20expenditure%20survey.R#L104
>
> (1) line 134 change the very last -9 to 9
> (2) line 137 between "pop100" and "intptlat" add an "hu100"
>
>
> summary file docs-
>
> http://www.census.gov/prod/cen2010/doc/sf1.pdf#page=18
>
>
>
> On Mon, Aug 3, 2015 at 11:55 AM, Keith S Weintraub <kw1958 at gmail.com>
> wrote:
>
>> Folks,
>>
>> I am using the UScensus2010 package and I am trying to figure out the
>> number of households per census block.
>>
>> There are a number of possible data downloads in the package but
>> apparently I am not smart enough to figure out which data-set is
>> appropriate and what functions to use.
>>
>> Any help or pointers or links would be greatly appreciated.
>>
>> Thanks for your time,
>> Best,
>> KW
>>
>> ______________________________________________
>> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
>> https://stat.ethz.ch/mailman/listinfo/r-help
>> PLEASE do read the posting guide
>> http://www.R-project.org/posting-guide.html
>> and provide commented, minimal, self-contained, reproducible code.
>>
>
>

	[[alternative HTML version deleted]]