Last data update: 2014.03.03
R: Obuchowski-Rockette analysis with Hillis improvements
ORHAnalysis R Documentation
Obuchowski-Rockette analysis with Hillis improvements
Description
Performs Obuchowski-Rockette analysis with Hillis improvements for the specified dataset.
Usage
ORHAnalysis(dataset, fom = "wJAFROC", alpha = 0.05,
covEstMethod = "Jackknife", nBoots = 200, option = "ALL")
Arguments
dataset
See DBMHAnalysis.
fom
See DBMHAnalysis.
alpha
See DBMHAnalysis.
covEstMethod
The method used to estimate the covariance matrix; can be "Jackknife"
, "Bootstrap"
or "DeLong"
,
the last assumes the Wilcoxon FOM is chosen, otherwise an error will result.
nBoots
The number of bootstraps (default is 200), used only if the "Bootstrap"
method is used to estimate
the covariance matrix.
option
See DBMHAnalysis.
Value
The return value is a list with following elements:
fomArray
Figures of merit array. See the return of FigureOfMerit.
msT
Mean square of the figure of merit corresponding to the the treatment effect.
msTR
Mean square of the figure of merit corresponding to the the treatment-reader effect.
varComp
Obuchowski-Rockette variance component and covariance estimates.
fRRRC
See DBMHAnalysis.
ddfRRRC
See DBMHAnalysis.
pRRRC
See DBMHAnalysis.
ciDiffTrtRRRC
See DBMHAnalysis.
ciAvgRdrEachTrtRRRC
See DBMHAnalysis.
fFRRC
See DBMHAnalysis.
ddfFRRC
See DBMHAnalysis.
pFRRC
See DBMHAnalysis.
ciDiffTrtFRRC
See DBMHAnalysis.
ciAvgRdrEachTrtFRRC
See DBMHAnalysis.
ciDiffTrtEachRdr
See DBMHAnalysis.
varCovEachRdr
Obuchowski-Rockette Variance and Cov1 estimates for each reader.
fRRFC
See DBMHAnalysis.
ddfRRFC
See DBMHAnalysis.
pRRFC
See DBMHAnalysis.
ciDiffTrtRRFC
See DBMHAnalysis.
ciAvgRdrEachTrtRRFC
See DBMHAnalysis.
Examples
retOR <- ORHAnalysis(rocData, fom = "Wilcoxon")
print(retOR)
Results
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(RJafroc)
Loading required package: tools
Loading required package: xlsx
Loading required package: rJava
Loading required package: xlsxjars
Loading required package: ggplot2
Loading required package: stringr
Loading required package: shiny
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/RJafroc/ORHAnalysis.Rd_%03d_medium.png", width=480, height=480)
> ### Name: ORHAnalysis
> ### Title: Obuchowski-Rockette analysis with Hillis improvements
> ### Aliases: ORHAnalysis
>
> ### ** Examples
>
> retOR <- ORHAnalysis(rocData, fom = "Wilcoxon")
> print(retOR)
$fomArray
Rdr - 0 Rdr - 1 Rdr - 2 Rdr - 3 Rdr - 4
Trt - 0 0.9196457 0.8587762 0.9038647 0.9731079 0.8297907
Trt - 1 0.9478261 0.9053140 0.9217391 0.9993559 0.9299517
$msT
[1] 0.004796171
$msTR
[1] 0.0005510306
$varComp
varCov
Var(R) 0.0015349993
Var(T*R) 0.0002004025
COV1 0.0003466137
COV2 0.0003440748
COV3 0.0002390284
Var(Error) 0.0008022883
$fRRRC
[1] 4.456319
$ddfRRRC
[1] 15.25967
$pRRRC
[1] 0.05166569
$ciDiffTrtRRRC
Treatment Estimate StdErr DF t Pr > t CI Lower
1 0 - 1 -0.04380032 0.02074862 15.25967 -2.110999 0.05166569 -0.0879595
CI Upper
1 0.0003588544
$ciAvgRdrEachTrtRRRC
Treatment Area StdErr DF CI Lower CI Upper
1 0 0.8970370 0.03317360 12.74465 0.8252236 0.9688505
2 1 0.9408374 0.02156637 12.71019 0.8941378 0.9875369
$fFRRC
[1] 5.475953
$ddfFRRC
[1] Inf
$pFRRC
[1] 0.01927984
$ciDiffTrtFRRC
Treatment Estimate StdErr DF t Pr > t CI Lower
1 0 - 1 -0.04380032 0.01871748 Inf -2.340075 0.01927984 -0.08048591
CI Upper
1 -0.00711473
$ciAvgRdrEachTrtFRRC
Treatment Area StdErr DF CI Lower CI Upper
1 0 0.8970370 0.02428971 Inf 0.8494301 0.9446440
2 1 0.9408374 0.01677632 Inf 0.9079564 0.9737183
$ciDiffTrtEachRdr
Reader Treatment Estimate StdErr DF t Pr > t CI Lower
1 0 0 - 1 -0.02818035 0.02551213 Inf -1.1045864 0.26933885 -0.07818322
2 1 0 - 1 -0.04653784 0.02630183 Inf -1.7693768 0.07683102 -0.09808848
3 2 0 - 1 -0.01787440 0.03120965 Inf -0.5727202 0.56683414 -0.07904418
4 3 0 - 1 -0.02624799 0.01729129 Inf -1.5179891 0.12901715 -0.06013829
5 4 0 - 1 -0.10016103 0.04405746 Inf -2.2734182 0.02300099 -0.18651207
CI Upper
1 0.021822507
2 0.005012792
3 0.043295388
4 0.007642316
5 -0.013809995
$varCovEachRdr
Reader Var Cov1
1 0 0.0006989006 3.734661e-04
2 1 0.0011060528 7.601598e-04
3 2 0.0008423434 3.553224e-04
4 3 0.0001505777 1.083399e-06
5 4 0.0012135668 2.430368e-04
$fRRFC
[1] 8.704
$ddfRRFC
[1] 4
$pRRFC
[1] 0.04195875
$ciDiffTrtRRFC
Treatment Estimate StdErr DF t Pr > t CI Lower
1 0 - 1 -0.04380032 0.01484629 4 -2.950254 0.04195875 -0.08502022
CI Upper
1 -0.00258042
$ciAvgRdrEachTrtRRFC
Treatment Area StdErr DF CI Lower CI Upper
1 0 0.8970370 0.02482994 4 0.8280981 0.9659760
2 1 0.9408374 0.01615303 4 0.8959894 0.9856854
>
>
>
>
>
> dev.off()
null device
1
>