harmonize: Harmonisation by calibration

In RJauslin/transportMatching: Statistical Matching using optimal transport and balanced sampling.

Description

Harmonisation by calibration

Usage

harmonize(X1, d1, id1, X2, d2, id2, totals)

Arguments

X1	dataset 1
d1	weights 1
id1	identifiers 1
X2	dataset 2
d2	weigths 2
id2	identifiers 2
totals

Examples

N <- 10000
X <- data.frame(x1 = rnorm(N,0,1), x2 = rnorm(N,0,1))

n1=1000
n2=3000
s1=srswor(n1,N)
s2=srswor(n2,N)
id1=(1:N)[s1==1]
id2=(1:N)[s2==1]
d1=rep(N/n1,n1)
d2=rep(N/n2,n2)
X1 = X[s1==1,]
X2 = X[s2==1,]

W <- harmonize(X1,d1,id1,X2,d2,id2)

colSums(W$w1*X1)
colSums(W$w2*X2)


############## if knowing the true totals

totals <- c(N,colSums(X))
W <- harmonize(X1,d1,id1,X2,d2,id2,totals)

colSums(W$w1*X1)
colSums(W$w2*X2)
colSums(X)

RJauslin/transportMatching documentation built on Dec. 18, 2021, 8:49 a.m.