Weighted variance function

PHOTO EMBED

Wed Sep 01 2021 00:49:08 GMT+0000 (Coordinated Universal Time)

Saved by [deleted user] #r

wt_var <- function(x, w, na.rm = FALSE) {
  if (na.rm) {
    na <- is.na(x) | is.na(w)
    x <- x[!na]
    w <- w[!na]
  }
  wm <- weighted.mean(x, w)
  sqrdev <- (x - wm)^2
  (sum(sqrdev*w))
}
content_copyCOPY

Calculates variance of a weighted discrete distribution