0% found this document useful (0 votes)
9 views6 pages

Understanding Sufficient Statistics

A sufficient statistic is a statistic that contains all the information needed for estimating a parameter, θ, from a random sample. It is defined such that the conditional distribution of the sample given the statistic does not depend on θ. The document also discusses jointly sufficient statistics and provides the factorization criterion for identifying sufficient statistics.

Uploaded by

Elelan
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
9 views6 pages

Understanding Sufficient Statistics

A sufficient statistic is a statistic that contains all the information needed for estimating a parameter, θ, from a random sample. It is defined such that the conditional distribution of the sample given the statistic does not depend on θ. The document also discusses jointly sufficient statistics and provides the factorization criterion for identifying sufficient statistics.

Uploaded by

Elelan
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

5.

Sufficiency
A sufficient statistic with respect to 𝜃 is a statistic 𝑇 = 𝑡(𝑋1 , 𝑋2 , … , 𝑋𝑛 ) that contains all the
information that is useful for the estimation of 𝜃.
Let 𝑋1 , … , 𝑋𝑛 be a random sample from some density, say 𝑓 (∙ ; 𝜃). We defined a statistic 𝑇 =
𝑡(𝑋1 , … , 𝑋𝑛 ) to be a function of the sample which is also a random variable. It condenses the 𝑛
random variables 𝑋1 , … , 𝑋𝑛 into a single random variable. It is interesting to see that if we lost any
“information” by this condensing process. The condensing can also be viewed another way.
Let 𝖝 denote the range of values that (𝑋1 , … , 𝑋𝑛 ) can assume.
For example,

1. Bernoulli distribution, the 𝖝 is a collection of vectors with components in the vector is


either 0 or 1.
2. Normal distribution, the 𝖝 is n-dimensional Euclidean space.

Now a statistic induces or defines a partition of 𝖝.

Partition of 𝖝 is a collection of mutually disjoint subsets of 𝖝 whose union is 𝖝.


For example if a sample of size 3, say, (𝑋1 , 𝑋2 , 𝑋3 ) is selected from Bernoulli distribution, then:
𝖝 = {(0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), (1, 0, 1), (1, 1, 0), (0, 1, 1), (1, 1, 1)}
Define a statistic 𝑇 = 𝑡(𝑋1 , 𝑋2 , 𝑋3 ) = 𝑋1 + 𝑋2 + 𝑋3 .
That is [𝑡(𝑥1 , 𝑥2 , 𝑥3 ) = 𝑥1 + 𝑥2 + 𝑥3 ]

Then 𝑡(∙ , ⋯ , ∙) takes on the values 0, 1, 2, and 3. The partition of 𝖝 induced by 𝑡(∙ , ⋯ ,∙) consists
of the four subsets.
{0, 0, 0},
{(1, 0, 0), (0, 1, 0), (0, 0, 1)},
{(1, 0, 1), (1, 1, 0), (0, 1, 1)},
{(1, 1, 1)}
Corresponding, respectively to the four values 0, 1, 2, and 3 of 𝑡(∙ , ⋯ , ∙).

Then the statistic is really a condensation of 𝖝.


In the above example, if we use the statistic 𝑇 = 𝑡(𝑋1 , 𝑋2 , 𝑋3 ) = 𝑋1 + 𝑋2 + 𝑋3 , we have only four
different values to worry about instead of the eight different points of 𝖝.
Several different statistics can induce the same partition. In the above example, 𝑡(𝑥1 , 𝑥2 , 𝑥3 ) =
6(𝑥1 + 𝑥2 + 𝑥3 )2 and 𝑡(𝑥1 , 𝑥2 , 𝑥3 ) = 𝑥1 2 + 𝑥2 2 + 𝑥3 2 induces the same partition as
𝑡(𝑥1 , 𝑥2 , 𝑥3 ) = 𝑥1 + 𝑥2 + 𝑥3 .
A sufficient statistic is a particular kind of statistic. It is a statistic that condenses 𝖝 in such a way
that no “information about 𝜃” is lost. The only information about the parameter 𝜃 in the density
𝑓 (∙ ; 𝜃) from which we sampled, is contained in the sample 𝑋1 , … , 𝑋𝑛 ; so, when we say that a
statistic loses no information, we mean that it contains all the information about 𝜃 that is contained
in the sample.

Definition 1: Sufficient statistic


Let 𝑋1 , 𝑋2 , … , 𝑋𝑛 be a random sample from the density 𝑓(∙ ; 𝜃), where 𝜃 may be a vector. A statistic
𝑆 = 𝑠(𝑋1 , 𝑋2 , … , 𝑋𝑛 ) is defined to be a sufficient statistic if and only if the conditional distribution
of 𝑋1 , 𝑋2 , … , 𝑋𝑛 given 𝑆 = 𝑠 does not depend on 𝜃 for any value 𝑆 of 𝑠.
The definition says that a statistic 𝑆 = 𝑠(𝑋1 , 𝑋2 , … , 𝑋𝑛 ) is sufficient if the conditional distribution
of the sample given the value of the statistic does not depend on 𝜃.
We say 𝑇 is a sufficient statistic if the statistician who knows the value of 𝑇 can do just as good a
job of estimating the unknown parameter 𝜃 as the statistician who knows the entire random sample.
Example:
Let 𝑋1 , 𝑋2 , 𝑋3 be a sample size 3 from the Bernoulli distribution. Consider the two statistics 𝑆 =
𝑠(𝑋1 , 𝑋2 , 𝑋3 ) = 𝑋1 + 𝑋2 + 𝑋3 and 𝑇 = 𝑡(𝑋1 , 𝑋2 , 𝑋3 ) = 𝑋1 𝑋2 + 𝑋3 . We will show that 𝑆 is
sufficient and 𝑇 is not sufficient.
First, we should identify 𝖝
𝖝 = {(0, 0, 0), (1, 0, 0), (0, 1, 0), (0, 0, 1), (1, 0, 1), (1, 1, 0), (0, 1, 1), (1, 1, 1)}
Then, the conditional distributions for the values of 𝖝, 𝑆 and 𝑇 are calculated.
The conditional densities given in the last two columns are routinely calculated. For instance,
𝑓𝑥1,𝑥2,𝑥3|𝑆=1 (0, 1, 0|1) = 𝑃[𝑋1 = 0; 𝑋2 = 1; 𝑋3 = 0|𝑆 = 1]
𝑃[𝑋1 = 0; 𝑋2 = 1; 𝑋3 = 0|𝑆 = 1] ( 1 − 𝑝 ) 𝑝( 1 − 𝑝 ) 1
= = =
𝑃 [ 𝑆 = 1] 3 3
( ) 𝑝( 1 − 𝑝 ) 2
1
𝑓𝑥1,𝑥2,𝑥3|𝑇=0 (0, 1, 0|0) = 𝑃[𝑋1 = 0; 𝑋2 = 1; 𝑋3 = 0|𝑇 = 0]

𝑃[𝑋1 = 0; 𝑋2 = 1; 𝑋3 = 0|𝑆 = 1] ( 1 − 𝑝 ) 𝑝( 1 − 𝑝 ) 𝑝
= = 3 2
=
𝑃 [ 𝑇 = 0] (1 − 𝑝) + 2𝑝(1 − 𝑝) 1+𝑝
Similarly, the conditional distributions for the other values of 𝖝, 𝑆 and 𝑇 can be calculated and
shown in the below table.
𝖝 Values of 𝑆 Values of 𝑇 𝑓𝑥1 ,𝑥2 ,𝑥3 |𝑆 𝑓𝑥1 ,𝑥2 ,𝑥3 |𝑇
(0, 0, 0) 0 0 1 1−𝑝
1+𝑝
(1, 0, 0) 1 0 1 𝑝
3 1+𝑝
(0, 1, 0) 1 0 1 𝑝
3 1+𝑝
(0, 0, 1) 1 1 1 1−𝑝
3 1 + 2𝑝
(1, 0, 1) 2 1 1 𝑝
3 1 + 2𝑝
(1, 1, 0) 2 1 1 𝑝
3 1 + 2𝑝
(0, 1, 1) 2 1 1 𝑝
3 1 + 2𝑝
(1, 1, 1) 3 2 1 1

The conditional distribution of the sample given the values of 𝑆 independent of 𝑝. Thus 𝑆 is a
sufficient statistic. But 𝑇 is not a sufficient statistic since the conditional distribution of the
sample given the values of 𝑇 depends on 𝑝.

It is easy to see that if 𝑓 (𝑡) is a one to one function and 𝑇 is a sufficient statistic, then 𝑓(𝑇) is a
sufficient statistic. In particular we can multiply a sufficient statistic by a non-zero constant and
get another sufficient statistic.
Definition I of a sufficient statistic is not very workable. First, it does not tell us which statistic is
likely to be sufficient, and, second, it requires us to derive a conditional which may not be easy,
especially for continuous random variables.

Definition 2: Sufficient statistic


Let 𝑋1 , 𝑋2 , … , 𝑋𝑛 be a random sample from the density 𝑓(∙ ; 𝜃). A statistic 𝑆 = 𝑠(𝑋1 , 𝑋2 , … , 𝑋𝑛 ) is
defined to be a sufficient statistic if and only if the conditional distribution of 𝑇 given 𝑆 does not
depend on 𝜃 for any statistic 𝑇 = 𝑡(𝑋1 , 𝑋2 , … , 𝑋𝑛 ).
For some problems, there will exists jointly sufficient statistics. Specially when there are more
than one unknown parameters. In this case we need to use more than one statistic to get sufficiency.
We consider 𝑘 statistics 𝑆𝑖 = 𝑠𝑖 (𝑋1 , 𝑋2 , … , 𝑋𝑛 ), 𝑖 = 1, 2, ⋯ , 𝑘. We say 𝑆1 , … , 𝑆𝑛 are jointly
sufficient statistics if the statistician who knows the values of 𝑆1 , … , 𝑆𝑛 can do just as good a job
of estimating the unknown parameter 𝜃 as the statistician who knows the entire random sample.
In this setting 𝜃 typically represents several parameters and the number of statistics, 𝑘, is equal to
the number of unknown parameters.

Definition 3: Jointly sufficient statistics


Let 𝑋1 , 𝑋2 , … , 𝑋𝑛 be a random sample from the density 𝑓 (∙ ; 𝜃). The statistics 𝑆1 , … , 𝑆𝑛 are defined
to be jointly sufficient if the conditional distribution of 𝑋1 , … , 𝑋𝑛 given 𝑆1 = 𝑠1 , … , 𝑆𝑟 = 𝑠𝑟 does
not depend on 𝜃.
The sample 𝑋1 , … , 𝑋𝑛 itself is always jointly sufficient since the conditional distribution of the
sample given the sample does not depend on 𝜃.
Theorem
If 𝑆1 = 𝑠1 (𝑋1 , … , 𝑋𝑛 ) , … , 𝑆𝑟 = 𝑠𝑟 (𝑋1 , … , 𝑋𝑛 ) is a set of jointly sufficient statistics, then any set of
one-to-one functions, or transformations, of 𝑆1 , … , 𝑆𝑛 is also jointly sufficient.

For example, if ∑ 𝑋𝑖 and ∑ 𝑋𝑖 2 are jointly sufficient, then 𝑋̅ and ∑(𝑋𝑖 − 𝑋̅)2 = ∑ 𝑋𝑖 2 − 𝑛𝑋 2 are
also jointly sufficient.

Factorization Criterion
It is difficult to use the definition to check if a statistic is sufficient or to find a sufficient statistic
or set of statistics are sufficient.
Factorization theorem (Single sufficient statistic)
Let 𝑋1 , … , 𝑋𝑛 be a random sample of size 𝑛 from the density 𝑓(∙ ; 𝜃), where the parameter 𝜃 may
be a vector. A statistic 𝑆 = 𝑠(𝑋1 , … , 𝑋𝑛 ) is sufficient if and only if the joint density of 𝑋1 , … , 𝑋𝑛
which is ∏𝑛𝑖=1 𝑓 (𝑥𝑖 ; 𝜃), factors as
𝑓𝑥1,⋯,𝑥𝑛 (𝑥1 , ⋯ , 𝑥𝑛 ; 𝜃) = 𝑔(𝑠(𝑥1 , ⋯ , 𝑥𝑛 ); 𝜃)ℎ(𝑥1 , ⋯ , 𝑥𝑛 )

= 𝑔(𝑠; 𝜃)ℎ(𝑥1 , ⋯ , 𝑥𝑛 ),
Where the function ℎ(𝑥1 , ⋯ , 𝑥𝑛 ) is nonnegative and does not involve the parameter 𝜃 and the
function 𝑔(𝑠(𝑥1 , ⋯ , 𝑥𝑛 ); 𝜃) is nonnegative and depends on 𝑥1 , ⋯ , 𝑥𝑛 only through the function
𝑠(∙ , ⋯ , ∙).
Factorization theorem (Jointly sufficient statistic)
Let 𝑋1 , … , 𝑋𝑛 be a random sample of size 𝑛 from the density 𝑓(∙ ; 𝜃), where the parameter 𝜃 may
be a vector. A set of statistics 𝑆1 = 𝑠1 (𝑋1 , … , 𝑋𝑛 ) , … , 𝑆𝑟 = 𝑠𝑟 (𝑋1 , … , 𝑋𝑛 ) is jointly sufficient if
and only if the joint density of 𝑋1 , … , 𝑋𝑛 can be factored as
𝑓𝑥1 ,⋯,𝑥𝑛 (𝑥1 , ⋯ , 𝑥𝑛 ; 𝜃) = 𝑔(𝑠1 (𝑥1 , ⋯ , 𝑥𝑛 ), ⋯ , 𝑠𝑟 (𝑥1 , ⋯ , 𝑥𝑛 ); 𝜃)ℎ(𝑥1 , ⋯ , 𝑥𝑛 )
= 𝑔(𝑠1 , ⋯ , 𝑠𝑟 ; 𝜃)ℎ(𝑥1 , ⋯ , 𝑥𝑛 ) ,
Where the function ℎ(𝑥1 , ⋯ , 𝑥𝑛 ) is nonnegative and does not involve the parameter 𝜃 and the
function 𝑔(𝑠1 , ⋯ , 𝑠𝑟 ; 𝜃) is nonnegative and depends on 𝑥1 , ⋯ , 𝑥𝑛 only through the functions
𝑠1 (∙ , ⋯ , ∙), ⋯ , 𝑠𝑟 (∙ , ⋯ , ∙).
The above two theorems give us a relatively easy method for judging whether a certain statistic is
sufficient or a set of statistics is jointly sufficient.

Example 1:
Let 𝑋1 , 𝑋2 , … , 𝑋𝑛 be a random sample from the Bernoulli density with parameter 𝜃. Find a
sufficient statistic.
𝑓 (𝑥 ; 𝜃) = 𝑃[𝑋 = 𝑥 ] = 𝜃 𝑥 (1 − 𝜃)1−𝑥 ; 𝑥 = 0, 1
∏𝑛𝑖=1 𝑓 (𝑥𝑖 ; 𝜃) = ∏𝑛𝑖=1 𝜃 𝑥𝑖 (1 − 𝜃)1−𝑥𝑖
𝑛 𝑛
= 𝜃 ∑𝑖=1 𝑥𝑖 (1 − 𝜃)𝑛−∑𝑖=1 𝑥𝑖 × 1
= 𝑔(𝑠(𝑥1 , ⋯ , 𝑥𝑛 ); 𝜃)ℎ(𝑥1 , ⋯ , 𝑥𝑛 )
By factorization criterion,
sufficient statistic 𝑠(𝑥1 , ⋯ , 𝑥𝑛 ) = ∑𝑛𝑖=1 𝑥𝑖 and ℎ(𝑥1 , ⋯ , 𝑥𝑛 ) = 1.

Example 2:
Let 𝑋1 , 𝑋2 , … , 𝑋𝑛 be a random sample from the normal density with mean 𝜇 and variance unity.
Find a sufficient statistic.

Example 3:
Let 𝑋1 , 𝑋2 , … , 𝑋𝑛 be a random sample from the normal density 𝜙𝜇,𝜎2 (∙). Here the parameter 𝜃 is a
vector of two components; that is 𝜃 (𝜇, 𝜎). Find the jointly sufficient statistics.
2
1 1 𝑥−𝜇
𝑓 (𝑥 ; 𝜇, 𝜎 2)
= 𝑒 − ⁄2( 𝜎 )
√2𝜋 𝜎
𝑥 −𝜇 2
1 −1⁄2( 𝑖 )
∏𝑛𝑖=1 𝑓 (𝑥𝑖 ; 𝜇, 𝜎 2)
= 𝑛
∏𝑖=1 𝑒 𝜎
√2𝜋 𝜎
𝑛⁄
1 2 −1⁄ 2 ∑𝑛 (𝑥 −𝜇)2
= (2𝜋𝜎2 ) 𝑒 2𝜎 𝑖=1 𝑖

𝑛⁄
1 2 −1⁄ 2 [∑𝑛 2 𝑛 2
𝑖=1 𝑥𝑖 −2𝜇 ∑𝑖=1 𝑥𝑖 +𝑛𝜇 ]
= (2𝜋𝜎2 ) 𝑒 2𝜎
𝑛
∑ 𝑥 2 ∑𝑛 𝑥 𝜇 𝑛𝜇 2
𝑛⁄ [ 𝑖=12 𝑖 + 𝑖=12 𝑖 − 2 ]
1 2
= (2𝜋 ) 𝜎 −𝑛 𝑒 2𝜎 𝜎 2𝜎

= ℎ(𝑥1 , ⋯ , 𝑥𝑛 )𝑔(𝑠1 (𝑥1 , ⋯ , 𝑥𝑛 ), 𝑠2 (𝑥1 , ⋯ , 𝑥𝑛 ); 𝜇, 𝜎 2 )


By factorization criterion,
∑𝑛𝑖=1 𝑥𝑖 2 and ∑𝑛𝑖=1 𝑥𝑖 are jointly sufficient statistics for 𝜇 and 𝜎 2 .

Common questions

Powered by AI

Sufficient statistics induce partitions of the sample space, where each partition corresponds to a particular value that the statistic can take. This partitioning is essential because it reduces the sample space into equivalence classes where each class contains the same information about the parameter θ. The concept of sufficiency ensures that estimations based on these partitions are as effective as using the entire sample space, as no information about θ is lost during the condensation process .

Examples of transformations that maintain the sufficiency of a statistic include one-to-one transformations like linear scaling (e.g., multiplying by a constant) or adding a constant. Non-one-to-one transformations, such as squaring a statistic that can take both positive and negative values, might not maintain sufficiency because they could result in losing distinguishing information needed for estimation. For instance, both X and -X square to the same value, potentially merging distinct parameter information .

For a normal distribution with parameters μ and σ², the jointly sufficient statistics are ∑Xi and ∑Xi². By the factorization criterion, the joint density factors into a product involving these statistics and another function that does not involve θ. Hence, these statistics efficiently summarize the sample for the purposes of estimating both μ and σ² jointly, as any transformation of these statistics also retains sufficiency .

The factorization criterion is a method used to determine if a statistic is sufficient. It states that a statistic S is sufficient for parameter θ if the joint probability density function of the sample can be factored into a product of two functions: one that depends on the sample only through the statistic and θ, and another that does not depend on θ. Specifically, the joint density can be expressed as f(x1,...,xn;θ) = g(s(x1,...,xn);θ)h(x1,...,xn), where h(x1,...,xn) is nonnegative and independent of θ .

Transformations of statistics can maintain sufficiency through one-to-one transformations. If a statistic T is sufficient for θ, any one-to-one function of T will also be a sufficient statistic. This implies that multiplying a sufficient statistic by a non-zero constant or applying any reversible transformation preserves sufficiency. This principle indicates that sufficiency is an intrinsic property of the statistic's relationship with the parameter, not the specific form of the statistic itself .

The statistic S = X1 + X2 + X3 is sufficient for a Bernoulli distribution because the conditional distribution of the sample given S is independent of the parameter θ, thus fulfilling the condition for sufficiency. However, T = X1X2 + X3 is not sufficient because the conditional distribution of the sample given T depends on θ, violating the condition of sufficiency which requires independence from the parameter .

A sufficient statistic is a function of the sample that contains all the information necessary to estimate a parameter θ. It condenses the sample data without losing any information about θ, meaning that the statistic can provide the same estimation of θ as the full dataset. The formal definition states that a statistic S is sufficient if the conditional distribution of the sample given the statistic does not depend on θ .

For a sample X1,...,Xn from a Bernoulli distribution with parameter θ, the joint probability density function is ∏(θ^Xi * (1-θ)^(1-Xi)) for i = 1 to n. This can be factored into θ^∑Xi * (1-θ)^(n-∑Xi), where ∑Xi is the sum of the sample values. By the factorization criterion, the statistic ∑Xi is sufficient for θ, as the factor g(∑Xi;θ) depends on θ and the statistic, and h(x1,...,xn) = 1 does not involve θ .

The significance of h(x1,...,xn) in the factorization theorem is that it ensures the part of the joint density function that does not depend on the parameter θ. It's crucial because it confirms the statistic involves the parameter of interest only through the function g(s(x1,...,xn);θ). This separation allows the statistic to encapsulate all information about the parameter θ while h(x1,...,xn) maintains the structural integrity of the density without involving θ .

A jointly sufficient statistic involves more than one statistic to retain all information on a parameter vector θ that has multiple components. When multiple parameters need estimation, jointly sufficient statistics ensure that all information about these parameters is retained, much like a single sufficient statistic retains all information about a single parameter. Formally, statistics S1,...,Sk are jointly sufficient if the conditional distribution of the sample given all S1,...,Sk does not depend on θ .

You might also like