Two sample problems. 1) Confidence intervals and tests for difference of two means:

November 19, 2017 | Author: Roger Rose | Category: N/A
Share Embed Donate


Short Description

1 Two sample problems Chapters: 17 and 19. Outline of presentation: 1) Confidence intervals and tests for difference of ...

Description

Two sample problems Chapters: 17 and 19. Outline of presentation: 1) Confidence intervals and tests for difference of two means: µ1 − µ2 2) Confidence intervals and tests for difference of two proportions: p1 − p2 Summary: all intervals of the form estimate ± multiplier × SE where SE is standard error of the estimate. All tests of the form estimate − target standard error 202

An example: 8 animals treated with a steroid, 10 controls. Measured body weight gain. Two samples: n1 = 8, n2 = 10. Question: does steroid change body weight gain? Data: two sample means: x ¯1 = 32.8, x ¯2 = 40.5. Two sample sds: s1 = 2.6, s2 = 2.5. Step 1: recognize that question has yes/no answer. So: test hypothesis.

203

Introduce notation. 1) µ1 is population mean weight gain of steroid treated animals. 2) µ2 is population mean weight gain of untreated animals. Frame null hypothesis: µ1 = µ2. Equivalent hypothesis: µ1 − µ2 = 0. Examine question to choose alternative: Alternative does not predict direction so: Alternative is µ1 6= µ2.

204

Develop test statistic: Numerator is measure of discrepancy: (¯ x1 − x ¯2) − target value of µ1 − µ2 This is just x ¯1 − x ¯2 Need standard error for difference. Basic formula: SE of difference of two independent quantities is: q

2 SE2 + SE 1 2

So: standard error of x ¯1 − x ¯2 is v u 2 2 uσ t 1 + σ2 n1 n2

Normally we must estimate this using v u 2 us t 1 + n1

s2 2 n2 205

This leads to the test statistic: x ¯1 − x ¯2 r t= 2 s2 1 + s2 n n 1

2

In our example get: 32.8 − 40.5

t=r

2.62 8

2.52

+ 10

= −6.35

How do we compute a P -value? Solution used in text: from t tables taking df = min{n1 − 1, n2 − 1} = 7 in our case. In table C for 7 df largest value of t is 5.408 corresponding to P = 0.001 so P < 0.001 From software P ≈ 0.00038. Conclusion: strong (very highly significant) evidence that steroid affects average weight gain. 206

Confidence interval for mean difference: µ1 − µ2 ? Same formula as always: v u 2 us ∗ x ¯1 − x ¯2 ± t t 1 + n1

s2 2 n2

In our example with 7 df get t∗ = 2.365 for 95% CI. Also get standard error: v u 2 us t 1 + n1

s2 2 = 1.212 n2

Degrees of Freedom; Controversy There is not universal agreement on how to do this test. Book gives two options and dismisses a third.

207

Option 1: Easy: use t statistic as above and df as above. Option 2: Satterthwaite’s approximation. Use t as above but  2 2 s s2 2 1 n1 + n2 df =  2 2  2 2 s s2 1 1 1 + n1−1 n1 n2 −1 n2 

Option 3: Assume σ1 = σ2 and use pooled estimate of standard error: v u u (n1 − 1)s2 + (n2 − 1)s2 1 2 t

n1 + n2 − 2

1 1 + n1 n2

!

and then take

df = n1 + n2 − 2

208

In our example pooling produces the standard error s

7 ∗ 2.62 + 9 ∗ 2.52 1 1 + = 1.207 16 8 10 



and t = −6.38 with df = 16. The P value becomes much smaller, however. From software P = 9.1 × 10−6 Option 2 gives df = 14.86 and P = 1.36×10−5. Commentary: 1) Software usually does option 3 by default. 2) Better software also produces Option 2. 3) In this case not much difference in conclusions. 209

Comparing two proportions. Example: two samples of praying mantis. Brown: 65; Green: 45. Of brown: 45 put on green leaves. Of green: 25 put on brown leaves. After 3 weeks: of the 45 brown on green leaves 26 still alive. Of the 25 green on brown leaves 16 still alive. Question: difference in survival rates? Common presentation of results: Contingency table.

210

Insect Type Status Brown Green Total Alive 26 16 42 Dead 19 9 28 Total 45 25 70 Model: Let X1 be surviving number of brown. Let X2 be surviving number of green. Each of X1, X2 is Binomial. Numbers of trials n1 = 45, n2 = 25. Population survival probabilities: p1, p2. Null hypothesis p1 = p2. Alternative p1 6= p2. Test statistic: pˆ1 − pˆ2 z=r   1 1 pˆ(1 − pˆ) n + n 1 2

Note: in denominator pˆ = (X1 + X2)/(n1 + n2) is overall success rate. Called pooled estimate. 211

In our case: 26 16 pˆ1 = = 0.578 pˆ2 = = 0.64 45 25 and 26 + 16 = 0.6 pˆ = 45 + 25 This gives −0.0622 z=r   = −0.51 1 + 1 0.6 × 0.4 45 25

Get P -value from normal tables: two sided. P = 0.61 Interpretation: not much evidence of a difference in survival rates. Ref: di Cesnola, A.P. (1904) Biometrika, 4, 58–59.

212

Confidence interval for p1 − p2: pˆ1 − pˆ2 ± z ∗

s

pˆ (1 − pˆ2) pˆ1(1 − pˆ1) + 2 n1 n2

Notice: No pooling. (In testing, pooling justified by null hypothesis.) Commentary: text recommends: add 1 to each Xi and 2 to each ni then do all arithmetic as above. Not standard. Improves coverage probability. Large sample methods not recommended unless all of n1p1, n1(1 − p1), n2p2, n2(1 − p2) large enough. Book recommends all be at least 10. Judged by all cell counts at least 10 in contingency table. 213

Matched pairs designs: instead of 2 independent samples, have 1 sample of pairs. Example: look back at cross-fertilization of peas example. Originally would have had 2 measurements for each parent. Data reduced by subtraction to 1 sample problem! Example: Pearson Lee data on father / son height. Consists of N = 1078 pairs. Denote: Fi father’s height and Si son’s height in ith pair. Problem: are sons taller than fathers? Idea: µ1 is population average height of sons. µ2 pop average height of fathers. (At point in time when data collected!) 214

Point of next piece: illustrate merit of matched pairs design. Treat the N = 1078 pairs as the population. Then µ1 = 68.68, µ2 = 67.69 , σ1 = 2.81 and σ2 = 2.74. In the population the variables F and S are correlated: ρ = 0.502 Consider two methods of comparing µ1 and µ2 based on sampling. Method 1: take two samples of size n1 = n2 = 9, one of Fathers, other of Sons. Method 2: take one sample of n = 9 pairs of Fathers and sons. 215

An explicit example: For Method I: drew following 2 independent samples Family # i 128 251 756 150 257 ..

Fi 70.01 68.32 65.24 69.52 64.07 ..

Family # i 635 574 564 778 160 ..

Si 78.25 70.70 69.20 69.12 70.82 ..

Drew total of 1000 samples of n = 9 fathers and 1000 samples of n = 9 sons.

216

For Method II I drew the following sample of pairs: Family # i 851 53 919 475 754 ..

Fi 69.07 65.83 65.68 64.68 64.34 ..

Si 78.36 67.07 67.68 66.79 69.23 ..

Si − Fi 9.29 1.24 2.00 2.11 4.88 ..

Repeated this 1000 times.

217

¯ − F¯. Here is a histogram of S Solid lines: two independent samples.

0.0

0.1

0.2

0.3

0.4

Dotted lines: sample of pairs.

−2

0

2

4

Difference in Heights (inches)

218

6

Numerical summary of this Monte Carlo experiment. Method 1 outcomes: F¯ 68.15 66.63 68.11 68.48 67.17 ..

¯ S 69.78 9.98 7.98 69.36 67.32 ..

¯ − F¯ S 1.63 3.35 -0.13 0.88 0.15 ..

Method 2 outcomes: F¯ 67.76 68.66 68.96 66.33 68.33 ..

¯ S 67.68 68.01 69.15 67.72 67.90 ..

¯ − F¯ S -0.07 -0.65 0.20 1.39 -0.43 ..

219

To compare: examine mean and sd of the last columns: Get Independent Mean SD 1.046 1.302

Matched Pair Mean SD 0.958 0.932

Major point: both means close to µ1 − µ2 = 0.997. But: SD for matched pairs is smaller. Formula for SE of difference of two independent means: v u 2 2 uσ t 1 + σ2 = 1.308 n1 n2

Formula for SE of difference in paired sample: s

σ12 + σ22 − 2ρ ∗ σ1σ2 = 0.925 n

Notice great match of theory to Monte Carlo. 220

Example problem: Does too much sleep impair intellectual performance. 10 subjects tested twice each. Once after two normal night’s sleep , Once after two nights of ‘extended sleep’. Data on test for vigilance: low scores are alert: Subject Normal Extended Diff

1 8 8 0

2 9 9 0

3 14 15 -1

4 4 2 2

5 12 21 -9

6 11 16 -5

7 3 9 -6

8 26 38 -12

9 3 10 -7

10 11 11 0

WARNING: I might put in a column of differences even if data are not paired.

221

Null: pop mean difference µN − µE in vigilance scores is 0. Alternative: µN < µE . ¯ −E ¯ = −3.8; s = 4.66. Summary statistics: N Test statistic: −3.8 − 0 √ t= = −2.58 4.66/ 10 One sided alternative. 9df.

P -value in lower tail.

P = 0.015 In tables best approx is 0.01 < P < 0.02.

222

What if: had used 20 subjects. 10 assigned to Normal, 10 to Extended at random? Could have presented same data (but probably without row ‘Subject’). Analysis: not paired, so 2 sample t test. Hypotheses unchanged! x ¯N = 10.1, sN = 6.81, x ¯E = 13.9, sE = 9.92 Two sample t statistic is 10.1 − 13.9

t=r

6.812 10

9.922

+ 10

= −1.00

which gives P = 0.172 In tables 0.15 < P < 0.2. Not significant.

223

Summary points: 1) for original description of experiment paired analysis right, two sample analysis wrong. (Only 10 subjects.) 2) Since two variables positively correlated paired design is better. 3) conclusion is that extra sleep does seem to worsen vigilance. 4) but if we had collected same data in unpaired design would have concluded no real evidence that extra sleep worsens vigilance.

224

Another example: Studying gopher tortoise burrows to see which are active. Two methods of evaluation of ‘active’ compared. Camera versus ‘experience’. Data: 151 burrows judged by ‘experience’. 107 rated active. 114 judged by cameras. 48 rated active. Problem: evaluation methods equivalent? Assume: burrows assigned to evaluation method at random.

225

If X1 is number judged active by experience then X1 is Binomial with n1 = 151, and some p1. We estimate pˆ1 = X1/n1 = 107/151 = 0.7086. Then X2 number judged active by camera is Binomial n2 = 114, pˆ2 = 48/114 = 0.44211 Null hypothesis: p1 = p2. Alternative: p1 6= p2. Pooled estimate of p assuming p1 = p2 is 107 + 48 pˆ = = 0.5849. 151 + 114 Test statistic: pˆ1 − pˆ2 z=r   = 4.70 1 + 1 0.5849(1 − 0.5849) 151 114

Get two sided P -value; less than 0.006 in Table A. 226

View more...

Comments

Copyright � 2017 SILO Inc.