pSceptical | R Documentation |
Computes sceptical p-values and z-values based on the z-values of the original and the replication study and the corresponding variance ratios. If specified, the sceptical p-values are recalibrated.
pSceptical(
zo,
zr,
c,
alternative = c("one.sided", "two.sided"),
type = c("golden", "nominal", "controlled")
)
zSceptical(zo, zr, c)
zo |
Numeric vector of z-values from original studies. |
zr |
Numeric vector of z-values from replication 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. |
alternative |
Either "one.sided" (default) or "two.sided". If "one.sided", the sceptical p-value is based on a one-sided assessment of replication success in the direction of the original effect estimate. If "two.sided", the sceptical p-value is based on a two-sided assessment of replication success regardless of the direction of the original and replication effect estimate. |
type |
Type of recalibration. Can be either "golden" (default), "nominal",
or "controlled". Setting |
pSceptical
is the vectorized version of
the internal function .pSceptical_
.
Vectorize
is used to vectorize the function.
pSceptical
returns the sceptical p-value.
zSceptical
returns the z-value of the sceptical p-value.
Leonhard Held
Held, L. (2020). A new standard for the analysis and design of replication studies (with discussion). Journal of the Royal Statistical Society: Series A (Statistics in Society), 183, 431-448. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1111/rssa.12493")}
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. \Sexpr[results=rd]{tools:::Rd_expr_doi("10.1214/21-AOAS1502")}
Micheloud, C., Balabdaoui, F., Held, L. (2023). Beyond the two-trials rule: Type-I error control and sample size planning with the sceptical p-value. https://arxiv.org/abs/2207.00464
sampleSizeReplicationSuccess
, powerReplicationSuccess
, levelSceptical
## no recalibration (type = "nominal") as in Held (2020)
pSceptical(zo = p2z(0.01), zr = p2z(0.02), c = 2, alternative = "one.sided",
type = "nominal")
## recalibration with golden level as in Held, Micheloud, Pawel (2020)
pSceptical(zo = p2z(0.01), zr = p2z(0.02), c = 2, alternative = "one.sided",
type = "golden")
## two-sided p-values 0.01 and 0.02, relative sample size 2
pSceptical(zo = p2z(0.01), zr = p2z(0.02), c = 2, alternative = "one.sided")
## reverse the studies
pSceptical(zo = p2z(0.02), zr = p2z(0.01), c = 1/2, alternative = "one.sided")
## both p-values 0.01, relative sample size 2
pSceptical(zo = p2z(0.01), zr = p2z(0.01), c = 2, alternative = "two.sided")
zSceptical(zo = 2, zr = 3, c = 2)
zSceptical(zo = 3, zr = 2, c = 2)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.