I have a large database of point data throughout all of Eastern Australia. Each point represents species presence/absence data and has attached about 18 predictor variables. To account for sampling intensity, I would like to make a 1km x 1km grid and count the frequency of the number of presence records (this would be my response variable) and create another column in the data frame that is the sum of all points within the cell (this would become an offset variable in a GLMM).
What is the best approach to this problem? I know of the Point to Raster tool, however I dont think it allows for keeping the other predictor variables. Aggregate Points is another option, however I am using ArcGIS 10.3 and I dont have access to the ArcGIS GeoAnalytics Server.