Hi Lucas,
The easiest way to do this is the Cross Validation geoprocessing tool. It takes a geostatistical layer as input, and it creates a derived output with all of the summary statistics as properties.
Quick Python code example:
outCV = arcpy.ga.CrossValidation('myGALayer')
outCV.rootMeanSquare
'294.63006397102004'