>>145
へー
> x <- rnorm(25)
> sum(abs(x - mean(x)))/length(x)
[1] 0.7459197
> sqrt(sum((x - mean(x))^2)/length(x))
[1] 0.9833284
確かにσは少し大きいね。勉強になった。
計算コストはσの方が圧倒的に有利だけどね。
> system.time(for(i in 1:1000){sum(abs(x - mean(x)))/length(x)})
ユーザ システム 経過
0.013 0.000 0.014
> system.time(for(i in 1:1000){sqrt(sum((x - mean(x))^2)/length(x))})
ユーザ システム 経過
0.004 0.000 0.003
このあたりが歴史的にσが好まれた理由だと思っている。