[R] Households per Census block

[Zack Almquist]

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]]