FDRpenalties {APCanalysis} | R Documentation |
FDRpenalties() generates the set of penalties used for the APC criterion so that the false discovery rate (FDR) is controlled at the specified level. The values of the penalties are estimated using a Monte Carlo simulation procedure and also depend on the run size of the experiment, the number of candidate effects and the maximum model size.
FDRpenalties(n, k = n - 1, m = min(n - 2, k), fdr = .1, reps = 50000, rnd = 3)
n |
The number of experimental runs in the study, i.e. the row dimension of the design matrix (for orthogonal 2-level designs n must be a multiple of 4). |
k |
The number of candidate effects under study, i.e., the column dimension of the design matrix (1 < k <= n-1). |
m |
The maximum size of the candidate models (0 < m < min(n-2, k)). |
fdr |
The level (0 < fdr < 1) at which the experimentwise error rate will be controlled (default is fdr = 0.1). |
reps |
The number of repetitions used by the Monte Carlo simulation algorithm which estimates the set of penalties (default is 50000). |
rnd |
The number of decimal places returned for the estimated penalties (default is rnd = 3). |
A vector of containing the m + 1 penalties for the APC procedure that controls the FDR at the specified level.
Penalties are estimated using Monte Carlo simulations and thus the estimates will not be exactly the same each time the function is run. The precision of the estimates can be increased by increasing the number of reps but the function will take longer to run. The amount of time needed to run this programme increases as the values of n, k and m increase. For larger experiments it may be necessary to reduce the number of reps.
Arden Miller and Abu Zar Md. Shafiullah
Miller A.E. (2005) _The analysis of unreplicated factorial experiments using all possible comparisons_ Technometrics, 47, 51-63, 2005.
apc
, IERpenalties
, EERpenalties
.
## Penalties for a 8-run experiment that has 5 candidate effects are generated. ## The maximum model size is set to 5 and a false discovery rate of .05 is used. FDRpenalties(n = 8, k = 5, m = 5, fdr = .05, reps = 12000)