Calculate a weighted distance between two vectors, using pairwise complete observations.

Simple function to count the fraction of different elements (in the same position) between two vectors of the same length, after removing elements from both vectors corresponding to positions that are NA in either vector. Distance is the probability for observing the matches and mismatches in two random patients.

Usage

vectorWeightedDist(x, y, k, l)

Arguments

x: a matrix
y: a matrix with the same number of columns as x
k: row in x to test for differences
l: row in y to test for differences

Value

Returns a numeric value, the log of the probability of observing the matches in x and y

Author

Levi Waldron, Markus Riester, Marcel Ramos

Examples


mymat1 <- matrix(rnorm(20), ncol = 5)
mymat1[1, 4] <- NA
mymat2 <- matrix(rnorm(20), ncol = 5)
vectorWeightedDist(mymat1, mymat2, 1, 2)
#> [1] 1