Added gis.R

So far, includes only majority_aggregate