14.03.2014 Views

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

Modeling and Multivariate Methods - SAS

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

684 Statistical Details Appendix A<br />

Power Calculations<br />

Solving for Lb, a scalar, given an F for some α-level, like 0.05, <strong>and</strong> using a t as the square root of a<br />

one-degree-of-freedom F, making it properly two-sided, gives<br />

( Lb) LSV = t α ⁄ 2<br />

s LX'X ( ) – 1 L'<br />

For L testing some β i , this is<br />

b i<br />

LSV<br />

= t α ⁄ 2<br />

s (( X'X) – 1 ) ii<br />

which can be written with respect to the st<strong>and</strong>ard error as<br />

b i<br />

LSV<br />

=<br />

t α ⁄ 2<br />

stderr( b i<br />

)<br />

If you have a simple model of two means where the parameter of interest measures the difference between<br />

the means, this formula is the same as the LSD, least significant difference, from the literature<br />

LSD = t α ⁄ 2<br />

s<br />

-----<br />

1 1<br />

+ -----<br />

n 1<br />

n 2<br />

In the JMP Fit Model platform, the parameter for a nominal effect measures the difference to the average of<br />

the levels, not to the other mean. So, the LSV for the parameter is half the LSD for the differences of the<br />

means.<br />

Computations for the LSN<br />

The LSN solves for n in the equation:<br />

nδ<br />

α = 1 – ProbF --------------- 2<br />

, dfH,<br />

n – dfR – 1<br />

dfHσ 2<br />

where<br />

ProbF is the central F-distribution function<br />

dfH is the degrees of freedom for the hypothesis<br />

dfR is the degrees of freedom for regression in the whole model<br />

σ 2 is the error variance<br />

δ 2 is the squared effect size, which can be estimated by<br />

SS(H) -------------<br />

n<br />

When planning a study, the LSN serves as the lower bound.

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

Saved successfully!

Ooh no, something went wrong!