data('rockland')
data('orange')
data('nrcsd')
o_and_r <- rbind(orange, rockland)
o_and_r <- o_and_r %>%
geo_filter(nrcsd) %>%
geo_trim(nrcsd)
adj <- adjacency(o_and_r)
seam_rip(adj, o_and_r, 'county', c('071', '087'))
Run the code above in your browser using DataLab