This functions works quite equally as the wtd.var
function
from the Hmisc
package.
Usage
wtdVar(x , weights , na.rm = FALSE)
Arguments
- x
numeric vector
- weights
a numeric vector of non-negative weights
- na.rm
set to FALSE
to suppress checking for NAs. If TRUE
, NAs are removed
from x
as well as from weights
prior to variance estimation.
Examples
x <- c(50, 1, 25)
w <- c(1, 4, 1)
wtdVar(x, w)
#> [1] 417.7667