<div dir="ltr">Hello all,<div><br></div><div>i am fairly new to R and just started using the ACS package.</div><div><br></div><div>I have been playing with the examples in the link here:</div><div><br></div><div><a href="http://eglenn.scripts.mit.edu/citystate/category/code/">http://eglenn.scripts.mit.edu/citystate/category/code/</a><br></div><div> </div><div>I have adapted the code to look at another table (foreign born individuals) by county  in Louisiana:  5 year ACS ending in 2011: B05006</div><div>In the example the code ultimately calculates a proportion which is then plotted on a map. I however want to create a dataset with the  estimate of the populations (no standard errors) of a few countries of origin (Liberia :HD01_VD111 and Nigeria: HD01_VD112 and Sierra Leone: HD01_VD112). also i would like to create a new variable say &quot;combined&quot; which is the sum of each of these variables.<br></div><div>I cannot seem to figure out how to create a dataset with these variables. i am unsure of how to manipulate the code:</div><div><br></div><div><div><br></div><div>us.pub.fb=divide.acs(numerator=us.fb[,111], </div><div>                 denominator=us.fb[,1], method=&quot;proportion&quot;)</div><div>pub.fb.est=data.frame(county=geography(us.pub.fb)[[1]], </div></div><div><br></div><div>to not calculate a proportion but to simply list the variables of interst in a dataframe. below is the full code:</div><div><br></div><div>any suggestions on how to accomplish this would be appreciated.</div><div><br></div><div>Thanks in advance for your help!</div><div>Nevin</div><div><br></div><div>PS..the full code below calculates the proportion:</div><div><div><br></div><div>us.county=geo.make(state=&quot;LA&quot;, county=&quot;*&quot;)</div><div><br></div><div><br></div><div>us.fb= acs.fetch(geography=us.county, table.number=&quot;B05006&quot;)</div><div><br></div><div><br></div><div>us.pub.fb=divide.acs(numerator=us.fb[,111], </div><div>                        denominator=us.fb[,1], method=&quot;proportion&quot;)</div><div>pub.fb.est=data.frame(county=geography(us.pub.fb)[[1]], </div><div>                         percent.pub.trans=as.numeric(estimate(us.pub.fb)))</div><div># this next step is all for Louisiana!</div><div>pub.fb.est$county=gsub(&quot;Parish&quot;, &quot;County&quot;, pub.fb.est$county)</div><div># clean up county names and find the states</div><div>pub.fb.est$state=gsub(&quot;^.*County, &quot;, &quot;&quot;, pub.fb.est$county)</div><div>pub.fb.est$county=gsub(&quot; County,.*&quot;, &quot;&quot;, pub.fb.est$county)</div></div></div>