View source: R/SsPowerGivenJK.R
SsPowerGivenJKOrVarCom | R Documentation |
Power given J, K and Obuchowski-Rockette variance components
SsPowerGivenJKOrVarCom(
J,
K,
KStar,
effectSize,
VarTR,
Cov1,
Cov2,
Cov3,
Var,
alpha = 0.05,
analysisOption = "RRRC"
)
J |
The number of readers in the pivotal study |
K |
The number of cases in the pivotal study |
KStar |
The number of cases in the pilot study |
effectSize |
The effect size |
VarTR |
The modality-reader OR variance component |
Cov1 |
The OR Cov1 covariance |
Cov2 |
The OR Cov2 covariance |
Cov3 |
The OR Cov3 covariance |
Var |
The OR pure variance term |
alpha |
The size of the test (default = 0.05) |
analysisOption |
Specify the random factor(s): "RRRC", "FRRC", "RRFC" |
The variance components are obtained using St
with method = "OR"
.
A list containing the estimated power and associated statistics for the specified random factor(s).
dataset <- dataset02 ## the pilot study
KStar <- length(dataset$ratings$NL[1,1,,1])
VarCom <- St(dataset, FOM = "Wilcoxon",
method = "OR", analysisOption = "RRRC")$ANOVA$VarCom
VarTR <- VarCom["VarTR",1]
Cov1 <- VarCom["Cov1",1]
Cov2 <- VarCom["Cov2",1]
Cov3 <- VarCom["Cov3",1]
Var <- VarCom["Var",1]
ret <- SsPowerGivenJKOrVarCom (J = 5, K = 100, KStar = KStar,
effectSize = 0.05, VarTR, Cov1, Cov2, Cov3, Var, analysisOption = "RRRC")
##cat("RRRC power = ", ret$powerRRRC)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.