simMRMC {iMRMC}

Simulate an MRMC data set

Description

This program simulates observations from one set of readers scoring one set of cases. It produces one modality and one truth state of ROC data following Roe1997_Acad-Radiol_v4p298 and Roe1997_Acad-Radiol_v4p587. In order to produce an entire ROC data set, please use sim.gRoeMetz.

Usage

simMRMC(simMRMC.config)

Arguments

simMRMC.config

[list] of simulation parameters:

• modalityID [character] label modalityID

• readerIDs [factor] the ID of each reader

• caseIDs [factor] the ID of each case

• mu [numeric] mean

• var_r [numeric] variance of random reader effect

• var_c [numeric] variance of random case effect

• var_rc [numeric] variance of random reader by case effect

Details

The simulation is a linear model with one fixed effect and three normally distributed independent random effects corresponding to readers, cases, and an interaction between the two.

L.rc = mu + readerEffect.r + caseEffect.c + readerXcaseEffect.rc

Value

L [data.frame] with nC*nR rows of 4 variables

• L$modalityID [factor] determined by input modalityID

• L$readerID [factor] determined by input readerIDs

• L$caseID [factor] determined by input caseIDs

• L$score [numeric] R.r + C.c + RC.rc

  • r = 1,2,...,nR

  • c = 1,2,...,nC

  • R.r ~ N(0,var_r)

  • C.c ~ N(0,var_c)

  • RC.rc ~ N(0,var_rc)