14.03.2014 Views

Basic Analysis and Graphing - SAS

Basic Analysis and Graphing - SAS

Basic Analysis and Graphing - SAS

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Chapter 2 Performing Univariate <strong>Analysis</strong> 83<br />

Statistical Details for the Distribution Platform<br />

Binomial<br />

The Binomial option accepts data in two formats: a constant sample size, or a column containing sample<br />

sizes.<br />

pmf: n p for ; x = 0,1,2,...,n<br />

x<br />

x ( 1 – p) n – x 0 ≤ p ≤ 1<br />

E(x) = np<br />

Var(x) = np(1-p)<br />

where n is the number of independent trials.<br />

Note: The confidence interval for the binomial parameter is a Score interval. See Agresti (1998).<br />

Beta Binomial<br />

This distribution is useful when the data is a combination of several Binomial(p) distributions, each with a<br />

different p. One example is the overall number of defects combined from multiple manufacturing lines,<br />

when the mean number of defects (p) varies between the lines.<br />

The Beta Binomial distribution results from assuming that x|π follows a Binomial(n,π) distribution <strong>and</strong> π<br />

follows a Beta(α,β). The Beta Binomial has parameters p = α/(α+β) <strong>and</strong> δ =1/(α+β+1). The δ is a<br />

dispersion parameter. When δ > 0, there is over dispersion, meaning there is more variation in x than<br />

explained by the Binomial alone. When δ < 0, there is under dispersion. When δ =0, x is distributed as<br />

Binomial(n,p). The Beta Binomial only exists when n ≥ 2 .<br />

pmf:<br />

n<br />

Γ 1<br />

<br />

δ -- – 1 1 Γ x + p -- – 1<br />

Γ n – x + ( 1 – p) 1<br />

δ<br />

<br />

-- –<br />

δ 1 <br />

<br />

--------------------------------------------------------------------------------------------------------------------------<br />

x<br />

Γ p 1 -- –<br />

δ 1 1 Γ ( 1 – p)<br />

-- – 1 Γ 1 n + -- – 1<br />

δ<br />

δ <br />

p 1 – p<br />

for 0 ≤ p ≤1<br />

; max( –------------------- ,–--------------------<br />

) ≤ δ ≤ 1 ; x = 0,1,2,...,n<br />

n – p – 1 n – 2 + p<br />

E(x) = np<br />

Var(x) = np(1-p)[1+(n-1)δ]<br />

where<br />

Γ ( . )<br />

is the Gamma function.<br />

Remember that x|π ~ Binomial(n,π), while π ~ Beta(α,β). The parameters p = α/(α+β) <strong>and</strong> δ =1/(α+β+1)<br />

are estimated by the platform. To obtain estimates of α <strong>and</strong> β, use the following formulas:<br />

αˆ<br />

1 – δˆ<br />

= pˆ -----------<br />

<br />

<br />

δˆ<br />

βˆ<br />

( 1 – pˆ<br />

) 1 – δˆ<br />

= -----------<br />

<br />

<br />

δˆ

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!