Performing linear trend analysis within shapefile using ArcGIS Desktop?

Question

I have a shapefile (or several shapefiles, one per year, that I was going to join together to get all of the data in one attribute table) with an attribute for surface water area for each year. So like this: area_2000, area_2001, ..., area_2015; each year with about 2 million objects/rows.

I was hoping to do a linear trend time series analysis in ArcGIS Desktop to get the slope, p-value, etc. as additional attributes for each object. I have used ArcGIS Desktop before to do regression analysis, but this time I want to do time series to get the trend in area for each object.

Is there a tool that would support this?

Answer 1

This question has no GIS component, but anyway. I ran this field calculator expression:

import numpy as np
x=[2001,2002,2003,2004,2005]
def getSlope(y):
 a,b=np.polyfit(x, y, 1)
 return a
#----------------
getSlope([ !Y_2001!, !Y_2002!, !Y_2003!, !Y_2004!, !Y_2005!])

on the field "SLOPE" in below table:

To get a slope of a trend. Look at numpy functions closely.