RJafroc (version 1.3.2)

SsPowerGivenJKOrVarComp: Power given J, K and Obuchowski-Rockette variance components

Description

Power given J, K and Obuchowski-Rockette variance components

Usage

SsPowerGivenJKOrVarComp(
  J,
  K,
  KStar,
  effectSize,
  varTR,
  cov1,
  cov2,
  cov3,
  varEps,
  alpha = 0.05,
  option
)

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 treatment-reader OR variance component

cov1

The OR cov1 covariance

cov2

The OR cov2 covariance

cov3

The OR cov3 covariance

varEps

The OR pure variance term

alpha

The size of the test (default = 0.05)

option

The desired generalization ("RRRC", "FRRC", "RRFC", "ALL")

Value

A list object containing the estimated power and associated statistics for each desired generalization.

Details

The variance components are obtained using StSignificanceTesting with method = "ORH".

Examples

# NOT RUN {
dataset <- dataset02 ## the pilot study
KStar <- length(dataset$NL[1,1,,1])
VarComp <- StSignificanceTesting(dataset, FOM = "Wilcoxon", method = "ORH", option = "RRRC")$varComp
varTR <- VarComp$varTR
cov1 <- VarComp$cov1
cov2 <- VarComp$cov2
cov3 <- VarComp$cov3
varEps <- VarComp$var
ret <- SsPowerGivenJKOrVarComp (J = 5, K = 100, KStar = KStar,  
   effectSize = 0.05, varTR, cov1, cov2, cov3, varEps, option = "RRRC")
    
cat("RRRC power = ", ret$powerRRRC)

  
# }