I got the toolbox to work in 10.1 by updating the script using comments and links provided by @RyanDalton.
If you use the code below within the existing toolbox created by Ryan Dalton (http://resources.arcgis.com/gallery/file/geoprocessing/details?entryID=A0CD1751-1422-2418-882E-001EE0DC0D35) it should work in 10.1. I tested the tool and it worked for me. The only thing I have noticed is that the case field has to be a string and cannot be too long. Plenty of ways to workaround that though. Enjoy!
#-------------------------------------------------------------------------------
# Name: MedianStats.py
# Purpose: Calculate the Median value from a data layer or table, and save
# the results in a separate table. Allows you to calculate the median
# value for a case. Updated for use in 10.1
#
# Author: P. Johnson modified a script written by R.Dalton for use in 10.1.
# R. Dalton's script was modified version of a script by C.Denninger.
# C. Denninger's script was a modified version of a script by Chris Snyder
#
# Modified: July 3, 2013
#-------------------------------------------------------------------------------
#!/usr/bin/env python
def uniqify(seq):
# Not order preserving
uset = set(seq)
return list(uset)
import sys, arcpy, string, os
# Load required toolboxes...
inputTable = arcpy.GetParameterAsText(0)
statField = arcpy.GetParameterAsText(1) #this needs to be a numeric field
caseFields = arcpy.GetParameterAsText(2) #this needs to be a double or float
outTable = arcpy.GetParameterAsText(3)
arcpy.AddMessage("Analyzing median values for Table = " + inputTable + ":")
arcpy.AddMessage(" Statistics field = " + statField + " & Case field = " + caseFields)
# Needed outTable Schema = OID,caseField,Frequency,Med_statField
outFolder = os.path.dirname(outTable)
outFile = os.path.basename(outTable)
arcpy.CreateTable_management(outFolder, outFile, "", "")
caseList = string.split(caseFields,";")
if len(caseList) >1:
#arcpy.AddMessage("more than 1 case fields")
for cI in caseList:
##Output table field name
outCaseField = cI[:10]
##Only add the case field if it was selected
if outCaseField <> "":
arcpy.AddField_management(outTable, outCaseField, "TEXT", "", "", "30", "", "NULLABLE", "NON_REQUIRED", "")
else:
##Output table fields
outCaseField = caseFields[:10]
##Only add the case field if it was selected
if outCaseField <> "":
arcpy.AddField_management(outTable, outCaseField, "TEXT", "", "", "30", "", "NULLABLE", "NON_REQUIRED", "")
##Add the following fields in all cases
outCountField = "Count"
arcpy.AddField_management(outTable, outCountField, "LONG", "", "", "", "", "", "NON_REQUIRED", "")
outMedField = "MED_" + statField[:5]
arcpy.AddField_management(outTable, outMedField, "DOUBLE", "", "", "9", "", "", "NON_REQUIRED", "")
##Delete the "Field1" field created by default
arcpy.DeleteField_management (outTable, "Field1")
if caseFields == "":
# Read the values from the Source Table
valueList = []
searchRows = arcpy.SearchCursor(inputTable)
for searchRow in searchRows:
searchRowValue = searchRow.getValue(statField)
if searchRowValue == None:
pass #don't add a null value to the list!
else:
valueList.append(searchRowValue)
del searchRow
del searchRows
valueList.sort()
listLength = len(valueList)
if listLength == 0:
print "Every value was null! Exiting script..."; sys.exit()
elif listLength % 2 == 0: #even, so get the mean of the 2 center values
##medianValue = (valueList[listLength / 2]) + valueList[listLength / 2 - 1]) / 2.0
medianValue = valueList[listLength / 2]
else: #odd, so it's easy!
medianValue = valueList[listLength / 2]
#Write the values to the output table
listCount = listLength
insertRows = arcpy.InsertCursor(outTable)
insertRow = insertRows.NewRow()
insertRow.setvalue(outCountField, listCount)
insertRow.setvalue(outMedField, medianValue)
insertRows.InsertRow(insertRow)
del valueList
del listCount
del medianValue
elif caseFields <> "":
# Read the values from the Source Table
cList = []
caseRows = arcpy.SearchCursor(inputTable)
caseRow = caseRows.next()
while caseRow:
cFList = ""
for cF in string.split(caseFields,";"):
##arcpy.AddMessage(cF)
caseRowValue = caseRow.getValue(cF)
##arcpy.AddMessage(caseRowValue)
if caseRowValue == None:
pass #don't add a null value to the list!
else:
cFList = cFList + ";" + caseRowValue
cFList = cFList.strip(";")
##arcpy.AddMessage("list of data=" + cFList)
cList.append(cFList)
caseRow = caseRows.next()
del caseRow
del caseRows
caseVals = uniqify(cList)
caseVals.sort()
for c in caseVals:
#arcpy.AddMessage(c)
whereClause = ""
vaList = string.split(c,";")
for i in range(len(vaList)):
wf = caseList[i]
wv = vaList[i]
wherei = "\"" + wf + "\"='" + wv + "'"
whereClause = whereClause + " AND " + wherei
#arcpy.AddMessage(whereClause)
whereClause = whereClause.strip(" AND ")
arcpy.AddMessage(whereClause)
# Read the values from the Source Table
valueList = []
searchRows = arcpy.SearchCursor(inputTable,whereClause)
searchRow = searchRows.next()
while searchRow:
searchRowValue = searchRow.getValue(statField)
if searchRowValue == None:
pass #don't add a null value to the list!
else:
valueList.append(searchRowValue)
searchRow = searchRows.next()
del searchRow
del searchRows
valueList.sort()
listLength = len(valueList)
if listLength == 0:
print "Every value was null! Exiting script..."; sys.exit()
elif listLength % 2 == 0: #even, so get the mean of the 2 center values
##medianValue = (valueList[listLength / 2]) + valueList[listLength / 2 - 1]) / 2.0
medianValue = valueList[listLength / 2]
else: #odd, so it's easy!
medianValue = valueList[listLength / 2]
#Write the values to the output table
listCount = listLength
insertRows = arcpy.InsertCursor(outTable)
insertRow = insertRows.newRow()
for i in range(len(vaList)):
outCaseField = caseList[i]
outCaseVar = vaList[i]
insertRow.setValue(outCaseField, outCaseVar)
insertRow.setValue(outCountField, listCount)
insertRow.setValue(outMedField, medianValue)
insertRows.insertRow(insertRow)
del valueList
del listCount
del medianValue
