[R] Households per Census block

[Keith S Weintraub]

I had to download a bunch of stuff but I got it mostly working.

Unfortunately using the alternative method I get the following:

> housing<-CensusAPI2010(c("H0010001"), state.fips=state.fips, level = c("block"), key, summaryfile = c("sf1"))
Error in file(con, "r") : cannot open the connection
In addition: Warning message:
In file(con, "r") : cannot open: HTTP status was '400 Bad Request’

I have a feeling that this is not a problem with the API.

Thanks for your help,
KW







> On Aug 3, 2015, at 2:12 PM, Zack Almquist <almquist at umn.edu> wrote:
> 
> 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.
> 
>