ase(x) pase(q,k) qase(p,k)
An approach to this problem which motivates the use of robust estimators of scale is as follows: think of the estimated effects as a sample from a zero mean normal distribution (the null effects) contaminated by the non-null effects. Use robust methods to find an estimate of the scale of the null effects that is insensitive to the non-null effects. Then the estimated effects that are large compared to this scale estimate correspond to the non-null effects.
Dong (1993) proposes the adaptive standard error (ASE) as a robust scale estimator for this problem. Haaland and O'Connell (1994) studied the properties of this and several related tests. The ASE based test is recommended when there is a priori reason to believe that there will be only a few significant effects, say 0 to 3. However, the PSE (pseudo standard error -- Lenth, 1989) is recommended as an all around good test for identifying significant effects in a saturated fractional factorial design.
The value of the ASE is included in the fac.aov object created in the standard analysis of a fractional factorial design in S+DOX. The reference distribution is used to provide approximate p-values in the summary procedure and to draw a cut-off line for significant effects on the pareto and half-normal plots. The estimated ASE is equal to 1/slope of the line through the null effects on the half-normal plot. Tests based on the ASE are also used in the empirical Bayes plot.
Dong, F. (1993). "On the Idendification of Active Contrasts in Unreplicated Fractional Factorials." Statistica Sinica 3, 209-217.
Haaland, P. D. and M. A. O'Connell (1994), "Inference for Effect Saturated Fractional Factorials", to apear in Technometrics.
Lenth, R. V. (1989), "Quick and Easy Analysis of Unreplicated Fractional Factorials", Technometrics, 31, 469-473.
Zahn, D. A. (1975). "An Empirical Study of the Half-Normal Plot." Technometrics 17, 201-211.
buffer.fac <- fac.aov(buffer.df) buffer.fac$ase$ ase(buffer.fac$feffects$) qase(.95,15) summary(buffer.fac,method="ase") pareto(buffer.fac,method="ase") qqnorm(buffer.fac,method="ase") ebplot(buffer.fac,method="ase")