On Kelvin Lawrence's excellent gremlin guide, there's a section on how to compute the standard deviation for a list of property values.
mean=g.V().hasLabel('airport').values('runways').mean().next()
count = g.V().hasLabel('airport').count().next()
g.withSideEffect("m",mean).
withSideEffect("c",count).
V().hasLabel('airport').values('runways').
math('(_ - m)^2').sum().math('_ / c').math('sqrt(_)')
How would you condense this down into a single query?