Reply to thread

Message: [QUOTE="Paul Sanchez, post: 4464561"] Don't know if this will meet your needs, it's not fancy but it works for me... class SimpleStats attr_reader :n, :sampleMean, :min, :max def initialize reset end def reset @n = 0 @min = 1.0 / 0.0 @max = -@min @sumsq = @sampleMean = @min / @min end def newObs(datum) x = datum.to_f @max = x if !@max || x > @max @min = x if !@min || x < @min @n += 1 if (@n > 1) delta = x - sampleMean @sampleMean += delta / @n @sumsq += delta * (x - @sampleMean) else @sampleMean = x @sumsq = 0.0 end end def sampleVariance @sumsq / (@n - 1) end def stdDeviation Math::sqrt sampleVariance end end [/QUOTE]