## SIMULTANEOUS CAUSAL INFERENCE

## FOR

## MULTIPLE TREATMENTS

VIA

## SUFFICIENCY

Yannis G. Yatracos

Cyprus U. of Technology

August 5, 2014

e-mail: yannis.yatracos@cut.ac.cy

## Summary

Some units from a population receive the same treatment that is different from treatments available for other reservoir populations. The minimal sufficient statistic s for the pre-treatment x -covariates's distributions in the populations is the coarsest balancing score. s is used to select matching units for simultaneous causal comparisons of multiple treatments.Necessary and sufficient conditions on the posterior distribution of the treatment variable (given x ) determine whether a statistic is either sufficient or minimal sufficient for the x-covariates' distributions. Results in the literature are thus extended. Strong ignorability of treatment assignment given s x ( ) is also established. Consequently, the expected treatments' differences given s x ( ) are shown to be simultaneously unbiased for the average causal effects of all treatments' differences. The existing statistical theory for s and its estimates support their use in causal inference.

Some key words: Causal Inference, Coarsest Balancing Score, Generalized Propensity Scores, Generalized Linear Models, Matching, Minimal Sufficient Statistic, Propensity Function

Running head: S-matching for Simultaneous Causal Inference

## 1 Introduction

When reservoir populations receive each a different treatment, the minimal sufficient statistic s of the pre-treatment x -covariates' distributions is used to select matching units for simultaneous causal comparisons of the treatments. Strong ignorability of treatment assignment given s x ( ) is established and the expected treatments' differences given any s -value are shown to be simultaneously unbiased for the average causal effects of all treatments' differences. Criteria are provided to obtain either s or a balancing score from the posterior distribution, q t x , ( | ) of the treatment variable, T, given the x -covariates. The results in Imai and van Dyk (2004) are extended, providing balancing scores for a larger class of q t x ( | )-models.

For two treatments, t = 1 , 2 , Rosenbaum and Rubin (1983) propose the scalar propensity score e x ( ) to balance the pre-treatment covariates, x ( ∈ R p ) , of the n units in the treatment groups; e x ( ) is the conditional probability of receiving, say, treatment 1 given x. It is stated therein that e x ( ) is the coarsest balancing score and it is showed, among others, that if treatment assignment and the potential units' responses to treatments, r i (1) and r i (2) , i = 1 , . . . , n, are conditionally independent given x, then the difference between the sample treatments' means given e x ( ) is unbiased for the average causal effect E r { (2) -r (1) } ; E denotes expectation over the whole population.

R. Bahadur recognized that e x ( ) is equivalent to the likelihood ratio of the x -populations' densities which is minimal sufficient (Rubin and Thomas, 1996, p. 250). As we found recently, the minimal sufficient statistic s is used with more than two treatments a) for dimension reduction theory, in particular when propensities do not exist (Nelson and Noorbaloochi, 2009), and b) in causal inference for each pair of treatments assuming that strong ignorability

of treatment assignment holds (Noorbaloochi, Nelson and Asgharian, 2010, p. 12, lines -2, -1, p. 13, lines 1-4). This assumption and the use of the x -covariates distributions to obtain s constitute two of the differences with this work.

For more than two treatment levels Joffe and Rosenbaum (1999) study causal effects using a small number of balancing linear functions of x. For multi-valued categorical treatments Imbens (2000) introduced the generalized propensity score P T ( = t x | ) and used it to estimate average causal effects for treatments' pairs but not for simultaneous causal comparisons of all treatments; see Imbens (2000, page 709, lines 15 to 17 and lines -4 to -1). For general treatment regimes with any type of t -values, Imai and Van Dyk (2004) introduce the propensity function e ψ ( ·| x , ) that has the form of q t x , ( | ) and assume that for all x -values it depends on x only through a unique, finite dimensional parameter θ ψ ( x ) ( Assumption 3 , p. 856); ψ is a known parameter. They show that e ψ ( ·| x ) (i.e. θ ψ ( x )) is balancing score ( Result 1 , p. 856), not necessarily the coarsest , and use it for simultaneous causal comparisons. Their Assumption 3 , that does not hold when q t x ( | ) belongs to a general exponential family, and θ ψ 's dimensionality constitute two of the differences with this work where s is not necessarily finite-dimensional and is obtained without additional assumptions. Corollary 3.1 shows also that Assumption 3 can be weakened to allow for more q t x ( | ) models and Result 1 holds automatically.

An estimate of s may not be sufficient and the same also holds for estimates of the propensity score, the generalized propensity score and the propensity function. This has been neglected so far in the Causal Inference literature that has not its own tools to confirm the balancing property unlike s ; see section 4.

The findings explain clearly what 'matching' means. Units from different populations receiving each a different treatment form a matching group

when they provide the same information for the x -covariates' distributions. Such groups can be used in causal comparisons, for example, to determine the 'right' dose for a new drug, by examining simultaneously the expected response differences E r t { ( 2 ) -r t ( 1 ) | s x ( ) = s 0 } , E r t { ( 3 ) -r t ( 2 ) | s x ( ) = s 0 } , . . . , E { ( r t ( k ) -r t ( k -1 ) | s x ( ) = s 0 } for different doses' levels t 1 &lt; t 2 &lt; . . . &lt; t k .

The framework is presented in section 2. The main theoretical results and s for generalized linear models are in section 3. In section 4 some directions are given for s -matching's implementation in practice. The proofs are in the Appendix.

## 2 Causal inference framework and assumptions

For a random vector U use p U ( u ) to denote its density (but also its probability). When random vector V is also available use p u v ( | ) to denote the conditional density of U given V. Let T denote the treatments and let T be the treatment variable with values t in T and prior density π T . Treatment t is used in selected units of population P t having balanced x -covariates with respect to T . The units in P t have covariates x ∈ C P ( t ) ⊂ R p and unless otherwise stated it is assumed that C P ( t ) = C , t ∈ T . Let p x t ( | ) denote the x -covariates' density of units in P t and let D T = { p x t , ( | ) t ∈ T } ; p X ( x ) is the marginal density of the x -covariates. The notation p x t ( | ) does not mean necessarily that p is the same density with the parameter t changing, t ∈ T , but simply denotes the covariates' distribution in P t . Use q t x ( | ) to denote T 's density (or probability) given the x -covariates. For unit i, r i ( ) is the response for treatment t t and the potential outcomes R is the set { r i ( ) t , t ∈ T , for i = 1 , . . . , n } . Conditional independence of x and y given z is denoted by x ⊥ ⊥ y z | (Dawid,1979). The expression 'covariates u, v match' means that the units with these covariates

match.

Assumption 1 (Stable Unit Treatment Value Assumption ( SUTVA ), Rubin, 1980, 1990) The distribution of potential outcomes for one unit is assumed to be independent of potential treatment status of another unit given the observed covariates.

Assumption 2 (Strong ignorability of treatment assignment given x, Rosenbaum and Rubin, 1983)

- (i) R and T are conditionally independent given x : R⊥ ⊥ T | x, and
- (ii) for every t ∈ T , 0 &lt; p t x ( | ) (or equivalently 0 &lt; p x t ( | )).

Recall that b x ( ) is a balancing score if the conditional distribution of x given b x ( ) is the same for all treatment values, i.e.

$$p ( x | t, b ( x ) ) = p ( x | b ( x ) ), \text{ for all } t \in \mathcal { T }.$$

From (1), thinking of t as parameter value for the distribution of x it follows that b x ( ) is a sufficient statistic for the family D T = { p x t ( | ); t ∈ T } .

## 3 Matching and Causal Inference with s

In this section the minimal sufficient statistic s x ( ) is assumed known. This is possible for various models. The results are also applicable for large samples when s is estimated.

## 3.1 Vector valued s

The first result, obtained directly from statistical theory, extends Nelson and Noorbaloochi (2009, p. 619, Theorem 1) and justifies the use of s for countably finite or countably infinite treatments. It also indicates that for multiple treatments s will be often vector valued.

Proposition 3.1 Assume that T consists of either countably finite or countably infinite treatments and that the covariates' distributions D T have all common support. Let T k = { t 1 , . . . , t k } ⊂ T and let

$$s ( x ) = s ^ { ( 1 ) } ( x ) ^ { 1 } = \left ( \frac { p ( x | t _ { 2 } ) } { p ( x | t _ { 1 } ) }, \frac { p ( x | t _ { 3 } ) } { p ( x | t _ { 1 } ) }, \dots, \frac { p ( x | t _ { k } ) } { p ( x | t _ { 1 } ) } \right )$$

When s is sufficient for T (i.e. for D T ) it is also minimal sufficient.

Remark 3.1 When D T k = { p x t , t ( | ) ∈ T } k do not have common support, s x ( )'s dimension depends on the x -values (Lehmann and Casella, 1998, p. 70, Theorem 9.1).

The s -Matching Rule for covariates: Match u to v when s u ( ) = s v . ( )

The s -matching rule can be used for any number of units using their covariates.

The next proposition shows that s -matching is not changed when

/negationslash

$$s ^ { ( j ) } ( x ) = \left ( \frac { p ( x | t _ { 1 } ) } { p ( x | t _ { j } ) }, \dots, \frac { p ( x | t _ { j - 1 } ) } { p ( x | t _ { j } ) }, \frac { p ( x | t _ { j + 1 } ) } { p ( x | t _ { j } ) }, \dots, \frac { p ( x | t _ { k } ) } { p ( x | t _ { j } ) } \right ), j \neq 1, \quad ( 3 )$$

is used instead of s = s (1) in (2).

Proposition 3.2 If s u ( ) = s v , ( ) then s ( j ) ( u ) = s ( j ) ( v , ) j &gt; 1 .

Without loss of generality s x ( ) = s (1) ( x ) is used in this section.

Propositions 3.1 and Remark 3.1 indicate clearly that with several treatments s is not expected to be scalar. Nelson and Noorbaloochi (2009, p. 619) point out that s x ( ) may be infinite dimensional thus contradicting the existence of a finite dimensional parameter θ in Assumption 3 (Imai and Van Dyk, 2004). In Noorbaloochi, Nelson and Asgharian (2010, p. 8, lines 1, 2) it is also

1 In s (1) ( x , ) (1) indicates the denominator is p x t ( | 1 ) .

mentioned that ' In general, there is no univariate propensity score .' This is indirectly confirmed for several data sets with high dimensional x -covariates for which the scalar propensity score deteriorates more as x 's dimension increases (King et al., 2011, p. 18).

We revisit an example in Rosenbaum and Rubin (1983, p. 47) when the number of treatments k is larger than 2.

Example 3.1 Let p x t ( | ) be a polynomial exponential family distribution,

$$p ( x | t ) = h ( x ) \exp \{ P _ { t } ( x ) \}, \ t = 1, \dots, k,$$

with P t ( x ) a degree m polynomial. Then, the statistic

$$\left ( \ln \frac { p ( x | t _ { 2 } ) } { p ( x | t _ { 1 } ) }, \dots, \ln \frac { p ( x | t _ { k } ) } { p ( x | t _ { 1 } ) } \right ) = \left ( P _ { 2 } ( x ) - P _ { 1 } ( x ), \dots, P _ { k } ( x ) - P _ { 1 } ( x ) \right ) \\ = \left ( Q _ { 1 } ( x ), \dots, Q _ { k - 1 } ( x ) \right )$$

is equivalent to the minimal sufficient statistic (2) with Q x i ( ) a degree m polynomial, i = 1 , . . . , k -1 .

## 3.2 Causal inference framework and s

To obtain s using the likelihood ratios in (2) the densities of the covariates in D T = { p x t , ( | ) t ∈ T } have either to be known or to be estimated. This may not be possible in practice. Results are now presented to determine s with the causal inference framework and without using D T , simply from the conditional density q t x ( | ) of T given the x -covariates. The first result involves D T but it is used to prove subsequent results.

Proposition 3.3 (see, e.g. Chen, 2010, Ch. 6) Let D T be the family of the x -covariates densities. Assume that there exist function s ∗ ( x ) such that for

any covariates x 1 and x 2 the ratio p x ( 1 | t ) p x ( 2 | t ) is constant as function of t if and only if s ∗ ( x 1 ) = s ∗ ( x 2 ) . Then, s ∗ ( x ) is minimal sufficient.

The tool to determine s via q t x ( | ) is the decomposition

$$p ( x | t ) = q ( t | x ) \cdot p _ { X } ( x ) \cdot \pi _ { T } ^ { - 1 } ( t ),$$

that leads to T 's posterior factorization criterion (PFC) and the coarsest balancing score criterion (CBSC).

Proposition 3.4 Let X and T be random vectors in Euclidean spaces with densities, respectively, p X and p T and with conditional densities p x t ( | ) and q t x . ( | ) Then, a) (Posterior Factorization Criterion) s x ( ) is sufficient statistic for D T if and only if

$$q ( t | x ) = g _ { 1 } ( s ( x ), t ) \cdot g _ { 2 } ( t ) \cdot g _ { 3 } ( x ) \ \forall \ x, \ t.$$

b) (Coarsest Balancing Score Criterion) Assume that for any x 1 and x 2 the ratio q t x ( | 1 ) q t x ( | 2 ) is independent of t if and only if s ∗ ( x 1 ) = s ∗ ( x 2 ) . Then, s ∗ is minimal sufficient statistic for D T .

Proposition 3.4 is used to derive directly previous results in the literature.

Example 3.2 (The propensity score, Rosenbaum and Rubin, 1983) The treatments T = { 1 2 , } and the propensity score e x ( ) = q (1 | x . ) For q t x ( | ) it holds

$$q ( 1 | x ) = e ( x ), \ q ( 2 | x ) = 1 - e ( x ).$$

From (6) and Proposition 3.4 a ) s x ( ) = e x ( ) is sufficient statistic. Since the ratio q t x ( | 1 ) q t x ( | 2 ) is independent of t for all x , 1 x 2 if and only if

$$\frac { q ( 1 | x _ { 1 } ) } { q ( 1 | x _ { 2 } ) } = \frac { q ( 2 | x _ { 1 } ) } { q ( 2 | x _ { 2 } ) } \quad \Leftrightarrow \quad e ( x _ { 1 } ) = e ( x _ { 2 } ),$$

from Proposition 3.4 b ) e x ( ) is minimal sufficient. The same result is obtained via Proposition 3.1.

Example 3.3 (The propensity function, Imai and van Dyk, 2004) For various kinds of treatments t ∈ T , the propensity function e ψ ( ·| x ) = q ψ ( ·| x ) depends on x only through the unique, finite dimensional parameter θ ψ ( x )( Assumption 3); ψ is known parameter. Therefore, q ψ ( t x | ) has form (5) with g 3 ( x ) = 1 for every x and from Proposition 3.4 a ) s x ( ) = θ ψ ( x ) is sufficient.

Proposition 3.4 extends the results in Imai and Van Dyk (2004) by weakening their Assumption 3 to accommodate more q t x ( | )-models.

Corollary 3.1 Let q ψ ( t x | ) be the density of T given the x -covariates; ψ are known parameters. Assume that there are functions θ 1 ,ψ ( x , ) g 1 ,ψ , g 2 ,ψ , g 3 ,ψ such that

$$q _ { \psi } ( t | x ) = g _ { 1, \psi } ( \theta _ { 1, \psi } ( x ), t ) \cdot g _ { 2, \psi } ( t ) \cdot g _ { 3, \psi } ( x ) \ \forall \ x, \ t.$$

Then, a) θ 1 ,ψ ( x ) is a balancing score, and

b) θ 1 ,ψ ( x ) is the coarsest balancing score when for every x , 1 x , 2 the ratio

$$\frac { g _ { 1, \psi } ( \theta _ { 1, \psi } ( x _ { 1 } ), t ) } { g _ { 1, \psi } ( \theta _ { 1, \psi } ( x _ { 2 } ), t ) }$$

is independent of t if and only if θ 1 ,ψ ( x 1 ) = θ 1 ,ψ ( x 2 ) .

The minimal sufficient statistic s for DT is now determined when T 's posterior is a generalized linear model.

Example 3.4 Assume the treatment variable T with values in R d is modeled given the x -covariates ( ∈ R p ) with a generalized linear model in canonical form, i.e.

$$q ( t | x ) = \exp \{ t ^ { \prime } b ( x ) + c ( x ) + d ( t ) \} ;$$

t ′ denotes t 's transpose and n ∈ R . d From Proposition 3.4 a), b x ( ) is sufficient statistic for D T . Since the ratio

$$\frac { q ( t | x _ { 1 } ) } { q ( t | x _ { 2 } ) } = e x p \{ t ^ { \prime } [ b ( x _ { 1 } ) - b ( x _ { 2 } ) ] + c ( x _ { 1 } ) - c ( x _ { 2 } ) \}$$

is independent of t for all t, x 1 , x 2 if and only if b x ( 1 ) = b x ( 2 ) , from Proposition 3.4 b) b x ( ) is also minimal sufficient.

There are different forms b x ( ) can have. For example, when b takes real values,

$$b ( x ) = \phi ( \sum _ { j = 0 } ^ { p } \beta _ { j } x _ { j } ), \quad b ( x ) = \sum _ { j = 0 } ^ { K } \beta _ { j } \phi _ { j } ( x ) ;$$

x 0 = 1 , x j is x 's j -th coordinate, β j ∈ R, the functions φ, φ j are assumed to be known and real valued, j 's values are according to the corresponding sum.

When b x ( ) is known it can be used for matching units from different populations. When b x ( ) is not known, it has to be estimated with ˆ ( b x ) that is used for matching.

Remark 3.2 For the generalized linear model (8) Assumption 3 in Imai and Van Dyk (2004) does not hold because q t x ( | ) depends on x via c x ( ) also.

## 3.3 Simultaneous causal comparisons

The key result allowing for simultaneous causal comparison of several treatments follows, establishing strong ignorability of treatment assigment given s x . ( )

Proposition 3.5 Under Assumption 2, for the responses R and the treatment variable T = t it holds

$$p \{ t, \mathcal { R } | s ( x ) = s \} = p \{ t | s ( x ) = s \} \cdot p \{ \mathcal { R } | s ( x ) = s \}. \quad \quad ( 1 0 )$$

Proposition 3.5 suggests simultaneous causal comparisons using s x ( ) to balance subpopulations for all treatments and obtain unbiased estimates of the average treatment effects.

/negationslash

Proposition 3.6 Suppose that treatment assignment is strongly ignorable (Assumption 2) and that a value s 0 of s x ( ) is randomly sampled from the population of units with covariates x ∈ C . Units receiving treatments t i and t j are sampled with s -value for their covariates equal s , 0 i = j. Then, the expected difference in response for the units chosen is the expected treatment effect at s x ( ) = s . 0 The mean of such pair differences over all s x ( ) -values is unbiased for the average treatment effect E r t { ( i ) -r t ( j ) } and the same holds, concurrently given s x , ( ) for any number of average treatment effects.

## 4 Implementation

For the s -matching's implementation there are practical issues some of which depend on the data to be analyzed and the assumptions on the x -covariates models. Among these issues a) s should be determined, b) when likelihood ratios have to be estimated, the curse of dimensionality problem should be addressed, and c) the dimensionality of s may be reduced if there is no much loss of information. Some directions for a) c) -follow.

Known theorems in statistics (see, e.g., Lehmann and Casella, 1998) allow to obtain the minimal sufficient statistic s. Proposition 3.1 is used with a small number k of treatments to define a minimal sufficient statistic s ∗ for the corresponding distributions, D T k , and then show that s ∗ is sufficient for all the x -covariates distributions, D T , i.e. s = s . ∗ For example, if p x t ( | ) follows a normal distribution with mean t and known variance (say) 1, the minimal

/negationslash sufficient statistic, s , ∗ is determined for the distributions D { t 1 ,t 2 } ; t 1 , t 2 are treatments, t 1 = t 2 . Neyman's factorization criterion shows that s ∗ is sufficient for all t -values, so it is minimal sufficient for D T . Alternatively, when p x t ( | ) belongs to a p -parameter exponential family in canonical form, Neyman's factorization criterion determines s x ( ) that is also minimal sufficient if the parameter space T contains an open, p -dimensional rectangle.

For the implementation of Proposition 3.1 in applications only a subpopulation ˜ P t of P t may be available. Let MP t denote the units to be matched from ˜ -subpopulation. P t Use s = s (1) in (2) to match a unit in MP t having covariates u with a unit from ˜ P r having covariates v m,r ∈ C P ( ˜ r ) , such that

$$v _ { m, r } = \arg \min _ { v \in \mathcal { C } ( \tilde { \mathcal { P } } _ { r } ) } | | s ( u ) - s ( v ) | | ^ { 2 }, \ r \in T - \{ t \} ; \quad \quad ( 1 1 )$$

|| · || is the usual Euclidean distance in R p and in v m,r the index m denotes 'matching' unit from sub-population ˜ P r . This approach is the nearest neighbor 1 : 1 matching with replacement and can be properly modified for 1 : k matching with or without replacement. For more information on matching methods and for optimal matching questions see, e.g., Rosenbaum (1989) and Stuart (2010).

Additional matching sets for MP t can be obtained using s = s ( j ) (or its estimates) in (11), j = 2 , . . . , k, and the decision maker can select the 'best' matching set, for example, that with the nearest means to the MP t covariates' means with respect to || · || or the sup-norm distance || · || ∞ .

When the form of the x -covariates densities in D T is not known and T = T k , the usual nonparametric estimation of each density in the ratios (2) is affected by the curse of dimensionality of the x 's. Rather than estimating separately each density one may use the approach adopted in Machine Learning for determining the ratio of the densities from the training and test data. The densities'

ratio is expressed as linear model with respect to a basis of functions. The coefficients are estimated using observations from the two populations according to a given method that usually ends with a convex minimization problem. For the description of the estimation methods in Machine Learning see, e.g., Sugiyama et al. (2007) and Nguyen et al. (2010). For independent x -samples from nonparametric models, one for each t ∈ T , the class of empirical distributions is minimal sufficient.

The use of s is supported by the existing statistical theory that allows i) to obtain approximate sufficient statistics and evaluate the approximation's error (see, e.g., Le Cam, 1964, Joyce and Marjoram, 2008), and ii) use a principal components transformation of its components to investigate whether a sufficient summary of smaller dimension exists (Nelson and Noorbaloochi, 2009).

## 5 Appendix

Proof of Proposition 3.1: It is direct consequence of Theorem 6.12, in Lehmann and Casella, 1998, p. 37 and related theorems therein.

Proof of Proposition 3.2: Since s u ( ) = s v , ( ) it holds

$$\frac { p ( u | t _ { i } ) } { p ( u | t _ { 1 } ) } = \frac { p ( v | t _ { i } ) } { p ( v | t _ { 1 } ) }, \ i = 2, \dots, k.$$

/negationslash

In (12), divide the i -th equality with the j -th equality, i = j, and invert the j -th equality to obtain

/negationslash

$$\frac { p ( u | t _ { i } ) } { p ( u | t _ { j } ) } = \frac { p ( v | t _ { i } ) } { p ( v | t _ { j } ) }, \, i \neq j, \, \text{or} \ s ^ { ( j ) } ( u ) = s ^ { ( j ) } ( v ).$$

Proof of Proposition 3.4: a) From Neyman's Factorization criterion is sufficient statistic if and only if

## s

$$p ( x | t ) = h _ { 1 } ( s ( x ), t ) \cdot h _ { 2 } ( x ) \ \forall \ x, \ t,$$

and from decomposition (4)

$$\Leftrightarrow \quad q ( t | x ) \cdot p _ { X } ( x ) \cdot \pi _ { T } ^ { - 1 } ( t ) = h _ { 1 } ( s ( x ), t ) \cdot h _ { 2 } ( x )$$

$$\Leftrightarrow \quad q ( t | x ) = h _ { 1 } ( s ( x ), t ) \cdot \pi _ { T } ( t ) \cdot h _ { 2 } ( x ) \cdot p _ { X } ^ { - 1 } ( x ).$$

Equality (5) follows with

$$g _ { 1 } ( s ( x ), t ) = h _ { 1 } ( s ( x ), t ), \quad g _ { 2 } ( t ) = \pi _ { T } ( t ), \quad g _ { 3 } ( x ) = h _ { 2 } ( x ) \cdot p _ { X } ^ { - 1 } ( x ).$$

Conversely, from (5) Neyman's Factorization criterion (13) is obtained via (4).

- b) From decomposition (4) the ratio

$$\frac { p ( x _ { 1 } | t ) } { p ( x _ { 2 } | t ) } = \frac { q ( t | x _ { 1 } ) } { q ( t | x _ { 2 } ) } \cdot \frac { p _ { X } ( x _ { 1 } ) } { p _ { X } ( x _ { 2 } ) }$$

is independent from t if and only if the ratio

$$\frac { q ( t | x _ { 1 } ) } { q ( t | x _ { 2 } ) }$$

is independent of t and this holds if and only if

$$s ^ { * } ( x _ { 1 } ) = s ^ { * } ( x _ { 2 } ).$$

Thus, from Proposition 4.1 s ∗ is minimal sufficient statistic.

Proof of Corollary 3.1: Both parts follow from Proposition 3.4.

Proof of Proposition 3.5: The proof follows the lines in Imai and Van Dyk (2004),

$$p \{ x, t, \mathcal { R } | s ( x ) = s \} = p \{ x, t | s ( x ) = s \} \cdot p \{ \mathcal { R } | x, t, s ( x ) = s \}$$

$$= p \{ t | s ( x ) = s \} \cdot p \{ x | t, s ( x ) = s \} \cdot p \{ \mathcal { R } | x, t, s ( x ) = s \} \\ = p \{ t | s ( x ) = s \} \cdot p \{ x | s ( x ) = s \} \cdot p \{ \mathcal { R } | x, s ( x ) = s \}.$$

The third equality is obtained using Proposition 3.1 and strong ignorability of treatment assignment given x (Assumption 2). It follows that

$$p \{ t, x, \mathcal { R } | s ( x ) = s \} = p \{ t | s ( x ) = s \} \cdot p \{ x, \mathcal { R } | s ( x ) = s \}$$

Integrating both sides of the last equation over the x 's for which s x ( ) = s, we obtain that given s x ( ) = s, R and T are independent.

Proof of Proposition 3.6: From Assumption 2,

$$E \{ r ( t _ { i } ) | s ( x ) = s, T = t _ { i } \} - E \{ r ( t _ { j } ) | s ( x ) = s, T = t _ { j } \}$$

$$= E \{ r ( t _ { i } ) | s ( x ) = s \} - E \{ r ( t _ { j } ) | s ( x ) = s \} = E \{ r ( t _ { i } ) - r ( t _ { j } ) | s ( x ) = s \}$$

and it follows that

$$E _ { s } \left [ E \{ r ( t _ { i } ) - r ( t _ { j } ) | s ( x ) = s \} \right ] = E \{ r ( t _ { i } ) - r ( t _ { j } ) \} ;$$

E s denotes expectation with respect to all values s of s x , ( ) x ∈ C .

## References

- [1] Agresti, A. (2002) Categorical Data Analysis. Wiley.
- [2] Chen, H. (2010) Notes in Statistical Inference, Chapter 6, http://www.math.ntu.edu.tw/ ∼ hchen/teaching/StatInference/notes/ch6.pdf
- [3] Dawid, A. P. (1979) Conditional independence in statistical theory (with discussion). JRSS B 41 , 1-31.

- [4] Imai, K. and Van Dyk, D. A. (2004) Causal inference with general treatment regimes: generalizing the propensity score. JASA 99 , 854866.
- [5] Imbens, G. W. (2000) The role of the propensity score in estimating dose-response functions. Biometrika 87 , 706-710.
- [6] Joffe, M. M. and Rosenbaum, P. R. (1999) Propensity scores. American Journal of Epidemiology 150 , 327-333.
- [7] Joyce, P. and Marjoram, P. (2008) Approximately sufficient statistics and Bayesian computation. Statistical applications in Genetics and Molecular Biology 7 : Iss. 1, Article 26.
- [8] King, G., Nielsen, R., Coberley, C. Pope, J. E. and Wells, A. (2011) Comparative effectiveness of matching methods for causal inference. At gking.harvard.edu/publications .
- [9] Le Cam, L. M. (1964) Sufficiency and Approximate Sufficiency. Ann. Math. Stat. 35 , 1419-1455.
- [10] Lehmann, E. L. and Casella, G. (1998) Theory of Point Estimation. Springer, New York.
- [11] McCullagh, P. and Nelder, J. A. (1989) Generalized Linear Models. Chapman and Hall, London.
- [12] Nelson, D. and Noorbaloochi, S. (2009) Dimension reduction summaries for balanced contrasts. J. Stat. Plan. and Inf. 139 , 617-628.

- [13] Nelson, D. and Noorbaloochi, S. (2013) Information preserving sufficient summaries for dimension reduction. J. Multiv. Anal. 115 , 347358.
- [14] Nguyen, X.-L., Wainwright, M. J. and Jordan, M. I. (2010) Estimating divergence functionals and the likelihood ratio by convex minimization. IEEE Trans. Inform. Theory 56 , 5847-5861.
- [15] Noorbaloochi, S., Nelson, D. and Asgharian, M. (2010) Balancing and elimination of nuisance variables. Intern. J. of Biostat. 6 , 2 , Article 6 .
- [16] Rosenbaum, P. R. and Rubin, D. B. (1983) The central role of the propensity score in observational studies for causal effects. Biometrika 70 , 41-55.
- [17] Rosenbaum, P. R. (1989) Optimal Matching for Observational Studies. JASA 84 , 1024-1032.
- [18] Rubin, D. B. and Thomas, N. (1996) Matching using estimated propensity scores: relating theory to practice. Biometrics 52 , 249264.
- [19] Rubin, D. B. (1980) Discussion of paper by D. Basu. JASA 75 , 591593.
- [20] Rubin, D. B. (1990) Comments on 'On the Application of Probability Theory to Agricultural Experiments. Essay on Principles. Section 9,' by J. Neyman, Stat. Science 5 , 472-480.
- [21] Stuart, E. A. (2010) Matching methods for causal inference: A review and a look forward. Stat. Sci. 25 , 1-21.

[22] Sugiyama, M., Nakajima, S, Kashima, H., Von B¨nau, P. , Kawanabe, u M. (2007) Direct importance estimation with model selection and its application to covariate shift adaptation. NIPS
Yannis G. Yatracos
Simultaneous causal inference for multiple treatments via sufficiency
Some units from a population receive the same treatment that is different
from treatments available for other reservoir populations. The minimal
sufficient statistic $s$ for the pre-treatment $x$-covariates's distributions
in the populations is the coarsest balancing score. $s$ is used to select
matching units for simultaneous causal comparisons of multiple
treatments.Necessary and sufficient conditions on the posterior distribution of
the treatment variable (given $x$) determine whether a statistic is either
sufficient or minimal sufficient for the x-covariates' distributions. Results
in the literature are thus extended. Strong ignorability of treatment
assignment given $s(x)$ is also established. Consequently, the expected
treatments' differences given $s(x)$ are shown to be simultaneously unbiased
for the average causal effects of all treatments' differences. The existing
statistical theory for $s$ and its estimates support their use in causal
inference.