library(PracTools)
library(dplyr)
g <- GeoDistPSU(Test_Data_US$lat,
Test_Data_US$long,
"miles",
100,
Input.ID = Test_Data_US$ID)
Test_Data_US <- inner_join(Test_Data_US, g$PSU.ID, by=c("ID" = "Input.file.ID"))
GeoMinMOS(lat = Test_Data_US$lat,
long = Test_Data_US$long,
geo.var = Test_Data_US$psuID,
MOS.var = Test_Data_US$Amount,
MOS.min = 200000)
Run the code above in your browser using DataLab