¿Cómo se calculan las ponderaciones en la función matchit() de R para la coincidencia completa?

Question

En el ejemplo que tengo, un pequeño número de sujetos tratados se empareja con un gran número de controles no tratados. He utilizado

match= matchit(treat ~ SEX + ageint + bAVAL1 + bAVAL2, method = "full", caliper=0.3,
               exact  = c("SEX", "ageint"), data = data)

Tiene sentido que los sujetos tratados tengan peso=1. ¿Por qué los múltiples sujetos de control emparejados con cada sujeto tratado tienen pesos>1?

     weights          subclass
1    1.0000000        1
2    2.0421053        1
3    2.0421053        1
4    2.0421053        1
5    2.0421053        1
6    2.0421053        1
7    1.0000000        2
8    2.5526316        2
9    2.5526316        2
10   2.5526316        2
11   2.5526316        2
12   1.0000000        3
13   2.5526316        3
14   2.5526316        3
15   2.5526316        3
16   2.5526316        3

Answer 1

Esto se explica en el documentación para matchit() .

[Cada unidad se asigna a una subclase, que representa el par que que forma parte (en el caso de la correspondencia k:1) o el estrato al que pertenece (en el caso de la concordancia exacta, la concordancia exacta ampliada, la concordancia completa o la subclasificación). o subclasificación). La fórmula para calcular las ponderaciones depende del argumento suministrado a estimand. Un nuevo estrato "puntuación de propensión" ( p ) se calcula como la proporción de unidades en cada estrato que está en el grupo tratado, y a todas las unidades de ese estrato se les asigna esa puntuación de propensión. Las ponderaciones se calculan entonces utilizando las fórmulas estándar para las ponderaciones de probabilidad inversa: para el ATT, las ponderaciones son 1 para las unidades tratadas y p/(1-p) para el control unidades; para el ATC, los pesos son (1-p)/p para las unidades tratadas y 1 para las unidades de control; para la ATE, las ponderaciones son 1/p para las unidades unidades tratadas y 1/(1-p) para las unidades de control.

...

En cada grupo de tratamiento, los pesos se dividen por la media del de las ponderaciones no nulas de ese grupo de tratamiento para que las ponderaciones sumen el número de unidades de ese grupo de tratamiento. número de unidades en ese grupo de tratamiento.