平均と標準偏差のスウィフト加速

Question

私は、スウィフトの平均値と標準偏差を計算するAccelerateを検討しています。

私は平均を行うことができます。標準偏差はどのようにすればよいですか？

let rr: [Double] = [ 18.0, 21.0, 41.0, 42.0, 48.0, 50.0, 55.0, 90.0 ] 

var mn: Double = 0.0 

vDSP_meanvD(rr, 1, &mn, vDSP_Length(rr.count)) 

print(mn) // prints correct mean as 45.6250 

// Standard Deviation should be 22.3155 

Answer 1

あなたは平均値と 二乗平均値（https://en.wikipedia.org/wiki/Standard_deviation#Identities_and_mathematical_propertiesとhttps://en.wikipedia.org/wiki/Algebraic_formula_for_the_varianceを比較）から標準偏差を計算することができる：あるいは（IOSの9.0以降またはMacOSの10.11以降）

import Accelerate 

let rr: [Double] = [ 18.0, 21.0, 41.0, 42.0, 48.0, 50.0, 55.0, 90.0 ] 

var mn: Double = 0.0 // mean value 
vDSP_meanvD(rr, 1, &mn, vDSP_Length(rr.count)) 

var ms: Double = 0.0 // mean square value 
vDSP_measqvD(rr, 1, &ms, vDSP_Length(rr.count)) 

let sddev = sqrt(ms - mn * mn) * sqrt(Double(rr.count)/Double(rr.count - 1)) 

print(mn, sddev) 
// 45.625 22.315513501982 

、vDSP_normalizeDを使用します。

var mn = 0.0 
var sddev = 0.0 
vDSP_normalizeD(rr, 1, nil, 1, &mn, &sddev, vDSP_Length(rr.count)) 
sddev *= sqrt(Double(rr.count)/Double(rr.count - 1)) 

print(mn, sddev) 
// 45.625 22.315513501982