DVHmetrics-package: Analyze Dose-Volume Histograms and Check Constraints

Description Details Author(s) References Examples

Description

Functionality for analyzing dose-volume histograms (DVH) in radiation oncology: Read DVH text files, calculate DVH metrics, gEUD, BED, EQD2, NTCP, TCP, show DVH diagrams, check and visualize quality assurance constraints for the DVH. Includes web-based graphical user interface.

Details

Package: DVHmetrics
Type: Package
Version: 0.3.7
Date: 2016-11-22
License: GPL (>= 2)

Author(s)

Daniel Wollschlaeger and Heiko Karle

Maintainer: Daniel Wollschlaeger <[email protected]>

References

For a solution that also reads files in DICOM-RT format, see the RadOnc package: https://CRAN.R-project.org/package=RadOnc.

Examples

1
2
showDVH(dataMZ[[1]])
checkConstraint(dataMZ, "D1CC < 10Gy")

Example output

   patID structure     constraint  observed compliance      deltaV    deltaVpc
1   P123    AMYOCL D1CC_GY < 10GY 25.606184      FALSE  58.7000000  5870.00000
2   P123    AMYOCR D1CC_GY < 10GY 33.427681      FALSE  55.2289870  5522.89870
3   P123   AOVALVE D1CC_GY < 10GY 25.484937      FALSE  11.5000000  1150.00000
4   P123    AVNODE D1CC_GY < 10GY 24.761131      FALSE  36.6000000  3660.00000
5   P123     HEART D1CC_GY < 10GY 33.511656      FALSE 579.1902905 57919.02905
6   P123   MYOCARD D1CC_GY < 10GY 33.430537      FALSE 299.4555356 29945.55356
7   P123 PULMVALVE D1CC_GY < 10GY 23.690379      FALSE   2.8000000   280.00000
8   P234    AMYOCL D1CC_GY < 10GY  1.716810       TRUE  -1.0000000  -100.00000
9   P234    AMYOCR D1CC_GY < 10GY  3.928471       TRUE  -1.0000000  -100.00000
10  P234   AOVALVE D1CC_GY < 10GY  2.709391       TRUE  -1.0000000  -100.00000
11  P234    AVNODE D1CC_GY < 10GY  2.025139       TRUE  -1.0000000  -100.00000
12  P234     HEART D1CC_GY < 10GY  4.059292       TRUE  -1.0000000  -100.00000
13  P234   MYOCARD D1CC_GY < 10GY  3.929008       TRUE  -1.0000000  -100.00000
14  P234 PULMVALVE D1CC_GY < 10GY  1.370736       TRUE  -1.0000000  -100.00000
15  P345    AMYOCL D1CC_GY < 10GY  6.135375       TRUE  -1.0000000  -100.00000
16  P345    AMYOCR D1CC_GY < 10GY 43.841965      FALSE  46.6646599  4666.46599
17  P345   AOVALVE D1CC_GY < 10GY 23.429635      FALSE   3.5425994   354.25994
18  P345    AVNODE D1CC_GY < 10GY 23.158001      FALSE   3.7795000   377.95000
19  P345     HEART D1CC_GY < 10GY 44.093171      FALSE 180.8211505 18082.11505
20  P345   MYOCARD D1CC_GY < 10GY 43.840794      FALSE  89.1757248  8917.57248
21  P345 PULMVALVE D1CC_GY < 10GY 14.356940      FALSE   0.7857872    78.57872
      deltaD  deltaDpc     dstMin  dstMinRel   ptMinD    ptMinV
1  15.606184 156.06184 14.1106499  460.90311 23.33669  5.609031
2  23.427681 234.27681 20.6222809  969.73434 28.20000 10.697343
3  15.484937 154.84937 11.5000000 1150.00000 10.00000 12.500000
4  14.761131 147.61131 12.2123321  416.28357 21.48093  5.162836
5  23.511656 235.11656 21.6479338  625.38304 30.72493  7.253830
6  23.430537 234.30537 21.3561083  717.05647 30.11632  8.170565
7  13.690379 136.90379  2.8000000  280.00000 10.00000  3.800000
8  -8.283190 -82.83190  1.0000000  100.00000 10.00000  0.000000
9  -6.071529 -60.71529  1.0000000  100.00000 10.00000  0.000000
10 -7.290609 -72.90609  1.0000000  100.00000 10.00000  0.000000
11 -7.974861 -79.74861  1.0000000  100.00000 10.00000  0.000000
12 -5.940708 -59.40708  1.0000000  100.00000 10.00000  0.000000
13 -6.070992 -60.70992  1.0000000  100.00000 10.00000  0.000000
14 -8.629264 -86.29264  1.0000000  100.00000 10.00000  0.000000
15 -3.864625 -38.64625  1.0000000  100.00000 10.00000  0.000000
16 33.841965 338.41965 30.1154429 1120.57472 37.95302 12.205747
17 13.429635 134.29635  3.4463258  335.58309 10.78458  4.355831
18 13.158001 131.58001  3.6544427  353.90072 10.91125  4.539007
19 34.093171 340.93171 30.9597340  992.92160 39.32432 10.929216
20 33.840794 338.40794 30.1302775 1093.10739 38.07749 11.931074
21  4.356940  43.56940  0.7461268   71.43809 10.21533  1.714381
Warning messages:
1: In getVolume(pm$valRefNum, type = "absolute", unitRef = pm$unitRef,  :
  max dose is less than requested dose
2: In getVolume(pm$valRefNum, type = "absolute", unitRef = pm$unitRef,  :
  max dose is less than requested dose
3: In getVolume(pm$valRefNum, type = "absolute", unitRef = pm$unitRef,  :
  max dose is less than requested dose
4: In getVolume(pm$valRefNum, type = "absolute", unitRef = pm$unitRef,  :
  max dose is less than requested dose
5: In getVolume(pm$valRefNum, type = "absolute", unitRef = pm$unitRef,  :
  max dose is less than requested dose
6: In getVolume(pm$valRefNum, type = "absolute", unitRef = pm$unitRef,  :
  max dose is less than requested dose
7: In getVolume(pm$valRefNum, type = "absolute", unitRef = pm$unitRef,  :
  max dose is less than requested dose
8: In getVolume(pm$valRefNum, type = "absolute", unitRef = pm$unitRef,  :
  max dose is less than requested dose

DVHmetrics documentation built on Nov. 22, 2017, 5:04 p.m.