See here about the Mean Center of the US Population https://en.wikipedia.org/wik
ID: 3824905 • Letter: S
Question
See here about the Mean Center of the US Population https://en.wikipedia.org/wiki/Mean_center_of_the_United_States_population
See http://www.cs.iit.edu/~cs104/labs/Lab10/populationByState.xls for the 2010 census population by state also including the latitude and longitude for each state.
STNAME: state name
POPULATION: 2010 Census population tabulated for the state
LATITUDE: latitude coordinate for the center of population for the state
LONGITUDE: longitude coordinate for the center of population for the state
Write a Matlab script called populationMean to read the file and calculate the Mean Center of the US Population for the 50 states. Output the latitude and longitude of the mean center.
Then use Matlab code to add in these territories and see how the mean center moves:
District of Columbia population=601723 latitude=38.91027 longitude=-77.014468
Puerto Rico population=3725789 latitude=18.280422 longitude=-66.346513
Look up the latitude and longitude of both answers and see which city is nearest to the mean center.
Explanation / Answer
data(boston) lw <- nb2listw(boston.soi) can.sim <- spdep:::can.be.simmed(lw) env <- new.env(parent=globalenv()) assign("listw", lw, envir=env) assign("can.sim", can.sim, envir=env) assign("similar", FALSE, envir=env) assign("verbose", FALSE, envir=env) assign("family", "SAR", envir=env) eigen_setup(env) get("similar", envir=env) do_ldet(0.5, env) rm(env) env <- new.env(parent=globalenv()) assign("listw", lw, envir=env) assign("can.sim", can.sim, envir=env) assign("similar", FALSE, envir=env) assign("verbose", FALSE, envir=env) assign("family", "SAR", envir=env) assign("n", length(boston.soi), envir=env) eigen_pre_setup(env, pre_eig=eigenw(similar.listw(lw))) do_ldet(0.5, env) rm(env) env <- new.env(parent=globalenv()) assign("listw", lw, envir=env) assign("can.sim", can.sim, envir=env) assign("similar", FALSE, envir=env) 38 do_ldet assign("family", "SAR", envir=env) assign("n", length(boston.soi), envir=env) Matrix_setup(env, Imult=2, super=FALSE) get("similar", envir=env) do_ldet(0.5, env) rm(env) env <- new.env(parent=globalenv()) assign("listw", lw, envir=env) assign("n", length(boston.soi), envir=env) assign("can.sim", can.sim, envir=env) assign("similar", FALSE, envir=env) assign("family", "SAR", envir=env) spam_setup(env) get("similar", envir=env) do_ldet(0.5, env) rm(env) env <- new.env(parent=globalenv()) assign("listw", lw, envir=env) assign("n", length(boston.soi), envir=env) assign("similar", FALSE, envir=env) assign("family", "SAR", envir=env) LU_setup(env) get("similar", envir=env) do_ldet(0.5, env) rm(env) env <- new.env(parent=globalenv()) assign("listw", lw, envir=env) assign("n", length(boston.soi), envir=env) assign("similar", FALSE, envir=env) assign("family", "SAR", envir=env) LU_prepermutate_setup(env) get("similar", envir=env) do_ldet(0.5, env) rm(env) env <- new.env(parent=globalenv()) assign("listw", lw, envir=env) assign("similar", FALSE, envir=env) assign("family", "SAR", envir=env) cheb_setup(env, q=5) get("similar", envir=env) do_ldet(0.5, env) rm(env) env <- new.env(parent=globalenv()) assign("listw", lw, envir=env) assign("n", length(boston.soi), envir=env) assign("similar", FALSE, envir=env) assign("family", "SAR", envir=env) set.seed(12345) mcdet_setup(env, p=16, m=30) get("similar", envir=env) do_ldet(0.5, env) rm(env)