Sensitivity Analysis for Sample Size Planning From the AIPE Perspective for a Pearson Correlation
Source:R/ss_aipe_r_sensitivity.R
ss_aipe_r_sensitivity.RdQuantifies how much misspecification of the population Pearson
correlation distorts an AIPE-based sample size plan. On each
replication the function draws an n-row sample from a
bivariate normal distribution with correlation true_rho and
computes the sample correlation and its Fisher's \(Z\) CI, the
interval ss_aipe_r plans for and
correlations_test reports. Because the back-transformed
width is largest at \(\rho = 0\) and shrinks as \(|\rho|\) grows,
a planning value whose magnitude overstates the population
correlation yields realized intervals wider than planned, and the
summary rows report by how much.
Usage
ss_aipe_r_sensitivity(
true_rho = NULL,
estimated_rho = NULL,
width,
specified_N = NULL,
conf_level = 0.95,
assurance = NULL,
G = 1000,
print_iter = FALSE,
save = FALSE,
filename = "ss_aipe_r_sensitivity_result.csv"
)Arguments
- true_rho
Population Pearson correlation; must lie in \((-1, 1)\).
- estimated_rho
Planning value of the correlation passed to
ss_aipe_r; supply this orspecified_Nbut not both.- width
Desired full width of the CI on the correlation.
- specified_N
Sample size to evaluate (incompatible with
estimated_rho).- conf_level
Confidence level (default
0.95).- assurance
Optional assurance probability passed to
ss_aipe_r.- G
Number of Monte Carlo replications (default 1000).
- print_iter
Logical. Print iteration index per replication.
- save
Logical. If
TRUEwrite per-replication results tofilename.- filename
Path used when
save = TRUE.
Value
A data.frame with rows for the realized correlation,
the interval width, the proportion of intervals at or below
width, tail-specific and overall non-coverage of
true_rho, and the input echoes, including assurance (present only when an
assurance was supplied).
References
Bonett, D. G., & Wright, T. A. (2000). Sample size requirements for estimating Pearson, Kendall and Spearman correlations. Psychometrika, 65(1), 23–28. doi:10.1007/BF02294183
Kelley, K., & Maxwell, S. E. (2003). Sample size for multiple regression: Obtaining regression coefficients that are accurate, not simply significant. Psychological Methods, 8(3), 305–321. doi:10.1037/1082-989X.8.3.305
Maxwell, S. E., Delaney, H. D., & Kelley, K. (2027). Designing experiments and analyzing data: A model comparison perspective (4th ed.). Routledge.
See also
ss_aipe_r, ss_aipe_partial_r_sensitivity
design_consequences for what a chosen design delivers:
power, the Type S (sign) and Type M (exaggeration) errors of the
significance filter, and the expected confidence interval width.
Other AIPE sample size planning:
ss_aipe_c_sensitivity(),
ss_aipe_cliff_delta(),
ss_aipe_cliff_delta_sensitivity(),
ss_aipe_composite_sem(),
ss_aipe_equivalence_r(),
ss_aipe_equivalence_r_sensitivity(),
ss_aipe_equivalence_smd(),
ss_aipe_equivalence_smd_sensitivity(),
ss_aipe_icc(),
ss_aipe_icc_sensitivity(),
ss_aipe_indirect_effect(),
ss_aipe_indirect_effect_sensitivity(),
ss_aipe_mixed_effects_sensitivity(),
ss_aipe_omega_squared(),
ss_aipe_omega_squared_sensitivity(),
ss_aipe_partial_r(),
ss_aipe_partial_r_sensitivity(),
ss_aipe_pcm_sensitivity(),
ss_aipe_r(),
ss_aipe_reliability_sensitivity(),
ss_aipe_semipartial_r(),
ss_aipe_semipartial_r_sensitivity()
Author
Ken Kelley kkelley@nd.edu
Examples
# Reduced replications and a wide target interval keep this fast.
set.seed(113)
ss_aipe_r_sensitivity(
true_rho = 0.30, estimated_rho = 0.30, width = 0.40,
G = 50, print_iter = FALSE
)
#> term value
#> mean_r 0.289
#> median_r 0.278
#> sd_r 0.0867
#> mean_ci_width 0.399
#> median_ci_width 0.404
#> sd_ci_width 0.0219
#> pct_ci_less_w 0.44
#> pct_ci_miss_low 0
#> pct_ci_miss_high 0
#> total_type_I_error 0
#> total_N 81
#> true_rho 0.3
#> estimated_rho 0.3
#> width 0.4
#> conf_level 0.95
#>
#> Confidence level: 95%