Power given J, K and Obuchowski-Rockette variance components

Power given J, K and Obuchowski-Rockette variance components

SsPowerGivenJKOrVarCom(
  J,
  K,
  KStar,
  effectSize,
  VarTR,
  Cov1,
  Cov2,
  Cov3,
  Var,
  alpha = 0.05,
  analysisOption = "RRRC"
)

Arguments

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"

Value

A list containing the estimated power and associated statistics for the specified random factor(s).

Details

The variance components are obtained using St with method = "OR".

Examples

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)