effectSizeReplicationSuccess {ReplicationSuccess} | R Documentation |
The minimum relative effect size (replication to original) to achieve replication success with the sceptical p-value is computed based on the result of the original study and the corresponding variance ratio.
effectSizeReplicationSuccess(
zo,
c = 1,
level = 0.025,
alternative = c("one.sided", "two.sided"),
type = c("golden", "nominal", "controlled")
)
zo |
Numeric vector of z-values from original studies. |
c |
Numeric vector of variance ratios of the original and replication effect estimates. This is usually the ratio of the sample size of the replication study to the sample size of the original study. |
level |
Threshold for the calibrated sceptical p-value. Default is 0.025. |
alternative |
Specifies if |
type |
Type of recalibration. Can be either "golden" (default),
"nominal" (no recalibration), or "controlled". "golden" ensures that for an
original study just significant at the specified |
effectSizeReplicationSuccess
is the vectorized version of
the internal function .effectSizeReplicationSuccess_
.
Vectorize
is used to vectorize the function.
The minimum relative effect size to achieve replication success with the sceptical p-value.
Leonhard Held, Charlotte Micheloud, Samuel Pawel, Florian Gerber
Held, L., Micheloud, C., Pawel, S. (2022). The assessment of replication success based on relative effect size. The Annals of Applied Statistics. 16:706-720. doi:10.1214/21-AOAS1502
Micheloud, C., Balabdaoui, F., Held, L. (2023). Assessing replicability with the sceptical p-value: Type-I error control and sample size planning. Statistica Neerlandica. doi:10.1111/stan.12312
sampleSizeReplicationSuccess
, levelSceptical
po <- c(0.001, 0.002, 0.01, 0.02, 0.025)
zo <- p2z(po, alternative = "one.sided")
effectSizeReplicationSuccess(zo = zo, c = 1, level = 0.025,
alternative = "one.sided", type = "golden")
effectSizeReplicationSuccess(zo = zo, c = 10, level = 0.025,
alternative = "one.sided", type = "golden")
effectSizeReplicationSuccess(zo = zo, c = 10, level = 0.025,
alternative = "one.sided", type = "controlled")
effectSizeReplicationSuccess(zo = zo, c= 2, level = 0.025,
alternative = "one.sided", type = "nominal")
effectSizeReplicationSuccess(zo = zo, c = 2, level = 0.05,
alternative = "two.sided", type = "nominal")