A list consists of the following
integer vectors f, h
and integers NL, NI, C
.
-
f
Non-negative integer vector specifying number of false alarms associated with each confidence level. The first component corresponding to the highest confidence level.
-
h
Non-negative integer vector specifying number of Hits associated with each confidence level. The first component corresponding to the highest confidence level.
-
NL
A positive integer, representing Number of Lesions.
-
NI
A positive integer, representing Number of Images.
-
C
A positive integer, representing Number of Confidence level.
Contents:
A single reader and single modality case
------------------------------------------------------------------------------------------------------
NI=57,NL=269 |
confidence level |
No. of false alarms |
No. of hits |
In R console -> |
c |
f |
h |
----------------------- |
----------------------- |
----------------------------- |
------------- |
definitely present |
3 |
0 |
21 |
probably present |
2 |
7 |
4 |
questionable |
1 |
36 |
3 |
---------------------------------------------------------------------------------------------------
* false alarms = False Positives = FP
* hits = True Positives = TP
Note that in FROC data, the confidence level means present (deseased, positive) case only. Since each reader marks their suspicous location only and it generate the hits and false alarms for his confidenc level representing that lesion is present.
In the absent case, reader does not mark any locations and hence, the absent cofidence level does not relate this dataset.
Note that the first column of confidence level vector c
should not be specified. If specified, will be ignored , since it is created by c <-c(rep(C:1))
automatically in the program and it does not refer from user input data even if it is specified explicitly, where C
is the highest number of confidence levels.
So you should check the compatibility of your data and the program's generating new confidence level vector by
a table which can be displayed by the function viewdata()
.
Note that The format for the above example data must be made by the following forms:
dat <- list(
c=c(3,2,1), #Confidence level
h=c(21,4,3), #Number of hits for each confidence level
f=c(0,7,36), #Number of false alarms for each confidence level
NL=60, #Number of lesions
NI=30, #Number of images
C=3) #Number of confidence level
This R object dat
can be passed to the function fit_Bayesian_FROC()
as the following manner fit_Bayesian_FROC(dat)
.