Download this file
1 2 3 4 5 6 7 8 9
w.sd = function (x, w) { sumw = sum(w) sumw.sq = sum(w)^2 w.sq = sum(w^2) denom = sum(w * ((x - mean(x))^2)) sqrt((sumw * denom) / (sumw.sq - w.sq)) }