ref.dist {unrepx} | R Documentation |
Analysis of effects from screening experiments
Description
These functions facilitate making inferences based on effect estimates in an unreplicated experiment, withn an underlying effect-sparsity model.
Usage
ref.dist(method, n.effects, nsets, save = TRUE)
eff.test(effects, method = "Zahn", pareto = TRUE, refdist, save = TRUE)
Arguments
method |
The method to use in determining the reference line, curve, and/or critical values. This must be the name of a provided pseudo-standard-error method (see |
n.effects |
Integer number of effects estimated. |
nsets |
The number of complete-null samples of size |
save |
Logical value. If |
effects |
Vector of observed effects to be tested against the reference distribution. |
pareto |
Logical value. If |
refdist |
A result of a previous call to |
Details
ref.dist
simulates samples of effects from the standard normal dstribution. For each sample, the pseudo standaerd error (PSE) of the effects (typically some kind of outlier-resistant estimate of the SD) is obtained via a call to PSE
with specified method. The absolute values are obtained as ratios of the simulated effects and the PSE, as well as and the maxima of these absolute
values. Quantiles and tail areas of these simulated distributions then form a reference for obtaining critical values and P values in testing an observed sample of effects.
eff.test
performs a traditional-style analysis for an observed sample of effects. It outputs the effects, PSE, ratios; and uses tail areas of the associated reference distribution to compute individual and simultaneous
values. The simultaneous
values implement a multiplicity correction for any type-I errors occurring among the tests.
Value
ref.dist
returns an object of class "eff_refdist"
– structurally, a list
with elements abst
(the absolute values of the simulated statistics),
max.abst
(the sample maxima of abst
), and sig
(a signature of the form method_n.effects
). There is a print
method for this class that displays a summary.
eff.test
returns a data.frame
containing the estimates, , and estimated P values as tail areas of
abst
and max.abst
from the reference distribution.
Author(s)
Russell V. Lenth
Examples
require("unrepx")
zahn15 <- ref.dist("Zahn", 15)
eff.test(pdEff, refdist = zahn15)