[R] Households per Census block

[Zack Almquist]

P.S. The US census has different "populations" (or worlds) so make sure the
housing variable you use is accessing the correct "world."

Best,

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

On Tue, Aug 4, 2015 at 2:12 PM, Zack Almquist <almquist at umn.edu> wrote:

> Hi Keith,
>
>
> On Tue, Aug 4, 2015 at 12:43 PM, Keith S Weintraub <kw1958 at gmail.com>
> wrote:
>
>> 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’
>>
>
> Sorry that is my bad, I didn't verify the variable name at (
> http://api.census.gov/data/2010/sf1/variables.html). This seems to work
> for me, as a quick test:
>
> housing<-CensusAPI2010(variables="H0060001", state.fips="06", level =
> c("tract"), key, summaryfile = c("sf1"))
>
> So the larger example:
>
>  ## 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("H0060001"), state.fips=state.fips, level =
> c("block"), key, summaryfile = c("sf1"))
>
> should work just fine.
>
> Best,
>
> -- Zack
> ---------------------------------------------------------
> Zack W.  Almquist
> Assistant Professor
> Department of Sociology and School of Statistics
> Affiliate, Minnesota Population Center
> University of Minnesota
>
>
>> 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.
>> >
>> >
>>
>>
>

	[[alternative HTML version deleted]]