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

View source: R/SsPowerGivenJK.R

SsPowerGivenJKOrVarComR Documentation

Power given J, K and Obuchowski-Rockette variance components

Description

Power given J, K and Obuchowski-Rockette variance components

Usage

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 treatment-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

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

Details

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

Value

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

Examples

dataset <- dataset02 ## the pilot study
KStar <- length(dataset$ratings$NL[1,1,,1])
VarCom <- StSignificanceTesting(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)


RJafroc documentation built on Nov. 10, 2022, 5:45 p.m.